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Charged Hadron Spectra in Au+Au Collisions at ^/s^ = 130 GeV 
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December 2001 

The collision of high energy heavy ions is the most promising laboratory for the study of 
nuclear matter at high energy density and for creation of the Quark-Gluon Plasma. A new 
era in this field began with the operation and first collisions of Au nuclei in the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory during 2000. This work 
concentrates on measurement of global hadronic observables in Au+Au interactions at a 
centre-of-mass energy of ^/s^ = 130 GeV, which mainly address conditions in the final state 
of the collision. The minimum bias multiplicity distribution, the transverse momentum {p±), 
and pseudorapidity (ry) distributions for charged hadrons (/i^, h^) are presented. Results 
on identified vr" transverse mass (m^) and rapidity (y) distributions are also discussed. The 
data were taken with the STAR detector with emphasis on particles near mid-rapidity. 

We find that the multiplicity density at mid-rapidity for the 5% most central interactions 
is dNf^- / dr]\rj=Q = 280 ± Igtat i 20syst, an increase per participant of 38% relative to pp 
collisions at similar energy. The mean transverse momentum is 0.508 it 0.012 GeV/c and 
is larger than in Pb + Pb collisions at lower energies. The scaling of the h~ yield per 
participant nucleon pair, obtained via a ratio of Au + Au to pp pj_ distributions, is a strong 
function of p±. The pseudorapidity distribution is almost constant within |r/| < 1. The vr" 
rapidity distribution is also flat around mid-rapidity in the region \y\ < 0.8, with a yield of 
pions for central collisions of dN^- / dy\y=o = 287 ± Igtat ± 21syst. However, the slope of the 
m± distributions is not the same for different rapidity bins, suggesting that boost-invariance 
is not fully achieved in the collisions. 
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Chapter 1 



Introduction 



Based on statistical Quantum Chromodynamics, our current expectation is that strongly 
interacting matter at extreme energy density (e ^ 1 GeV/fm^) and temperature (T > 
170 MeV) is found in a state where hadrons no longer exist as discrete entities The 
relevant degrees of freedom for such a system are those of the underlying partons, and we 
label this state the Quark-Gluon Plasma (QGP). ||3|-^ 

Such a state of matter is believed to be the one in which the early universe existed in a 
time scale ~ 10^^ s after the Big Bang and is also predicted to exist in the interior of 
neutron stars P,p^. 

Ultra-relativistic heavy ion collisions are the most promising tool for the creation of a 
QGP in the laboratory. The appearance and study of predicted signatures in such colli- 
sions has been the subject of intense theoretical and experimental work for more than two 
decades (for recent reviews one can refer to the proceedings of the 'Quark Matter' confer- 
ences |11-25|). In the Alternating Gradient Synchrotron (AGS) at the Brookhaven National 
Laboratory (BNL) collisions have been produced using beams from Si to Au with energies 
of 11-14 GeV per nucleon. The corresponding centre-of-mass energy for Au -|- Au collisions 
at the AGS is ^/s^ = 5 GeV. We use ^/s^ to denote the centre-of-mass energy per nu- 
cleon pair. The Super Proton Synchrotron (SPS) at the Conseil Europeen de la Recherche 
Nucleaire (CERN) has yielded results from O, S and Pb beams with energies of 60-200 GeV 
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per nucleon, this is a centre-of-mass energy of y^s^ = 17.2 GeV for Pb + Pb collisions. 
There has been considerable excitement in the field during the last year with the com- 



missioning of the Relativistic Heavy Ion Collider (RHIC) at BNL |Q. Dedicated experi- 
ments began taking data in the summer of 2000 with the highest colliding energy heavy ion 
beams available. RHIC ran during this period with Au beams at a center of mass energy of 
^/s^ = 130 GeV. Studying collisions at different energies helps to map the phase diagram 
of nuclear matter, where RHIC is expected to probe the region of high temperature and 
near-zero net-baryon density, a regime accessible by QCD lattice simulations. 

For all heavy ion experiments, global event observables have played an important role. 
It remains true that to understand specific plasma signatures one must first understand the 
global character of the reaction dynamics. The momentum distribution of the bulk of the 
particles measured in the detectors provide evidence mainly of the final phase of the system 
formed in the collision, also called the freeze-out phase. In order to obtain information about 
the early stages where plasma formation is expected to occur, one must use indirect methods 
such as hard probes which must be gauged to reference observables related to the global 
character of the reaction. The information obtained in the final state in the form of particle 
spectra is one of the main sources of global event information. In addition, global final 
state observables help provide limits on the possible evolution of the system at earlier times. 
All experiments at RHIC and elsewhere have some capability to measure global observables 
in order to characterize and classify events, compare results with other experiments, and 
perform systematic studies. Measuring the final state particle spectra is therefore a basic 
requirement for the study of the collision dynamics. 

The measurements of global observables for the first year of RHIC collisions are an 
exciting first topic of study in an energy regime where perturbative phenomena are expected 
to dominate. The contribution to the total charged particle multiplicity coming from hard 



processes (jets and mini-jets) at RHIC is estimated to be between 30 — 50% [|27H30||. Since 
particle production is a dominant feature of the collision, one of the first observables to 
study is the multiplicity of charged particles for each event. This quantity and related global 



observables such as the transverse energy {E±) and the energy deposition in the very forward 
region as measured typically by zero degree calorimeters (ZDCs) is related to the impact 
parameter of the collision. The more central the collision is, the more nucleons from both 
projectile nuclei will participate in the interaction, and hence the more secondary hadrons 
will be produced (large multiplicity and total E±). A geometrical correlation between the 
impact parameter b, E± and the energy reaching the zero degree calorimeter can be given in 
terms of the number of "wounded" nucleons ||3l|]. Central (high multiplicity) events have a 
special interest as it is for these events that the largest fraction of the incoming energy will be 
redistributed to new degrees of freedom. The next key observable is the distribution of the 
particles in momentum space. Rapidity and transverse momentum distributions allow one 
to address properties of the reaction dynamics such as the extent to which the nucleons are 
slowed down in the collision (stopping), the approach of the system to thermal equilibrium, 
and the shape of the emitting source of particles. 

The Solenoidal Tracker At RHIC (STAR) is one of the 4 experiments that partake in 
the RHIC program. The analyses presented here are based on data taken by STAR during 
the 2000 summer run. The pith of the experiment is a large acceptance cylindrical Time 
Projection Chamber (TPC) placed in a uniform solenoidal magnetic field for momentum 
determination. The TPC provides charged particle tracking in the mid-rapidity region with 
full azimuthal coverage and particle identification for low momentum particles. The STAR 
detector is thus ideally suited for the study of hadronic observables and will therefore focus 
on such measurements early on, although the STAR physics program includes other topics 
in addition to QGP physics. Since hadrons are the most copiously produced particles in 
the collision and n mesons comprise ^ 80% of the total hadron population, this work con- 
centrates on the study of charged hadron (/i~ and /i~^) and identified vr" meson production 
and momentum spectra for Au + Au collisions at ^/s^ = 130 GeV. The minimum bias 
h~ multiplicity distribution is presented. For different selections of event centralities, we 
present the pseudorapidity [rj) and transverse momentum {p±) distributions for h~ . We 
compare the p± distributions to those expected in similar energy pp collisions as reference. 
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We also discuss rapidity (y) and transverse mass {m±) distributions for identified tt~ . 

This work is organized in the following manner. It starts with a discussion on general 
aspects of Quark-Gluon Plasma physics and a review of hadronic particle production at the 
AGS and SPS (chapters |3| and followed by a description of the STAR detector (chapter 

An overview of the tracking strategy and detector calibration procedures is then given 
(chapter After a discussion of the detector simulation (chapter |^ and a layout of the 
analysis technique and corrections (chapter |8|), we present the results on hadron production 
and identified pion spectra and discuss their implications regarding the dynamics of the 
collision features at RHIC (chapters ^ and[l^). Finally, we present our conclusions (chapter 
0). 
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Chapter 2 

The Physics of the Quark-Gluon 
Plasma 



We give now an overview of the high temperature phase of QCD and the Quark-Gluon 



Plasma. For recent reviews on the subject, see e.g. |32-34], review articles in the proceedings 



of the 'Quark Matter' conferences [35,36] and the growing collection of books [p|,p|,|^,|37, 38 



2.1 Deconfinement in QCD 

At the fundamental level, strongly interacting matter is described by interactions of quarks 
through the exchange of gluons. The theory that describes these interactions, quantum 
chromo dynamics (QCD), has the remarkable properties that at large distances or small 
momenta g, the effective coupling constant as{q^) is large, and it decreases logarithmically 
at short distances or large momenta. This behaviour can be seen from perturbative QCD. 
The QCD Lagrangian is given by 

Cqcd = ih^{d^ - igA^)i^ - mi^ij - ^F^^Fr (2.1) 

where the colour potential is a 3 x 3 matrix (indicated by the circumflex symbol) and 
can be represented by a linear combination of the 8 Gell-Mann matrices: 
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(2.2) 



a=l 



This potential is introduced to make the Lagrangian invariant under local (space-dependent) 
rotations of the 3 colour components of the quark wavefunction ip. The eight-component 
field strength tensor expressed as 



where fabc are the antisymmetric structure constants for the Lie group SU(3). The product 
F^j^Fa'^ also remains invariant under a local colour gauge transformation. 

The perturbative quantization proceeds in QCD in a similar manner as in Quantum 
Electrodynamics. The quadratic terms in the lagrangian define free quark and gluon fields, 
described by propagators with the same form as those in QED for electrons and photons. 
The free propagators Do{q'^) are proportional to and if this were the only ingredient 

would lead to a colour force which would fall off like l/r. The main difference comes from 
the coupling of the gluon field to itself. This modifies the true gluon propagator, because one 
has to evaluate contributions arising from terms in the perturbation series corresponding to 
the virtual creation of a pair of coloured particles from the vacuum. We call this the 'vacuum 
polarization' function n(g^), and it turns out to be proportional to — ln(— g^). 

The higher order diagrams, in which the gluon interacts consecutively once, twice, three 
times with the vacuum polarization, etc. can be summed into a geometric series yielding 
the full propagator 





+ QfabcAnAI , 



(2.3) 



Doiq^) + i^Doiq^Mq^)Doiq^) + i^DoUDoUDo + ... 



l-n(g2)Do(52) 




\00 — ^iSp )(Xs 

127rln(-qV//^) 



(2.4) 
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where /i is a reference point introduced by renormalization, and Np counts the number of 



flavours with mass below So after renormalization, the second factor in Eq. UA acts 

as a momentum-dependent modification of the strong coupling constant. Combining it with 
as we obtain the "running" coupling constant 



"^^'^ ' " (ll-27V^/3)ln(-g2/A2) ^ M^JW) 
where A is a (dimensional) parameter introduced also by the renormalization process. We 
see that the running coupling exhibits a pole; in this approximation it is at = — A^. More 
sophisticated expressions for the gluon propagator indicate that the pole is really at = 0, 
and that as{q^) should behave as 1/q^ in the limit q^ —>■ 0. 

Converted into coordinate space, this means that as{r) grows like for large distances. 



corresponding to a linearly rising potential. The logarithm in Eq. |2.5| causes a gradual 
decrease of the coupling strength between colour charges at large momenta or small distances. 
It is this property that is known as ^asymptotic freedom'. This behaviour of the running 
coupling at large distances results in the confinement of quarks {i.e. isolated quarks are not 
observed in nature). 

It is important to note that the behaviour of the strong coupling constant outlined above 
is derived for interactions in vacuo. The usual point of comparison for measurements and 
calculations of the strong coupling constant is at the mass of the boson, Mz = 91 GeV/c^, 
where the world average is as{Mz) = 0.1172 it 0.0045 [^]. The typical initial momentum 
transfers at even the highest RHIC energies are significantly lower than Mz, so we are 
really talking about a running of the coupling with temperature. One must instead focus 
on obtaining an expression for the coupling constant using expressions for the propagators 
which have corrections due to the presence of a coloured medium Q. Using this modified 
propagator we again sum all diagrams with successive medium polarization functions n(q, T) 
to obtain an effective running coupling constant. Leaving aside the contribution of quark 
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loop diagrams, this yields 



^ " 47r{l -U{q,T)/q^) " llln(-g2/A2) - 48G((7/r) ^^'^^ 



where 



oo 



G(C) = / zdz /^^^ and 
7o - 1 

f{z) = 1/2Z + 1/8^3) In -1 + 1/2^2. (2.7) 

In the g — > limit, the polarization function remains finite, which means that the propagator 
eff'ectively contains a mass term, D{q, T) = 1/ {q^ — Il{q, T)), in a manner analogous to Debye 
screening in an electrolytic medium. This leads to the property that a test colour charge will 
cause a polarization of the charges in the coloured medium in the same way as electric charges 
in an electrolyte. In addition, an important property that follows from the temperature 
dependence of Eq. ^jo] is that as the temperature increases in QCD, the coupling becomes 
weak, falling logarithmically with increasing temperature. As a consequence, nuclear matter 
at very high temperature should not exhibit confinement. 

Another important property that arises from the study of the QCD Lagrangian is that 
of chiral symmetry breaking {i.e. the quarks confined in hadrons do not appear as nearly 
massless constituents, but instead possess a mass of a few hundred MeV that is generated 
dynamically). The expectation value (V'V') ~ —(235 MeV)^, usually called the quark con- 
densate, describes the density of qq pairs found in the QCD vacuum, and the fact that it is 
non- vanishing is directly related to chiral-symmetry breaking. In the limit of zero current 
quark mass, the quark condensate vanishes at high temperature, i.e. chiral symmetry is 
restored. It is this phase of QCD which exhibits neither confinement nor chiral-symmetry 
breaking that we entitle the Quark-Gluon Plasma. 

Since the quark condensate is zero in the high temperature phase and non-zero at low 
temperature, it therefore acts as an order parameter. This behaviour leads to the expectation 
that the change between the low-temperature and high-temperature phases should exhibit 
a discontinuity. For 2 (3) massless quark flavours, universality arguments predict a second- 



(first-) order phase transition |40,^]. The question of the order of the transition, or if there 
is a phase transition as opposed to a rapid cross-over, for QCD with the real values for the 
u, d and s quark masses is still the subject of current investigations |42,43|. 



2.1.1 Phase Diagram 
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Figure 2.1: Phase diagram of hadronic matter, showing hadron gas and quark-gluon plasma 
regions. The temperature T and the baryochemical potential /x^ data are derived from 
particle yield ratios. The solid curve through the data points represents the chemical freeze- 
out of hadronic-matter. 



One of the first pictures that can serve as a guide to the behaviour of QCD can be 
obtained from the simplified 'MIT bag model' Q]. We can identify two regions where we 
expect deconfinement to occur: when compression causes the hadrons to overlap signifi- 
cantly, reaching densities 3 to 5 times higher than those of ordinary nuclear matter (0.17 
fm~'^), or when the temperature of the medium exceeds some critical threshold. Figure 



2A\ shows the usual representation of the — fiB^ phase diagram, where /i_B is the Baryon 
chemical potential. It is normally represented with a continuous curve connecting the high 
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temperature transition region at /x^ = with the high baryon density region at T = 0. The 
figure shows the regions probed by the different beam energies from the SIS to RHIC. At 
high densities and near-zero temperature we expect a deconfined phase. This deconfined 
high-density phase is predicted to exist for example in the interior of neutron stars. The 
near-zero Baryon chemical potential and high temperature region is the one probed by the 
highest energy RHIC collisions. The region where /i^ = is believed to be the one in which 
the early universe existed, and is also accessible to numerical simulations of QCD on a lat- 
tice. The data points in the figure are from an analysis of particle ratios which is customarily 
employed to evaluate the degree of chemical equilibration observed in the final state. To 
roughly illustrate the region that is probed by RHIC in the diagram, preliminary analysis 



from mid-rapidity particle ratios at RHIC [45| indicate a value of /x^ ~ 0.050 GeV and of 



the chemical freeze-out temperature of T^h ~ 170 MeV. We discuss chemical and kinetic 
freeze-out in Chapter ^. It should be noted however, that we do not really know where the 
transition curve really lies, or even if there is a transition as opposed to a rapid cross-over 
at finite baryon number density. First principle numerical calculations of QCD, discussed 
in the next section, have provided us with guidance, but most results so far pertain to the 
/i_B = region. 

2.2 Lattice QCD Results 

From sophisticated numerical simulations of QCD on the lattice, we have gained much 
insight into the structure of the QCD at high temperature. 

The approach is based on the Feynman path integral: the aim is to calculate the action 
on the lattice and use it to evaluate expectation values of different observables. The first 
results used pure SU(3) gauge theory, sometimes called pure glue theory. A problem arises 
when one introduces the (discrete) fermion fields in a straightforward way, one obtains a 
'doubling' of the flavour spectrum in the quark sector. We get 2*^ = 16 copies of each quark 
species (in d = 4 dimensions). Different approaches were then introduced to incorporate 
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the fermion degrees of freedom in a way that reduces the doubhng problem. The approach 
proposed by Wilson |^] solves the flavour doubling essentially by giving the doublers a 
mass proportional to 1/a, where a is the lattice spacing, so they go away in the continuum 
limit. The disadvantage is that the mass introduces a chiral-symmetry breaking piece into 
the action. Another action known as the Kogut-Susskind E^, or also as the 'staggered' 



fermion action |4^, ig], has the advantage of preserving a part of chiral symmetry. However, 
it does not completely eliminate the flavour problem. The number of doublers is reduced to 
4. This leads for example to having 4^ — 1 = 15 pions instead of the usual 2^ — 1 = 3. In 
addition, the Wilson action yields results that in general need corrections of 0{a), where a 
is the lattice spacing, whereas the staggered fermion action is fine up to 0{a^). 

In the last decade, there have been significant developments that provide a solution to the 
doubling problem. The "domain wall" approach relies on introducing a fifth (fictitious) 
dimension such that the chiral zero modes live on 4D surfaces. If is the number of 
lattice spacings along the fifth dimension, in the limit A^s — > oo the left- and right-handed 
fields live in surfaces in opposite ends of the 5D lattice and do not mix, so we have exact 
chiral symmetry with no doublers. An alternative approach has been presented where chiral 
symmetry is preserved in the lattice if the lattice Dirac operator is of a certain form (see 



e.g. Ref. ||5l|). Although several constructions for the Dirac "overlap" operator exist, they 



all satisfy an identity, originally due to Ginsparg and Wilson ||52||, which guarantees having 
exact global chiral symmetries directly on the lattice. For a recent review of the overlap 
approach, see e.g. Refs. (53,54|. 



2.2.1 Critical Temperature and Energy Density 

Over recent years, thermodynamic calculations on the lattice have steadily been improved. 
This is partly due to the much improved computer resources, however equally important 
has been (and will continue to be) the development of improved discretization schemes, i.e. 
improved actions. 

Early calculations of the QCD transition temperature performed with standard Wilson 
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fermions |^,55] and staggered fermion actions |49| led to significant discrepancies of the 
results. These differences were greatly reduced based on improved Wilson fermions (Clover 
action) |19,^, 57|, as well as improved staggered fermions ||2[ 58|, and domain wall 



approaches. We shall not discuss these here, but only take the current results from improved 
staggered fermion actions for the following discussion. 

Typically, one can look at the dependence of the energy density e and pressure versus 
temperature where one expects a change due to the increase in the number of degrees of 
freedom {d.o.f ). In order to illustrate this point, we can obtain some semi-quantitative 
insight into the number of degrees of freedom and the energy density using the following 
simplified scenario. A massless non-interacting hadron gas is made up basically of pions, of 
which we have 3 types (vr^, vr^, vr") neglecting the resonances. From an ideal relativistic 
Bose gas at Temperature T we obtain the energy density 



d^p p AttT^ dx ^ 7r2 4 

-dx = —T , {2.i 



(27r)3eP/^-l (27r)3 - 1 30' 

where we have rescaled the momentum as x = p/T. This is the usual Stefan-Boltzmann 
relation (e = aT^). The energy density for the hadron gas is therefore simply 

e^G = 3 eg = 3 — . (2.9) 

For a QGP consisting of 2 massless quark flavours (u and d) at vanishing net baryon 
density (p = 0) we must sum the quark and the gluon contributions to the energy density. 



The gluon contribution is given by Eq. 2^ times the 16 gluonic d.o.f. (2spin x (3^ — l)colourJ 



The quark contribution is obtained from a similar integral for a Fermi gas 



(27r)3 eP/^ + 1 (27r)3 Jq + I 8 30' 

Multiplying by the number of (anti)quark d.o.f. (2gpin x 3coiour x Np) we obtain the energy 
density 
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eQGP = + 6iVF(6g + e,-) = 16^r^ + 12iV^^^r^ = (16 + '^Nf)'^T^ . (2.11) 



Figure shows recent results for the energy density e as a function of temperature 
in lattice QCD simulations with Np = 0, 2 and 3 light quarks as well as two light and 
a heavier (strange) quark (2+1 flavour QCD). The pressure p is shown in Figure ^ a for 
QCD with different number of flavours as well as for the pure SU(3) gauge theory. The 
curves clearly reflect the strong change in the number of degrees of freedom when going 
through the transition. In the high temperature limit (T > 1.5 Tc), we expect both e and p 



to asymptotically approach the Stefan-Boltzmann free gas limit Eq. |2.11| , indicated by the 
arrows in the figures. From the figure, it is evident that even at ATc the Stefan-Boltzmann 
limit is not reached. This has been taken as indication of a significant amount of interactions 
among partons in the high temperature phase, with only a logarithmic approach to the free- 
gas behaviour. 




Figure 2.2: Energy density e obtained from a numerical evaluation of QCD on the lattice |42|. 
e is divided by to exhibit the sudden rise of the number of thermally excited degrees of 
freedom near the critical temperature Tc. Arrows show the ideal gas values as given by 
Eq. 
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In addition, the dependence of Tc on the number of partonic degrees of freedom is clearly 
visible in Figure ^.3| a. It is therefore striking that p/psB is almost flavour independent when 
plotted in units of T/Tc as shown in Figure |2.3| b. 




Figure 2.3: The Pressure p in lattice QCD [^]. Left panel: Curves with Np = 0,2 and 3 
light quarks as well as the "(2+1)" case. Right panel: T axis scaled by the value of Tc for 
each curve. 

The most recently reported results on the value of Tc from lattice calculations are found 
to be, for one particular choice of actions (improved staggered fermion action) |Q 



Nf = 2 : Tc = 173 ± 8MeV 

Nf = 3 : Tc = 154 ± 8MeV . (2.12) 

The results also suggest that the transition temperature in (2+l)-flavour QCD is close to 
that of 2-flavour QCD. 

The behaviour of the different actions (improved staggered, Wilson, Clover, etc.) cur- 
rently studied in lattice QCD should show agreement in the vicinity of the phase transition. 
In this regime, the correlation lengths become large and cut-off effects in the calculations 
become less important. One can therefore compare calculations made with different actions. 
In particular, the recent results for the energy density at Tc yield 
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ec ~ (6 ± 2)rc^ = 0.70 ± 0.23 GeV/fml (2.13) 



as shown in Figure by the vertical line at T = Tc. It is important to stress that these 
values refer to initial energy densities, when the medium exists at the early stages of QGP 
evolution, and must be translated into a final energy density that can be measured in 
an experiment using detected particles. Bjorken introduced a relation to address this 
question. One can experimentally estimate the energy density achieved in a nucleus-nucleus 
collision via the energy measured in the central rapidity region, dE±/dy, divided by the 
effective interaction volume, (Aj- • r). In this case, we estimate the energy by a product 
of the particle multiplicity {dN / dy\y=o) , times the mean transverse mass of the particles, 
((mx)). 

.^^'7° ■<"">. (2.4, 

In this expression, the factor At is the transverse area of overlap in the collision. The 
quantity r is less accurately defined. It is normally interpreted as the parton formation 
time, i.e. the time needed to pass from the initial hadronic environment to the partonic 
degrees of freedom. Usually, this time is taken as 1 fm/c at SPS in order to compare different 
experiments. However, there is currently no real consensus as to what is the appropriate 
formation time r to use at RHIC, although if anything there are arguments that it should 
be smaller than 1 fm/c at high energies because it should take less time to equilibrate 
the system. An estimate of the energy density using this relation should then be at best 
considered a lower limit (with respect to this formation time). 

For reference, in the case of STAR, this relation can be calculated for a given centrality 
selection, we choose the 5% most central collisions, using the total charged multiplicity and 
the transverse momentum distribution for charged particles (Sec. and |9]^ respectively). 
The transverse area At can be computed in a geometrical model that reproduces the mea- 
sured multiplicity distribution (Sec. |4.1.1] ). Equation 2.14 is a very accessible experimental 
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observable, but there are additional caveats associated with its interpretation of which the 
formation time is but one example. We discuss the applicability and additional uncertainties 
associated with equation 2.14 in Section 



2.2.2 Heavy Quark Effective Potential 

In the lattice approach, an interesting quantity one can calculate is the effective potential 
between two heavy quarks, for it is the example which most cleanly illustrates the modifica- 
tion of the behaviour of the coloured fields at high temperature. Prom the previous section, 
we delineated the behaviour of the strong coupling constant: at large momenta or small 
distances it should be small and at small momenta or large distances the coupling rises as 
(asymptotic freedom). This leads to an effective potential that rises linearly with r, 
the distance between the coloured constituents, and can be calculated in the lattice. 



Figure shows the results of a lattice calculation for the heavy quark potential as a 
function of distance. We see the potential is very weak at small distances and the expected 
linear rise at large distances. The slope of the linear rise is usually called the string constant, 
since it can be thought of as the tension of a spring. 

One can test the idea of deconfinement at high temperature in lattice calculations in 
a straightforward way. As we raise the temperature, one can study the behaviour of the 
effective potential. We expect that near the critical temperature, the linear piece of the 
potential should be modified and weakened. At high temperature, the energy cost to create 
light qq pairs from the vacuum is reduced (this is related to the vanishing of the quark 
condensate and the restoration of chiral symmetry in the T > Tc phase of QCD). These light 
qq pairs then act as screening colour charges around the heavy quark pair, thus weakening 
the potential. 



Figure shows the heavy quark potential for different temperatures in the region near 
Tc. It is clear that there are important modifications to the strength of the potential at 
high temperatures. We see that as the temperature increases, the strength systematically 
decreases. The potentials in Figure also do not show the steep rise as in the quenched 
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Figure 2.4: The heavy quark effective potential as a function of distance r in lattice QCD. 
The left panel shows a quenched {Nj = 0) calculation using Wilson loops, in physical units. 
The right panel shows a comparison between a quenched calculation (red) and one including 
quark loop effects (blue) showing nice agreement. Vertical bars near the horizontal axis show 
the scale of the lattice spacing a, and intervals in physical units 0.1 fm wide for reference. 
Figures are from Refs. ||6l|| (a) and ||6^(b). 

case. This is an expected phenomenon that had proved elusive in lattice calculations. As 
the quark-antiquark pair separate, we expect the formation of a tube of flux, or string, 
which should break in the presence of light quark-antiquark pairs. The results in the figure 
show that the string does in fact break in the confining phase at non-zero temperature. The 
weakening of the potential at high temperature has important consequences for heavy quark 
bound states. In particular, it is this behaviour that is the basis for the concept of 3/ip 
suppression in a deconfined medium []6^ . 

We next need to address how we can measure the properties of excited hadronic and 
partonic matter in the laboratory and how can we hope to see a signal of QGP formation. 
We therefore proceed to discuss the experimental side of QGP physics. 
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Figure 2.5: The heavy quark effective potential for different temperatures, taken from 
Ref. |42|. The linear rise of the potential is weakened as one approaches the critical temper- 
ature. The solid curves show the Cornell potential V{r) = —a/r + ar with a = 0.25 it 0.05, 
which is used to normalize the finite temperature free energies at the shortest distance 
available, r = l/(4r). 
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Chapter 3 

Experimental Search for the QGP 



The use of beams of heavy elements allows us to distribute the incoming energy over an 
extended region of space, large compared to the size of one nucleon, with the hope of 
creating conditions that are suitable for the formation of a QGP. One can identify the 



following phases, shown schematically in Fig. 3.1 




Figure 3.1: 1-D Space-time picture of the evolution of an ultra-relativistic nuclear collision, 
distinguishing 3 regions: (a) Pre-equilibrium state, (b) thermalized QGP, (c) hadronization 
stage, (see text). The z axis denotes the beam direction, and the dashed lines indicate the 
colliding nuclei near the light-cone region. 

i) The interpenetration of the nuclei with partonic interactions at high energy. This 
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stage features the creation of e.g. high-p_L jets, cc pairs or other products of high momentum 
transfer scattering processes on the parton level. In addition, large cross-section soft nucleon 
nucleon scattering between the two highly Lorentz-contracted nuclei help redistribute a 
fraction of the incoming kinetic energy into other degrees of freedom. The small cross- 
section hard processes are used as experimental probes for the hot and dense zone. 

ii) The interaction of the particles of the system, driving it towards chemical and thermal 
equilibrium. Partons materialize out of the highly excited QCD field. If QGP forms, the 
quark mean free path at energy densities e = 2 GeV/fm^ is A 0.5 fm, and individual 
parton-parton scattering is expected to play a role in thermalizing the system during this 
early stage as the nuclear dimensions are larger than A. The interactions also originate the 
development of collective flow. Rapid expansion, mainly along the longitudinal direction, 
lower the temperature of the system eventually reaching the cross over temperature Tc. 
Direct photon signals from the QGP are generated from collisions of charged particles during 
the expansion stage, although there are also photon signals from the hot hadron gas stage. 
Final formation of charmonium states from the initial cc pairs happens also during 

this stage. 

Hi) The hadronization stage and the 'freeze-out' of the final state particles is then reached 
at T < Tc, as the system cools down so that there is not enough energy in each collision 
to further change the different species' populations or ratios {chemical freeze-out, T ~ 
160 — 170 MeV). Eventually, the energy is small enough and the system dilute enough 
such that the interactions cease and the momentum spectra do not change further {kinetic 
freeze-out, T ~ 120 MeV). To obtain information of the different stages, experimentally 
one must start from the measurement of final state particles. Global observables are useful 
to determine the initial conditions such as centrality, initial volume, and possibly energy 
density. 

There are some observables that also provide information from the early stages. Signals 
such as direct photons and dilepton pairs that originate during early times are interesting 
since they are little disturbed by the hadronic final state. The caveat in studying such signals 
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is that they typically have a much smaller cross section compared to hadronic observables. 

There are also a class of hadronic observables that are thought to be sensitive to early 
times. In particular, the measurement of an azimuthal anisotropy in the emission of particles 
(with respect to the reaction plane, i.e. the plane formed by the beam direction and the 
direction of a vector connecting the center of the two colliding nuclei) is one example. In 
non-central collisions, an initial spatial anisotropy would result in pressure gradients which 
drive the emission of particles, producing a modulation in the azimuthal distribution of 
particles with respect to the reaction plane. This effect is typically measured by a Fourier 
analysis of the azimuthal particle distribution, where the 2nd Fourier component is called 
elliptic flow (in reference to the picture of particles in a fluid moving under the influence 



of the initial pressure gradient) and is denoted as V2. The STAR experiment |64| already 
measured the charged particle elliptic flow signal at low p± , and from these results we already 
see evidence of significant collective behaviour during the early stages of collisions at RHIC 
energies. 

3.1 QGP Signatures 

Many different effects have been proposed as possible ways to detect experimentally the 
formation of the QGP state of matter. These range from the earliest and most naive studies 
which involved plotting mean p± {{p_l)) as a function of particle multiplicity as a quick way 
to look into the structure of the T — e phase diagram |^], to the more current searches for 
modification of particle properties (enhancement, suppression or medium-induced changes 
in mass or width), and the statistical studies of charge (or other observable) fiuctuations 
event-by-event. While some probes give information primarily about 'surface' effects, like 
the hadrons' p± distributions discussed here, strangeness production and particle interfer- 
ometry which reveal final state information; some others deal with deeper 'volume' effects 
which are sensitive to early times after the collision. These include most 'hard' processes. 
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i.e. those which involve large momentum transfers, for example charm and beauty produc- 
tion {J/tp, D mesons, T), jets and high-p^ particle production. Usually, the attention is on 
'volume' probes that address the features of the QGP itself, i.e. those probes which provide 
more direct information from the hot and dense phase of the reaction and are not influenced 
much by the high number of hadrons produced in the collision (e.g. by hadronic rescatter- 
ing). This is the ideal case, but the experience has been that one must carefully study the 
modification of the proposed signals by conventional nuclear means (i.e. non-QGP). It is 
also important to emphasize that all of the different clues must be investigated as a function 
of the associated particle multiplicity or equivalent probe giving information on global char- 
acteristics. One would like to understand the onset of all of the observed signals in terms of 
the different handles at our disposal. Experimentally we can vary the collision centrality (se- 
lecting on multiplicity or transverse energy), we can collide different nuclear species (a more 
controlled variation of system size) and we can vary the centre-of-mass energy. The experi- 
mental data will in turn help to constrain the theoretical model parameters and input {e.g. 
equation of state, expansion dynamics and collective flow, size and lifetime of the system 
for hydrodynamical models; and initial parton densities, parton mean-free path and cross 
section, nuclear shadowing of initial parton distributions, and amount of parton energy loss 
in the plasma for perturbative approaches). For example, using hadron and electromagnetic 
spectra from several SPS experiments, attempts have been made to constrain the equation 
of state and initial conditions in a hydrodynamical approach |66|. 



In the following we will briefly mention some of the experimental signals that have been 
proposed to probe the system created in heavy ion collisions. For a recent review on QGP 
signatures, see Ref. p^] . 

3.1.1 Charmonium Suppression 



The arguments for charmonium suppression were laid out in Section p.2.2| . Essentially, the 
weakening of the heavy quark effective potential with increasing temperature, or alternately 
viewed as the Debye screening of free colour charges in a QGP, is responsible for the breakup 
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of charmonium states |63|. Excited states of the cc system, such as ^pf and Xc, are easier to 
dissociate due to their larger radii, and are expected to dissolve just above Tc. The smaller 
J/ip becomes unbound at a higher temperature, T > 1.2Tc. Similar arguments apply to the 
dissociation of the heavier bb bound states, but they require much shorter screening lengths 
to dissolve, i.e. greater temperature and energy densities ||6^. The T state may dissolve 
above a temperature of ~ 2.5Tc, while the T' could also dissolve near Tc. 

"Anomalous" J ftp suppression has been reported by the NA50 collaboration for central 



Pb + Pb collisions at SPS which has been heralded as evidence for QGP formation ||68-70|. 
There are nuclear effects, such as the breakup of the J/ip hy hadronic comovers, which also 



suppress the measured cross section in nucleus-nucleus collisions |71|. This issue is still one 
of intense debate, with many journal publications devoted to the topic. There are already 
some model predictions on 3/ip suppression at RHIC and LHC energies based on normal 
nuclear effects, such as p± broadening and nuclear shadowing of the parton distribution 



functions |72|. Such an analysis lead for example to anti- shadowing at RHIC energies, i.e. 
the modifications introduced by cold nuclear matter enhance quarkonium production at 
RHIC (at LHC the opposite occurs). It has become clear that several mechanisms need to 
be considered in order to fully comprehend quarkonium production and suppression: the 
underlying physics description at different energies (hadronic at AGS and SPS, partonic 
at collider energies), nuclear effects such as shadowing, interaction with comovers, J/^ 
production due to the decay of higher mass resonances, and the signals from pA data. For the 
most recent developments, one can see for example the talks in the dedicated J/V' session at 



Quark Matter '01 |25|, where suppression from the comover mechanism |71|, a model based 



on E_\_ fluctuations for the most central events to explain the NA50 Pb + Pb data |73|, 



and even enhancement oi i/ijj production in deconfined matter at RHIC energies |Q were 
discussed. 
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3.1.2 Jet Quenching 

The colour structure of QCD matter can be probed by its effect on the propagation of a 
fast parton. The mechanisms are similar to those responsible for the electromagnetic energy 
loss of a fast charged particle in matter: energy may be lost either by the excitation of 
the penetrated medium or by radiation . The QCD analog of this effect indicates that the 
stopping power of the Quark-Gluon Plasma should be higher than that of ordinary matter 
|75|. This effect is called jet quenching |76,77] and has several consequences that could be 



observable in experiments. Of more direct relevance to particle spectra, a comparison of the 
transverse momentum spectrum of hadrons compared to appropriately scaled distributions 
from pp or pp collisions should show a suppression at high-px 4 GeV/c). In addition, 
a quark or gluon jet propagating through a dense medium will not only lose energy, it will 
also be deflected. This will destroy the coplanarity of the two jets with respect to the beam 
axis. The angular deflection in addition leads to an azimuthal asymmetry. One can then 
perform angular correlations among high p± particles to study the energy loss effects of the 
partons in the medium. 

3.1.3 Medium EfFects on Hadron Properties 

The widths and masses of the p, oj and (f) resonances in the dilepton pair invariant mass 
spectrum are sensitive to medium-induced changes, especially to possible drop of vector 
meson masses preceding the chiral symmetry restoration transition. The CERES data from 
S+Au and Pb+Au collisions at SPS showed an excess of dileptons in the low-mass region 0.2 
< M < 1.5 GeV/c^, relative to pp and pA collisions |^, 79|, which has been the subject of 
active discussion. Although the CERES data can be explained by a hydrodynamic approach 
assuming the creation of a QGP [^], alternative scenarios have also provided explanations. 
These have included for example microscopic hadronic transport models incorporating mass 
shifts of vector mesons, and calculations involving in-medium spectral functions (coupling 
the p with nucleon resonances) without requiring a shift in the p mass With the addition 
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of a TPC to the CERES experiment |82|, the resulting increase in resolution (and statistics) 
should help verify or falsify some of the conflicting hypotheses on the origin of the low-mass 
enhancement in the dilepton spectrum. 



3.1.4 Direct Photons and Thermal Dileptons 

The detection of radiation from a high temperature QGP would be an ideal signal to detect, 
as black body radiation is one of the most directly accessible probes of the temperature of 
a given system. In the quark-gluon phase, the gluon-photon Compton process gq jq is 
the most prominent process for the creation of direct (thermal) photons (with additional 
contributions from the qq — > jg annihilation process). Unfortunately, a thermal hadron gas 
with the Compton scattering reaction vrp 7p (and pion annihilation vrvr 7p) has been 
shown to 'shine' as brightly as a QGP (or even brighter still) |^. However, a clear signal 
of photons from a very hot QGP possibly formed at RHIC could be visible at transverse 



momenta in the range 2-5 GeV/c |84-|86||. However, it is also possible that flow effects can 
prevent a direct identification of the temperature and the slope of the p± distribution. WA98 



has observed a direct photon signal Pb + Pb collisions at SPS Comparing the results 
to pA data, they observe an enhancement for central collisions, suggesting a modification of 
the photon production mechanism. 

In addition, dileptons can also carry similar information as photons on the thermody- 
namic state of the medium {thermal dileptons). Since dileptons interact only electromag- 
netically, they can also leave the hot and dense reaction zone basically unperturbed. The 
difficulty of this type of signal is that one does not have a significant feature, such as a 
mass peak. One has to analyze a spectrum which is a convolution of several complicated 
backgrounds on top of the (small cross-section) signal. At CERN-SPS, the expectation is 
that the contribution of hadronic backgrounds to the dilepton spectrum will dominate over 
the QGP radiation. The main backgrounds are, at low masses: pion annihilation, resonance 
decays, ir — p interactions. At high masses, the Drell-Yan process dominates at SPS. At 
RHIC energies there is an additional charm contribution above 2 GeV/cP. There is only 
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a small window, 1 < M < 1.5 where the rates for a plasma (at very high temperatures, 
T 500 MeV) may be dominant. This signature has proved to be a difficult experimental 
observable, but there is a continued effort to improve the sensitivity of the measurements: 
a study of the p± dependence of various mass windows might perhaps help to disentangle 
the different contributions to the spectrum. 

3.1.5 Strangeness Enhancement 

In hadronic reactions, the production of particles containing strange quarks is normally 
suppressed due to the high mass of the s-quark (m^ ~ 60 — 170 MeV/c^) compared to 
u and d masses. In the presence of a QGP, the temperature of the order of the s-quark 
mass and the rapid filling of the phase space available for u and d quarks should favor the 
production of ss pairs in interactions of two gluons ||8^,^. This should be refiected in an 
enhancement of the production of multi-strange baryons and strange antibaryons if a QGP 
is formed compared as compared to a purely hadronic scenario at the same temperature. 
Important observables in this respect are the yields and ratios of strange hadrons (mesons, 
strange and multi-strange baryons and their antiparticles) which allow the determination 
of the relative strangeness equilibrium. To account for incomplete chemical equilibration, 
a strangeness fugacity 7s is introduced in a thermochemical approach. The particle ratios 
can be calculated assuming either a hadron gas scenario or a QGP and a comparison can 
be made of the values thus extracted in conjunction with other model parameters such as 
T and ^b- 

Because strange hadrons interact strongly, their final-state interactions must be modelled 
in detail before predictions and comparisons of strange particle yields can be done. It is also 
important to stress that an understanding of the enhancement mechanism present in pA 
collisions is crucial in order to interpret the signals in AA collisions. 

STAR is currently addressing several of these topics. The large acceptance of the de- 
tector coupled with precise tracking allows for the reconstruction of the decays of strange 



particles [Q. Studies of high-j3^ hadron spectra |91| and angular correlations |£4| are well 
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suited for STAR as the detector has full azimuthal coverage. Additional detector compo- 
nents for future runs, specifically the completion of a barrel electro-magnetic calorimeter, 
will permit studies of dilepton production and J/V' suppression. 

In all cases, the specific observables that are expected to be sensitive to deconfinement 
have to be correlated with the global characteristics of the collision in order to better un- 
derstand their systematics. For example, critical in the debate of the J/V' results has been 
the dependence of the effect on the measured transverse energy {E±) of the collision, an 
observable similar to multiplicity in that it is correlated with the collision geometry, and the 
determination of the number of participants (see Sec. [4.1.1D from the measured E±. It is 
therefore essential to understand the global observables involved in the systematics of any 
QGP signature. In the following chapter, we therefore turn our attention to these global 
hadronic observables, and their relationship to the collision geometry. 
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Chapter 4 



Global Observables and Charged 
Hadron Spectra 



The search for the new state of matter has not been an easy one. There appears to be 
no simple and unambiguous experimental signature of plasma formation, and one of the 
main lessons we have learned in this field is that an understanding of QGP formation and 
an elucidation of its properties will only come about through systematic studies. It is 
necessary to measure nucleus-nucleus (AA) collisions at various centre-of-mass energies, to 
use different beam species, and to make comparisons with reference data. These comparisons 
are especially important since from proton-proton (pp) collisions one can measure basic 
processes in a cleaner environment, and from proton- nucleus (pA) collisions one gains insight 
into the modification of the basic processes by the presence of normal nuclear matter. These 
are required to understand any signal in AA collisions. For a recent review on hadronic 



particle production in nucleus-nucleus collisions from SIS to SPS energies, see |92|. 

Through global distributions one gains insight into the 'kinetic freeze-out' stage of the 
system produced in the collision when hadrons no longer interact and their momenta no 
longer change. These final-state measures supply information that constrain the possible 
evolutionary paths of the system and can help establish conditions in the early, hot and 
dense phase of the collision. 
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Typically, the first information studied in heavy ion collisions comes from the observed 
particle distributions, both the transverse momentum {p±) and rapidity [y) distributions. 
The spectral shapes are intimately related to the underlying collision dynamics. The ex- 
pected behaviour of these distributions from scenarios consistent with a phase transition can 
be tested. While these may not be sufficient to completely answer the question as to whether 
a quark-gluon plasma is found, they are most certainly a necessary first step to provide con- 
sistency with any given scenario, be it QGP or other. The studies of the proposed QGP 
signatures so far, however, have not provided unambiguous evidence for quark-gluon plasma 
formation. It is the general belief that the proposed experimental programme of high energy 
heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Lab will yield key pieces of the puzzle. In the analysis presented here, we focus on charged 
particle distributions from the first collisions measured by the Solenoidal Tracker at RHIC 
(STAR) experiment. These measurements will serve as a baseline for studies of QGP sig- 
natures and as a guidance for theoretical models. We now give some background for the 
distributions that we will present here. 

4.1 Particle Multiplicity 

The negatively charged hadron {h~) particle multiplicity distribution {Nf^~) yields informa- 
tion on both the impact parameter and energy density of the collision. It is not possible to 
directly measure the impact parameter of the collision, so one must use an indirect measure. 
The event multiplicity is one of the observables that is correlated to the impact parameter. 
The idea is simple. Each of the nucleons in the nuclei that participate in the collision pro- 
duces (on average) a certain number of particles. We can calculate in a geometrical model 
the average number of nucleons that participate in the collision (A'^part) at a given impact 
parameter b. We can thus obtain a statistical mapping of {Nf^-) (A^part) {b)- The 
number of participants (A^part) is also called number of wounded nucleons and we will 
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use them interchangeably here. The scahng of the multiphcity with the number of partici- 
pants is typically thought of as a reflection of the particle production due to low momentum 
transfer (soft) processes. There are refinements to this model. At high energy, it is expected 
that there will be an increased particle production from large momentum transfer (hard) 
processes. Hard process cross sections in pA collisions, e.g. the p± distributions at very high 
p±, are found to be proportional to the number of elementary nucleon-nucleon collisions, 
which we call the number of binary collisions (A'^cou)- Some recent models ||9^ include for 
example the assumption that the particle production is derived from a linear combination 
of the soft and the hard processes, i.e. a linear combination of A'part and A'^coU- 

To make the distinction between the two quantities A^part and A'cou clear, since the two 
are related and sometimes lead to confusion, we describe the concepts in more detail. A^part 
refers to the number of nucleons that were hit, or that interacted in some sense, which is why 
they are also sometimes called wounded nucleons. For a head-on (6 = 0) Au -|- Au collision 
assuming the nucleus to be a hard sphere, or rather, to be a bag filled with hard spheres, 
we then find simply Apart = 197 -|- 197 = 394. There will be deviations from this, as one 
introduces a more realistic density profile for the Au nucleus. In addition, the distribution of 
nucleons in the nucleus is not always the same, there are volume fiuctuations and the nuclei 
have Fermi motion. This is typically introduced in the models by adding a parameter to 
represent the size of the fiuctuations. It is also common to find in the literature references 
to the number of participant pairs and comparisons made this way. If one normalizes the 
particle production "per participant pair", it is straightforward to compare to pp or pp data 
which can be thought of as the limit of 1 participant pair. Again, for a head-on Au -|- Au 
collision assuming hard spheres, the number of participant pairs is just A^part/2 = 197. 

AcoU refers to the number of elementary nucleon-nucleon collisions. It includes all par- 
ticipating nucleons (i.e. AcoU > A^part)- The difference can be thought of in the following 
simple picture. Let us follow one particular nucleon through the collision as if it were a 
billiard ball and do its accounting. If it does not interact at all, it does not count for either 
Apart or A'cou purposcs, and we call it a spectator nucleon. If it interacts, then we count 
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it once for iVpart purposes and that is the end of our A^part accounting using this nucleon. 
We of course also count it for the purposes of Nc^n at this point. The difference is that 
there are can still be other nucleons in its path, and if it interacts again, we increase our 
-^coU counter. We do this for every time our original nucleon collides, whereas our iVpart 
counter remains at 1. Since each nucleon in a nucleus can interact many times as it "punches 
through" the other nucleus, it is evident that A'^coU > -^part- The corresponding simple limit 
for A^coU for the hard sphere case is A'coii oc A^part^^^- This can be seen from the following 
argument. It is easy to see that A?part is proportional to A {e.g. Np^it = for central 
collisions). Every nucleon that counts for A?part must also count for ATpon- The number of 
additional collisions that count for A^coii have to do with the number of additional nucleons 
(from the target nucleus) that lie in the path of our original nucleon (from the projectile 
nucleus). As seen from the target nucleus, the path of the projectile nucleon is a straight 
line parallel to the beam axis. A nucleon that punches through the center of the target 
nucleus will cross a length 2R of the target nucleus, where R is the nuclear radius. If the 
nucleus is not in the center, the length will still be proportional to R. For the hard sphere, 
the nucleon density is constant and therefore the additional collisions that sum up to A^coii 
are then proportional to R. Since the nuclear radius is proportional to A^^^, we arrive at 

It is important to stress the limitations of our mapping from A^^- to impact parameter. 
First of all, it relies on the accuracy of the simplified model relationship between the number 
of participants and particle production. In addition, the relationship is statistical only, i.e. 
we cannot experimentally measure the impact parameter of a given event. We can only 
measure the multiplicity of each event, and for a given ensemble of events compute the 
mean multiplicity. Then we can relate this to an average number of participants, since for 
every event, even if we keep the number of participants fixed, there will still be multiplicity 
fluctuations. Furthermore, even keeping the impact parameter fixed does not fix the number 
of participants either, since we expect fluctuations in the initial configuration of the nucleons 
in the nucleus for every event due to Fermi motion as well as small variations in the size of 
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the nucleus. We therefore expect that such a statistical map will probably work best for the 
central collisions {i.e. we know that the highest-multiplicity events must come from central 
collisions) and be less reliable for peripheral collisions where fluctuations will dominate. 

The multiplicity distribution Nf^- belongs to the most global class of observables (or 
rather, of the most integral sort of observables, as we measure the cross section integrated 
over azimuth, over p± and for a wide pseudorapidity slice). The information that can be 
derived from the multiplicity distribution is basically related to whether there are significant 
deviations from the simple geometrical picture of the collision, commonly referred to as the 



Glauber model (see e.g. |38|). We now discuss the relationship between the nuclear geometry 



and the final multiplicity in this model in more detail. 
4.1.1 Glauber Model 

To clarify the meaning of the geometrical model, which can have different a priori assump- 
tions in different implementations in the literature, we discuss here the characteristics we 
consider. The starting point is to assume that the nuclei are composed of discreet and 
point-like nucleons. We distribute the nucleons according to the Woods-Saxon spherically 
symmetric density profile 

p{r) = . (4.1) 

1 + e c 

For the case of Au, the parameters are ro = 6.38 fm, pQ = 0.169 1/fm^ and c = 0.535 fm. 



obtained from eA scattering |g4|. The density is shown in Fig. With these parameters, 
we get a total number of nucleons in the Au nucleus of 

/ p{r) A-Kr"^ dr = 196.6 . (4.2) 
Jo 

It is also common in the literature to take a different value of po such that the integral over 
the density is normalized to unity. We can then interpret p{r) as the probability of finding 
a baryon in the volume element d^x = d'^bdz at the position {b,z). With this convention, 
when calculating quantities such as the number of binary collisions, we see that there will 
be factors of A every time this integral appears. 

32 




2 4 6 m- ''(f"^) 

Figure 4.1: The Woods-Saxon Density profile for the Au nucleus. 

The next assumption is that the interaction probability is just given by the pp cross 
section, neglecting effects like excitations and energy loss. At these energies, the pp and 
the pp cross sections are very similar in value, and can therefore be used interchangeably. 
Unfortunately, there are no measurements of either cross section, cjpp or a^p, at -y/i = 



130 GeV. The experiments UA5 m and UAl ^ at CERN-SPS have reported the pp 
cross section at \/i = 900 GeV. Although UA5 took data at ^ = 200 GeV, in Ref. |||| 
they only measure the ratio of the cross sections at the two different energies and use a 
parameterization given in jo^ (Eq. |4.3|) to get a value for the total cross section. Then, 
based on the ratio (Jei/crtot of elastic to total cross section, they obtain a value for cJinei- 



(jpp {p) = A + Bp-"" + C In^ (p) + D ln{p) . 



(4.3) 
(4.4) 



For Eq. O, the important part at energies above 100 GeV is the In^(s) term. The first 
two terms are used to describe the data at lower energies and the difference between pp and 
pp collisions, E is the beam energy. At ^/s = 130 GeV and above, the pp and pp cross 
sections are very similar. The parameterization found in the Particle Data Book ||3^ is 
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= 200 GeV 


= 130 GeV 


0"tot 
0"el 
"^inel 


52.40 mb 
10.66 mb 
41.74 mb 


49.26 mb 
8.91 mb 
40.35 mb 



Table 4.1: The pp total, elastic and inelastic cross sections obtained from the parameter- 



ization from the Particle Data Book |3G]. The errors on the values are on the order of 
1%. 



given in Eq. 4.4 in terms of the laboratory momentum p. They quote parameters for the 



total and for the inelastic cross section. The relevant numbers are found in Table ^Jj. At 
the energy of the RHIC 2000 run, we obtain a value of (Tpp(130 GeV) = 40.35 mb. For head 
on collisions 6 = 0, as we said before, in the hard sphere limit the number of participants 
will just be 2A where A is the mass number of the nucleus. We need a prescription to 
calculate the overlap at any given impact parameter. This is given by the nuclear overlap 
integral, Taa, which is a calculation of the overlap of the density profiles (in cylindrical 
coordinates, where the z direction is the beam direction) of two specific nuclei at a given 
impact parameter b: 



Taa(&) 



cfs dzi dz2 pi{s, zi) ■ P2{s - b, Z2) 



(4.5) 



We limit ourselves to the case of symmetric AA collisions. The coordinate system is shown 



schematically in Fig. 4^. For a given impact parameter, we calculate the product of the 
densities of each nucleus at a given point s and integrate over all space. We normalize the 
integral such that 

j TAA{b)dH = 1 . (4.6) 

Since d'^b is an element of area, Taa then has units of inverse area. With this definition, we 
can obtain the probability of having n interactions at a given impact parameter 

Pin, b)=(^^^^ (TAA(6)app)" (1 - Taa(&)^pp)^'"" • (4.7) 

The first terms takes care of the combinations of choosing n nucleons out of A^, the second 
term is the probability of having exactly n collisions and the third term is the probability 
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Coordinates for Overlap Integral 




Figure 4.2: Schematic of the colhsion geometry for the calculation of the overlap integral 



Taa, Eq. |4.5| . The area of overlap is the "football" shaped region in the middle of the two 



spheres. The zi and Z2 coordinates are perpendicular to the plane of the paper. 

of having exactly — n misses. The total hadronic cross section for Au + Au collisions is 
then found to be: 



0"AuAu 



j cfb 1 - (1 - TAA(&)fTpp)^' = 7.2 barn (4 



which can be read as 1 minus the probability of not having any collision (n = 0) at a given 
impact parameter (i.e. the probability of having at least one interaction at each b) integrated 
over all impact parameters. 

The mean number of binary collisions (A/'cou) and the mean number of participants 
(-/Vpart) are obtained from Taa at a given impact parameter as 

{Ncon){b) = app-A-A-TAA{b) 

(iVpart)(6) = 2Ay" dhTAAis) |l- (l - Taa(s - 6)app) ^ I . (4.9) 

Since the definition of the overlap integral takes care of counting interactions, it is not 
surprising that the number of binary collisions is simply proportional to Taa- In the above 
definition of A'^coll, the factor comes from our choice of normalization. 
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Mean Wounded Nucleons, 197Au+197Au Vsuu = 13G GeV 



I Mean Binary Collisions, 197Am-197Au vshn =130 GeV 
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Figure 4.3: (A^part ) and (A'^coii ) as a function of impact parameter. 



Fig. shows the resulting statistical relation between the impact parameter and iVcou 
(-^part) in the left (right) panel. For large impact parameter, both A^part and A'coii are close 
to zero. For the most central collisions, A'part is found to be ~ 394 as expected. We see that 
the shapes are similar, although the overall scale is the different. 



Indeed, Fig. shows the statistical relation between A^part and A'^coii as a function 
of impact parameter. The dotted curve corresponds to A'coU = (A'part/2)'^''^. The above 
definitions and relations are the basis of the geometrical model. We obtain a relationship 
between the A^part and A'coU as a function of impact parameter. These are still not directly 
measurable quantities in experiments. In a fixed target environment, one can try to esti- 
mate A^part by placing a hadron calorimeter to measure the forward-going energy, and thus 
approximately determine the spectator nucleons. In the RHIC environment, the closest we 
can get to such a scheme is to place calorimeters in the forward-backward regions, called 
zero degree calorimeters (ZDC's) at RHIC, which detect the spectator neutrons (see Section 



5.3). So there has to be another indirect step to model the particle production based on 



A^part, A'cou, and Taa- There is not a unique way to do this. We will describe the "eikonal" 



approach to particle production given in Ref. |98| which was used already in the analysis of 



PHOBOS data H. 
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Figure 4.4: (A'coU ) correlated to (A'part )• The curve is the Hmit iVcoU = (A^part/2)^/3 . 

We assume that each participant contributes to particle production, i.e. {Nf^~{b)) = 
Q{^pa.rt{b)) where g is a scale factor. This approach has worked well at low energies, but 
at higher energies the hard processes contribute to particle production as well. Therefore, 
in |p3| , the assumption is that particle production scales as a linear combination of A'part and 
A^coU- This model also attempts to incorporate the multiplicity measured in pp collisions. 
The particle production is then given by: 

{Nh^ {b)) = q-n^p {^) ((1 - x) ■ i^^E^ + ^ . (iV^„n(&))) (4.10) 

The idea is that x is a number between — 1 which gives the fraction of particle production 
that scales as ~ A'^coii- We also need the expected particle production per participant from 
pp collisions. It is measured at 200 GeV [|^,^. The number of charged particles in the 
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collision per unit of pseudorapidity at mid-rapidity is reported in Ref. |96| to scale as 



(iNch/(i?/(^/i) = -0.32 + 0.55 ln(^/GeV) , (4.11) 

which at 130 GeV yields Upp = dNh-/d'q\^=Q = 2.36/2 = 1.18 . 

The next assumption is to pick the statistical distribution to model the fluctuations 



about the mean value. One can choose e.g. a Poisson distribution |10C,101| or a Gaussian 



distribution |31,93|. We follow the Gaussian prescription here, we then require an additional 



parameter to get the variance for the Gaussian, which is taken as simply 

al = a-q^-{N^Ab)) . (4.12) 

The parameters of the model are thus a, x and q. 

We are finally in the position to calculate the multiplicity distribution, which is the actual 
experimental observable. This is done by convoluting the various Gaussian distributions 
obtained for each impact parameter, weighted by the appropriate interaction probability 
and integrating over all impact parameters: 

aaAuAu^^^_^ = / (1 - P{0,b)) -^L- e ^9 (4.13) 



Prom such a picture, we expect roughly the following behaviour as a function of impact 
parameter. The cross section is largest for very peripheral collisions (6 ~ 12 fm), dropping 
rapidly at first (6 ~ 12 — 10 fm) and then falling more slowly in the region of mid-central 
collisions (5 ~ 8 — 4 fm), eventually reaching a limit for the most central collisions (^ < 3 fm). 
A schematic curve of the cross section as a function of multiplicity obtained from the HIJING 
model [[102|, |103| is given in Figure [l.5| . The multiplicity is the lowest ordinate axis. A related 



experimental observable is the transverse energy E_\_ in the next axis. The percentage of 
the hadronic cross section is then given. The ordinate axes at the top of the figure are the 
non-measurable quantities which one can relate in the geometrical model to the multiplicity: 
A^part and the impact parameter h. In this model, we obtain for the 5% most central Au -|- Au 



events (A'part)/2 = 172 participant pairs |93] and (Acou) = 1050 using cxinei = 40.35 mb. 
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Figure 4.5: Schematic of a multiplicity distribution from the HiJiNG model, showing the 
relationship between the shape of the distribution and the underlying collision geometry 
(see text). Figure prepared by T. Ullrich. 

4.2 Kinematic Variables: p±, m±, y and 77 

The next observable after a simple measurement of the charged multiplicity is to look at 
the momentum distribution of particles. Particle spectra are often treated separately in the 
longitudinal and transverse directions. 

For the transverse direction one normally employs the transverse momentum {p±) or for 
identified particles the transverse mass 



m± 



(4.14) 



where m is the mass of the particle. 

It is convenient to treat longitudinal momenta using the rapidity 

y = - in 

2 \E-Vz 



(4.15) 



where E and "Pz are the energy and longitudinal momentum of the particle. The use of 
rapidity guarantees that the shape of the corresponding distribution is independent of the 
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Lorentz frame. Under a Lorentz transformation from a reference system i? to a system R' 
moving with velocity f3z with respect to R in the longitudinal direction, the rapidity y' in 
the R' frame is related to y in the R frame only by an additive constant: y' = y — yp, where 
is the rapidity of the moving frame 

For the incident energy of the Au beams at RHIC (7 = 70), the initial rapidity of each Au 
beam is ybcam = ±4.94, since the beams are symmetric. The rapidity of the centre-of-mass 
system, called mid-rapidity, is y = 0; i.e. the centre-of-mass reference frame is the same as 
the laboratory frame for the collider geometry. 

In the limit of a particle whose mass is much smaller than its momentum, the pseudora- 
pidity variable rj is often used. It is defined as 



For the case of negative hadron distributions, since no particle identification is performed, 
this is the variable of choice since one only needs to measure the angle 9 of the detected 
particle relative to the beam axis (polar angle in spherical coordinates, also called dip angle 
in the helix parameterization commonly used in tracking). Sometimes in the literature 
h~ distributions are also presented under the assumption that all particles are pions, e.g. 



1 104, 105 1 . For the present work, the h~ distributions are given in pseudorapidity since this 
is the actual measured observable, the identified vr distributions are given in rapidity. 

The momentum distributions are usually presented in terms of the invariant cross-section 

d^a d^a d'^a d'^a 

E = = = (4.18) 

d^p d(j)dyp±dp± 27rdyp±dp± 2TTdym±dm± 



The double differential in Eq. 4.18| is obtained from integration over (/>, and the last equality 



follows from the definition of transverse mass, Eq. |4.1^ . 

It is important to understand the difference in shape of the y and r/ distributions that 
arises simply from the change of variables. In particular, the Jacobian dy/dr] characterizes 
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the difference between a distribution given as (fiN/dp_\_dri and one given as d?N/dp±dy. 
From the relation Pz = m± sinhy = p± sinhry, the Jacobian at fixed p± is 

dy_ ^ p± cosh 7? ^\p[ 
dr] m_L cosh y E 

From this relation we again see that r/ approaches y for highly relativistic particles. The 
Jacobian in terms of the variables {p±, rj) and {m±, y) is 



dy m? 

■7r{'m±,y) = \ l 5— ,2 

or] y mj_'^ cosh y 



{m^,y) = A/l-— ^-Tir^ (4-20) 



^{P±,V) = (4.21) 
- ' (p_l/?ti)2 cosh^ r/ + 1 



We can then infer from Eq. 4.20| that in the region y ~ 0, there is a small depression in the 



pseudorapidity distribution dN/drj relative to dN/dy. At high energy, where dN/dy has a 
plateau shape, this leads to a small dip at mid-rapidity for dN/drj, with the dN/drj yield 



being smaller by a factor of approximately y^l — m? / {m±)'^ relative to dN/dy. 

4.3 Dynamics from the Kinematics 

To understand the dynamics of relativistic heavy-ion collisions, it is essential to have infor- 
mation on certain basic aspects of the collision dynamics. The basic idea is that particle 
spectra act as kinematic probes. By analyzing the rapidity and transverse momentum distri- 
butions of the particles produced in the collision, one can study the energy density e, pressure 
P, and entropy density s of the hadronic matter formed in the collision, as a function of 



temperature T and baryochemical potential fiB [[106|| . Prom Lattice QCD, the expectation is 
to observe a rapid rise in the effective number of degrees of freedom, expressed by the ratios 
e/T^ or s/T^ over a small range of temperatures |107|. Experimental observables which are 



thought to be related {e.g. in hydrodynamics) to T, s, and e are the average transverse mo- 
mentum {p±), the hadron rapidity distribution dN/dy [ [108 , 10E|, and the transverse energy 



dEj_/dy, respectively [£5|. (The average transverse momentum, or the related slope of the 
transverse momentum distributions, actually reflect not just the temperature in such models 
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but also the transverse expansion of the system.) If there is a rapid change in the effective 
number of degrees of freedom due to deconfinement in the medium, the first hope was that 
one might see vestiges of this effect in a plot of (p^) as a function of dN/dy, which would 
be expected to show a rise, then a characteristic saturation of {p±) while the mixed phase 
persists and then a second rise when the underlying matter undergoes a structural change to 
its coloured constituents This simple picture however has several caveats. We discuss 
here a few of them. We already stated that it is much too naive to identify {p±) or the 
inverse slope parameter of a transverse momentum spectrum directly with the temperature 
of the system. It is also simplistic to assume that the exponential shape arises in an iden- 
tical manner as that of a Boltzmann gas that reaches thermal equilibrium through a series 



of internal collisions of the particles in the system. Statistical thermal models IllO] have 



even been applied to e+e" collisions in fits to hadronic particle spectra [111|. The observed 
spectra for e"'"e~ collisions certainly do not come about because of a hadronic rescattering of 
the final state particles, but are a the result of a sampling of the available phase space, i.e. 



the particles are already born into equilibrium |112|. Furthermore, if the system formed in 
AA reaction does in fact thermalize through collisions, the final state particles observed in 
the experiments, and in particular the momentum spectra, will only reflect the coolest phase 
of the evolution of the system. Stated another way, if thermalization occurs the effective 
temperature extracted from the spectrum (even if we assume we can take care additional 
collective effects such as flow which modify the p± distributions) will probably not carry 
information of the hot initial phase, smearing any structure in a {p±) vs. dN/dy plot. In 
addition, a visible flat structure in the 'T vs. e' diagram necessitates a signiflcant duration 
of a mixed phase, an effect that probably requires the presence of a strong first-order phase 
transition. However, lattice simulations currently favour a more smooth cross over, perhaps 
a second order transition. In this respect, critical phenomena in the form of increased event- 
by-event fluctuations are possibly a more robust observable with respect to the existence 
of a phase transition. Even in this case, fluctuation studies would necessitate probing the 
region near a possible critical point in the T — phase diagram, and current efforts to 
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shed light as to the presence and location of a critical point indicate it is in the large fiB 



region [113|. 



Rapidity and transverse momentum distribution also allow us to address properties of 
the particle emitting source. The degree to which the incoming nuclei are stopped by the 
collision is reflected in the rapidity distributions of produced particles as a shift with respect 
to beam rapidity. At sufficiently high energies for example, a picture due to Bjorken |59|, 



assumes that the mid-rapidity region should undergo an idealized hydrodynamic longitudinal 
expansion. The charged particles found at mid-rapidity would be mainly produced particles, 
the energy of the incoming nuclei would be so great that the collision among the nucleons 
would be insufficient to stop the nuclei, so the incoming nuclei would essentially go right 
through each other. The incoming baryons are then found very close to their initial rapidity. 
More importantly, the rapid longitudinal expansion would have as a result that the rapidity 
distribution of produced particles should be flat around mid-rapidity. 

This contrasts the picture found at AGS energies, where a significant amount of stopping 
is observed and the mid-rapidity region is net-baryon rich. This is reflected in the rapidity 
distribution of charged particles which is peaked at mid-rapidity, as are the net-proton 
distributions [p — p). Full stopping as in the model of Landau Hl08| is expected then to 
work at low energies, and the rapidity distribution of produced particles at mid-rapidity is 



found to more closely follow a Gaussian shape ||llj-[ll6|. Since the observed width of the 
distribution is significantly narrower than that observed for lighter systems, this has been 
understood as evidence for strong baryon stopping. 

Prom the rapidity distributions, we can therefore obtain information for example, as to 
whether the source is spherically symmetric or elongated, whether it is static or expanding, 
and the degree of longitudinal and transverse expansion. The momentum distributions, 
when paired with elliptic and transverse flow studies, can also test if there is significant 
collective behaviour in the system and approach the question of whether the system reaches 
thermal equilibrium. Experimentally, one can test the hypothesis of boost invariance by 
measuring rj and y distributions. If the system is boost invariant to a certain extent and 
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dN/dy is flat as a function of rapidity in some phase space range, the h pseudorapidity 



distribution should be well approximated by Eq. 4.21 . The rapidity distribution of identified 



particles should simply be constant as a function of y. The slope parameter of vr distributions 



have also been found to show a significant rapidity dependence at lower energies ||114| - |117| . 
We can measure the rapidity dependence of the tt spectra at RHIC and further test the 
hypothesis of boost invariance. These are of course not sufficient conditions to establish 
boost invariance, but they are necessary in case it does hold. 

4.4 Overview of Transverse Momentum Spectra 

The transverse momentum spectra possess many interesting features in AA collisions. The 
use of the transverse mass is sometimes preferred because experimentally the cross- 
section {X/pi_)[da / dpx_) of a given particle species is better described by an exponential 



in m_L rather than in p_\_ [118|. There are contributions to the spectrum that come from 
the various physics processes of interest. In order to extract correct information about the 
collision, it is necessary to take these into account. Even now, our understanding of all 
the features of the spectra is incomplete. A brief overview of the main features of the p± 
distributions will be discussed. 

A general feature emerging from measurements of transverse momentum distributions 
in central collisions of heavy nuclei is that the invariant distributions d'^N/mj_dm±dy are 
approximately exponential, i.e. d'^N/m±dm±dy oc exp(— m_|_/T). One approach to analyze 
the p± distributions has been to treating the system as a hadronic gas, and in particular 
to use the measured {p_i_) to estimate the freeze out temperature T of the gas and its 
transverse flow velocity These also have inverse slope constants T increasing linearly 

with the mass of the particle under study. This has been interpreted as evidence for collective 
transverse flow: since p± = mj^^, if there is a common flow velocity superimposed on 
the random thermal motion of particles ^t^ermal^ ^.j^g slope constant T, which is proportional 
to {p±) increases linearly with m. It is important to note, however, that in these models. 
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T and are coupled: a higher T can compensate for a lower and vice versa. This 
approach relies only on the analysis of p± spectra. It must also be noted that in a plot of 
T vs. mass, it is difficult to make meaningful comparisons of slope parameters when they 
are obtained from different phase space regions, and there is the danger of obtaining a very 
different slope simply by fitting in a different part of the p± spectrum. Therefore, rather 
than concentrating on the slope parameters, more realistic models incorporate additional 
observables such as two-particle correlation data |119|| and attempt to describe both the 
observed p± distributions and the observed correlations to obtain T and 

The slope of the p± distribution is also seen to increase in going from pp to AA. The 
spectrum of the pion, due to its low mass, is not as affected as the spectra of heavier particles 
by a given collective flow velocity. Thus, the pions are good probes for studying thermal 
properties at freeze-out. These features of spectra can be measured at RHIC energies, and 
the linear scaling of spectra with the mass of the particle checked. The scaling could very 
well be reduced given that RHIC energies are a significant leap with respect to SPS energies, 
and the collision process can occur in a short enough time scale such that particles involved 
might not have time to develop a significant fiow velocity. On the other hand, an increase 
in the radial fiow (and also of the elliptic fiow ^2) would signal an significant amount of 
collective behaviour. Uncertainties remain, and they will only be clarified once we measure 
spectra and compare them with the various scenarios: from superposition of pp collisions 
governed by measured hadronic cross sections; to relativistic microscopic models based on 
string formation and string fragmentation (where a string is represented by the color fiux 
tube from a quark and a diquark. In longitudinal exchange string formation, the excitation 
of the string originates from a stretching of the original partons in the hadron caused by a 
large longitudinal momentum exchange from the hadronic collision. HIJING, FRITIOF and 
RQMD use this scheme.); to hydrodynamics, which can be thought of the limit where the 
mean free path of the constituents is zero and the system is a fiuid. 

It is also known that there are significant deviations from a purely exponential spectrum. 



both at low and high-p^ | 120 |. The high-p^ region (at SPS, p± > 1.5 GeV/c) shows 
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an enhancement present already in pA collisions known as the Cronin effect |121]. This 
enhancement shows a strong dependence on both target and projectile mass, and a weak 
rapidity and y/s dependence. It is thought to be partly due to multiple low-momentum 
scattering at the partonic level inside nuclear matter. Perturbative QCD calculations for 



pA have been done [122| and agree reasonably with data. However, a complete quantitative 
understanding of this effect has not yet been achieved. The CERES collaboration measured 
spectra in pA collisions, found a similar effect when comparing to pp, and reached similar 
conclusions. Their S+Au data indicated that the onset of this enhancement occurred at 
lower p± and that the overall enhancement was larger. This was taken as evidence that in 
AA collisions, there is an increased number of scattering processes per parton compared to 
pA collisions, leading to increased thermalization of the collision system. In addition, the 
shape of the high p± spectrum was very close to being exponential. These two observations 
suggested that a large degree of thermalization is reached very early for central S+Au 
collisions at ^/s^ ~ 20 GeV. 

The low-p_L region shows an enhancement as well. This enhancement, however, shows 
different systematic behaviour than the high-p^ excess. The underlying physics that give rise 
to this enhancement are different than in the high-p^ region. In fixed target experiments, the 
low-p± enhancement was observed to have a strong target dependence but weak projectile 
dependence, as well as a strong rapidity dependence. These effects are not yet completely 
understood. There have been several attempts at explaining the enhancement (at least 
qualitatively). Resonance decays are thought to contribute significantly in this region [|123| , 
|124|| , and their contribution to the low-p± enhancement would also be consitent with the 
observed strong dependence on target and rapidity. Produced particles can also rescatter in 
the target nucleus, and the ensuing cascade would yield additional soft pions and resonances. 
In addition, at the lowest p± one should expect that the treatment of a radiating gas based 
on Boltzmann statistics must break down, and one should really use either Bose-Einstein or 
Fermi-Dirac distributions to fit the data. For pions, a Bose-Einstein fit will give a natural 
low-p± enhancement compared to a Boltzmann fit to the spectrum. Other processes can 
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also contribute to the low-p_L region, but our understanding of the enhancement is by no 
means complete at this point. Studying this region at RHIC energies will certainly help to 
shed light on the subject. 

To fully understand the non-hadronic effects that may arise in an AA collision, studying 
the particle distributions, in order to obtain basic information about the collision dynam- 
ics, is a necessary step. Several approaches have been proposed to describe the observed 
transverse momentum spectra, and to obtain different types of information from them. In 
the absence of any collective behaviour, the spectra should be a simple geometrical super- 
position of nucleon-nucleon scatterings. A deviation of the particle distributions observed 
in AA collisions from those observed in pp or pp collisions is already strong evidence for 
some form of collective behaviour. We discuss two models that have almost diametrically 
opposite assumptions. 

In hydrodynamic models, the colliding system is treated as a fluid, and can be thought 
of as the limit of zero mean free path among the constituent partons. The transverse mo- 
mentum distributions are obtained from a transverse flow velocity profile of the component 
particles of the fluid. The profile will depend in general on the initial conditions (central- 
ity, initial energy density) as well as the underlying equation of state (including or not a 
phase transition) as a function of time. This models typically yield the strongest signals for 
collective effects such as elliptic flow an longitudinal expansion. They should work best at 
low transverse momentum, and their signals can be tested for example by measuring the 
centrality dependence of the elliptic flow signal. 

In partonic cascade models, the energy deposition in the collision is broken up into scat- 
tering of the constituent partons of the colUding nuclei. They rely on the framework of 
perturbative QCD, following the perturbative interactions among partons until they reach 
thermal equilibrium. These models predict a very rapid thermalization, mainly due to radia- 
tive energy degradation, which yields a transverse momentum distribution that is exponen- 
tial to a high degree. In all these theoretical models, the shape of the particle distribution 
yields information to the underlying collision dynamics. Since they rely on perturbative 
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calculations, they should work best for hard processes such as jets and high-p^ momentum 
spectra. 

There are, of course, other observables that need to be studied to provide further insight 
into the full picture of the reaction process. These include studies of strangeness produc- 
tion, HBT correlation analyses, and fluctuations of observables on an event-by-event basis. 
These and other questions are also being addressed in STAR. The particle distributions are, 
however, a very basic building block of the collision picture, and one of the first that can 
and has to be addressed in the initial stage of the experiment. 

In previous heavy-ion experiments at lower energies, the study of the particle spectra 
has helped to provide guidance into interesting aspects of the collision dynamics, and has 



also provided a means to compare results with model calculations [125]. The usefulness of 
studies with particle spectra is best illustrated with a few examples. 

4.5 Previous Studies Using m± Spectra 



At the BNL-AGS, E877 pT4|| measured hadron spectra from Au+Au colhsions at 10.8 GeV/c 
per nucleon. They measured proton and pion m± distributions as a function of y and 
performed exponential fits to extract inverse slope constants from the spectra. For both pions 
and protons, the spectra showed deviations from the simple exponential, e.g. an excess for 
pions was observed at m_L — m,r < 0.2 GeV/c^, with the excess increasing systematically as 
one approached mid-rapidity. They compared the Au -|- Au inverse slope constants extracted 
from the measured m± spectra with a lighter system, Si-|-Al. The slope parameter in the 
Au -|- Au case systematically increased when going from the fragmentation region to central 
rapidities. The data indicated that the maximum the observed slope parameter would be 
reached around mid-rapidity (e.g. the measured inverse slopes for protons at y/uheam = 0.7 
were T^s ~ 150 MeV for Si+Al and T^s ~ 250 MeV for Au + Au). This could be interpreted 
as being due to a larger collective transverse flow component in Au-hAu compared to that 
of Si-hAl. The proton rapidity distribution was also measured for both systems, being much 
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narrower in Au+Au. This provided information as to the degree of stopping in the colhsion. 
The greater stopping observed in the Au+Au system was interpreted as being a consequence 
of the smaller surface to volume ratio as well as the increased average number of rescatterings. 
E877 also measured pion spectra in Au+Au collisions. The rapidity distributions were well 
matched by calculations from RQMD, an event generator based on hadronic rescattering 
using a mean field approximation ||126|] . The m_L spectra, and specifically the low-m^ region 
{m± —uItt < 0.2 GeV/c^) showed an enhancement above the pure exponential. In addition, 
the deviations were seen to be systematically larger for vr" than for vr"'". This asymmetry 
was also observed to systematically decrease as a function of rapidity, an effect seen also 



by the E866 collaboration [127|. This was interpreted as due to the different Coulomb 
potentials seen by the different charge types at freeze-out. (Coulomb interactions are also 
relevant to the interpretation of particle correlation studies.) The study of the influence 
of the Coulomb effect on the shape of the particle spectra provided E877 with a different 
approach to determine the spacetime particle distribution at freeze-out. 

At the CERN-SPS, several experiments have measured charged hadron spectra. NA49 
|119|| in central Pb+Pb collisions measured mj_ spectra of charged pions. They also measured 
two-particle correlation functions to deduce parameters from the collision dynamics. As 
mentioned earlier, in a typical fit to an m± distribution, one can trade lower values of 
^flow higher values of T. NA49 applied a model including longitudinal and transverse 
expansion [128,126] using both the h~ and deuteron transverse momentum spectra as well 



as results from two-particle correlations. They were thus able to constrain the region for the 
freeze-out temperature T and the transverse fiow velocity to T = 120 ± 12 MeV and 
^flow _ 0.55+0.12 [|119|| . Comparison of the charged hadron spectra to predictions of a simple 
hydrodynamical model [125|, yields similar values for T and assuming a quark- gluon 
plasma first order phase transition at a temperature ~ 160 MeV (where Tc was fixed 
by the lattice results then available). Under these assumptions, the data favor a relatively 
low kinetic freeze-out temperature compared to the critical temperature, consistent with the 



picture from Fig. 3.1 
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CERES 1 104 1 in Pb+Au collisions measured m± spectra of identified pions and negative 
hadrons as well as 'net proton' — h~) spectra. Fourier analysis studies of the azimuthal 
charged particle distribution with respect to the reaction plane were performed. The am- 
plitude of the first harmonic, vi, is associated with directed flow. The second harmonic 
coefficient, ^2, is called elliptic flow. The pseudorapidity dependence of these two coeffi- 
cients was studied as a function of centrality, giving strong evidence for collective behaviour 
in these collisions. In addition, a comparison of the 'net proton' spectrum from CERES 
with random-walk models allowed one to conclude that these models were not suitable to 
describe the flow-like features observed in the transverse momentum spectra. 

For many of these measurements, a large acceptance detector such as STAR is ideal. We 
now proceed to discuss the STAR experimental setup for the RHIC summer 2000 run. 
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Chapter 5 

The STAR Experiment 



The main thrust of the STAR detector is dedicated to the heavy ion colhsion programme. 
It must cope with the large event multiphcities associated with central heavy ion collisions, 
and with the large interaction rates expected at RHIC. The collider was designed to deliver 
a luminosity of £ = 2 x 10^^ cm'^ s'^ for Au + Au colhsions at ^/s^ = 200 GeV. STAR 
is also designed with several other physics programmes in mind. It has the capability to 
measure the products from ultra-peripheral heavy ion collisions in order to study coherent 
photon and pomeron interactions from the intense electric fields originated by the Au nuclei 



|130| , 131 ]. There will also be pp and pA collisions at RHIC to obtain information on nuclear 
parton distribution functions and for heavy ion reference data. In addition, there is a 
strong programme of polarized pp collisions starting in 2001 which aims to measure the 



contributions from the gluon (and from sea-quarks) to the total spin of the proton ||132|| . 

This presents a challenge, since STAR has to deal with a wide dynamic range in several 
areas. For heavy ion collisions, the number of particles to be reconstructed in the detector 
vary from less than 10 for peripheral Au -|- Au events to as many as ~ 1000 for central 
collisions. The interaction rates are large: RHIC reached its goal for the first year of 
delivering 10% of the design luminosity. The collider is designed to operate with 60 bunches 
in each ring. The time between beam crossings assuming all 60 bunches are filled and taking 
into account the perimeter of the RHIC rings (3.83 Km) is — 0.21 ^s. 
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The rates will be even larger for polarized pp collisions (2 x 10'^^ cm^^ Coupled 
to the large amount of information recorded in the TPC, STAR must handle a large rate to 
tape (20 Mb/s for the 2000 run, double the number for 2001). As this would be the first year 
of RHIC operations, this also presented an opportunity: a wide range of available energies 
would now fall under experimental scrutiny. The excitement of getting the first glimpse 
at the data taken by the experiments was felt throughout the community, and was well 
matched with the struggle to understand the systematics of the different detector systems. 

The STAR experimental setup for the year 2000 run centered around charged particle 
tracking using a large acceptance Time Projection Chamber (TPC). The trigger detectors 
were a pair of hadron calorimeters placed in the very forward/backward direction and a 
barrel of scintillator slats surrounding the TPC. In addition, a Ring Imaging Cerenkov 
detector (RICH) placed at mid-rapidity provides particle identification for high momentum 
particles. The STAR experiment with the main detector subsystems which are in place for 



year 2001 data is illustrated in Fig. p.l| . In place for the year 2000 were the Central Trigger 
Barrel, the Zero Degree Calorimeters, the magnet, the TPC and the RICH (the ZDC and 
the RICH are not shown, as the ZDC is far from the detector center and the RICH is located 
under the TPC and obstructed from view in the figure). 

5.1 Magnet 

Of critical importance to the experiment, the STAR magnet is designed to provide a very 
uniform field parallel to the beam direction. It is constructed as a large solenoid which 
must house the main subsystems of STAR: the TPC, central trigger barrel, and RICH. The 
strength of the field can be tuned from up to 0.5 T, the strength used for the 2000 run 
was set to = 0.25 T. Having a uniform field is desirable for it greatly simplifies the track 
model used in the pattern recognition in the offline analysis: all charged tracks follow a 
helical trajectory to first order. The field uniformity at full strength was better than ±40 
Gauss in the radial direction and ±1 Gauss in the azimuthal direction. The magnetic field 
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Figure 5.1: The STAR Detector 



was mapped before TPC installation to a precision of 1-2 Gauss for all components of the 
field. This allows a calculation of the distortion effects on tracks due to field inhomogeneities 
to a precision of ~ 200 — 300 fim. 

5.2 STAR Time Projection Chamber 

The main tracking detector for STAR is a large Time Projection Chamber (TPC) with 
complete azimuthal acceptance. With the magnetic field of 0.25 T, the p± acceptance for 
charged particles starts is p± > 50 MeV/c. Particle identification is achieved by measuring 
the ionization energy loss applying topological cuts, or reconstructing invariant masses 
of pairs of particles. 
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The TPC tracking volume is 4 meters in diameter and 4.2 meters long. It is shown 
schematically in Fig. 5^. The inner radius of the tracking volume starts at 50 cm. The 




outer radius subtends a pseudorapidity interval of ±1 unit and the corresponding range for 
the inner radius is ±2 units. Thus, the useful tracking volume extends out to ~ ±1.5 units 
due to a requirement of having at least 10 padrow hits in order to reasonably reconstruct a 
track. For the analysis discussed in this work, we focused on the tracks crossing the entire 
tracking volume to reduce systematic effects from variations in acceptance. The tracking 



volume is split in two along the beam direction as shown in Fig. |5^. The two halves 
are divided by a high voltage cathode located at the center of the TPC (z=0), labelled 
"high voltage membrane" in the figure, but called simply the central membrane among the 
collaboration. The TPC gas is a mixture of 90% Ar and 10% methane (PIO). Part of the 
considerations to choose PIO were that the gas should not attenuate the drifting electrons 
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and it must be pure enough to prevent other modes of electron loss due to attachment on 
oxygen and water molecules. This means that the oxygen concentrations should be kept 
below a hundred parts per million. These stringent standards have already been met by PIO 



in the past ||133|| 



The path of a track crossing 150 cm of this mixture is equivalent to 1.17% of a radiation 
length. The tracking volume is surrounded by an electrostatic field cage which is built with 
11.5 mm wide rings to divide the voltage evenly, achieving a uniform decrease from -31 
kV at the central membrane to V at the ground wires. Secondary electrons produced by 
charged tracks ionizing the gas will drift in the E field away from the central membrane 
to the nearest endcap. The field cages are very thin: 0.62% radiation lengths for the inner 
field cage and 1.26% radiation lengths for the outer field cage. The drift velocity of the 
electrons was measured to be 5.44 it 0.01 cm//is, where the 0.01 here denotes the size of 
variations over a period of several days. For a given run, the drift velocity was monitored to 
a precision of ~ 0.001 cm//is. For such a large TPC, the gas and magnetic field strength are 
also important to maintain the diffusion of the drifting electrons to acceptable levels. The 
transverse diffusion in PIO at full field (0.5 T) is ax ~230 ^um/ y^cm or equal to a width of 
3.4 mm after the drift; and this sets the scale for the dimensions of the pads under the anode 
wires, described below. Similarly, the longitudinal diffusion is ctl ^ 360 /im/^cm. For a 
cluster that drifts the full length of the TPC, this translates into a longitudinal width a = 
0.52 cm , corresponding to a drift time of ~95 ns. This determines the electronic sampling 
rate of 10 MHz, or every 100 ns, in order for a simple three point tracking and clustering 
algorithms to work well. 

The drifting electrons are amplified by a grid of wires on each end of the TPC and the 
pulses are read out on small pads placed behind the anode wires. Each endcap of the TPC 
is segmented in the azimuthal direction like the face of a clock into 12 pad arrays, called 



sectors. This is illustrated Fig. |5.3| which shows the end view of the TPC. The beam {z axis) 



goes through the center of the TPC, perpendicular to the page in this view. 



The geometry of one of the sectors is shown in Fig. 5.4. Each sector is divided into two 
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Figure 5.3: End view of the STAR TPC, where we see the 12 sectors covering the full 
azimuthal range. 

parts, an inner and an outer sector. The inner sector has 1,750 pads, each pad measures 
11.5 mm in the radial direction and 2.85 mm in the tangential direction. The inner sector 
pads were made as small as practical in order to provide closer space points for an improved 
two-track resolution. The size is limited by the diffusion limit of the TPC. They are grouped 
into 13 pad rows located between 60 cm and 116 cm from the beam line, measured from the 
center of the padrow. The outer sector has 3,940 pads, and the corresponding dimensions 
are 19.5 mm (radial) and 6.2 mm (tangential). The outer sectors spans 32 pad rows between 
127.2 cm and 189.2 cm radius. The large outer sector pads completely cover the area under 
the anode wires. The gap between adjacent pads in the same row is 0.5 mm for both the 
inner and the outer sector. 

The wire geometry is depicted in Fig. ^ where we see a cross-sectional view of the wire 
layers and of the pad plane. There are 3 layers of wires terminating each side of the chamber. 
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Figure 5.4: Schematic of a full TPC sector with pads. The first 13 rows comprise the inner 
sector, and rows 14-45 are the outer sector pads. 



The first layer is the gating grid. The wires of the gating grid are arranged to terminate 
the field cage voltage and it is therefore used to prevent ionization from reaching the anode 
wires of the chamber when no trigger is detected. The gating grid also helps minimize the 
space charge effects due to the presence of positive ions with very low mobility. The second 
layer is the ground plane, 6 mm below the gating grid, which captures the ions from the 
amplification region. The third layer is the anode plane, 2 mm (4 mm) below the ground 
plane in the inner (outer) sector, where gas amplification is achieved. The induced signal is 
picked up on the pad plane located 2 mm (4 mm) under the anode wires in the inner (outer) 
sector. In order to keep the signal to noise ratio constant at 20:1 for both pad sizes, the 
anode voltages were set to achieve a gain of 1100 on the outer sector and 3000 on the inner 
sector. 

Combining the inner and outer sector, there are a total of 5690 pads per sector which 
corresponds to a total of 136,560 channels for all 24 TPC sectors. The signal measured 



on the pads is amplified using custom CMOS integrated circuits |134]. These consist of a 
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Figure 5.5: Side view of the inner and outer sector of the STAR TPC showing the wire 
geometry. All distances are shown in mm. 



pre-amplifier and shaper circuit followed by a Switched Capacitor Array (SCA) and a 10 
bit Analog-to-Digital Converter (ADC). These are responsible for sampling the arrival of 
the electrons into at most 512 time bins and to digitize the signal on the pads. Including 
the segmentation in the time (drift) direction, the TPC can be thought of as a large 70 
megapixel 3D digital camera. The position of the particle along the drift direction is then 
reconstructed by converting from time bin to position by knowing the drift velocity. The 
position resolution along the drift direction for the summer 2000 run was on the order of 500 
/xm. The front-end electronics are bundled into groups of 16, with the pre-amp, shaper and 
buffer in one chip, and the SCA and ADC in another chip. Two groups of 16 channels make 
up each front-end board, making a total of 32 channels and 4 chips per board. The data 
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are sent to the Data Acquisition system (DAQ) over optical fibers. The full TPC DAQ rate 
is designed to be read out once per second for central Au + Au collisions, which produces a 
rate to tape of ~20 Mbytes/s. 

For a more complete description of the STAR TPC see | 135| , 136|. More detail on the 



TPC operation will be presented in Chapter ^ where we describe the physical processes 
of the TPC that were simulated for Monte Carlo studies. For more information on drift 
chambers, the standard reference is |133|. 



5.3 Trigger Detectors 

The trigger detectors are an array of scintillator slats (CTB) arranged in a barrel at the outer 
diameter of the TPC; and two hadronic calorimeters (ZDCs) at ±18 m from the detector 
center and at close to zero degrees relative to the beam axis. The ZDC units subtend an 
angle of ~ 2.5 mrad from the interaction point. The calorimeters are designed to measure 
neutrons emitted from nuclear fragments of the spectator matter, i.e. matter that did not 
interact in the collision. In contrast to fixed target experiments, where one can measure all 
spectator nucleons in the forward direction, at RHIC the proton trajectories will be deviated 
by the dipole magnets in the beam line. The same dipole magnets that are used to steer 
the beam head on at each interaction region will bend the charged fragments away from the 
ZDCs. Protons of the same momentum as the Au ions are actually deviated by a larger 
angle, since they have a larger charge-to-mass ratio. Ergo the ZDC signal is produced by 
neutrons. Nevertheless, they still provide possibly the best determination of the collision 
centrality at RHIC, and are used for triggering (along with some other observable based on 
multiplicity) in all four experiments. This can serve as a standard for comparisons among 
the different results. 

During the Summer 2000 run, RHIC delivered collisions between Au nuclei at ^s^^ = 
130 GeV. The data presented here are from a minimum-bias sample, triggered by a coin- 
cidence of signals above threshold in both ZDCs with the RHIC beam crossing. The ZDC 
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threshold was set to ensure efficient detection of single spectator neutrons. The efficiency 
of the ZDC coincidence trigger for central events was measured using a high-threshold CTB 
trigger. The trigger efficiency was found to be above 99% for the entire range of multiplic- 
ities reported in this study. The CTB was used to trigger on central events, as the signal 
is correlated to the multiplicity at mid-rapidity. The threshold for the central trigger using 
the CTB was set to obtain the events with the 15% highest CTB signals. The maximum 
luminosity achieved in 2000 was ~ 10% of the design luminosity. This translates into a 
hadronic interaction rate oi R = ctauAu x £ ~ (7.2 x 10^ mb) x (0.02 mb~^s~^) = 144 Hz. 
For higher luminosities, the trigger will certainly play an important role. 

5.4 RICH 

A Ring Imaging Cerenkov detector (RICH) was placed outside the CTB. The RICH covers a 
smaller area, 1 m x 1 m and is designed to provide high precision velocity measurements for 
enhanced particle identification at high momentum. The momentum information is provided 
by the TPC. One then uses this information to search for Cerenkov photons in 3 different 
annulus regions. For tracks with a perpendicular trajectory to the RICH, the regions will be 
circles of different radii depending on the particle's velocity, which for a known momentum 
uniquely determines the mass of the particle. In the presence of the magnetic field, the 
trajectory of the tracks will not be normal to the RICH detector, so the Cerenkov photons 
will not be found in a circular region, but one can still determine the bounds of the search 
region analytically. The RICH detector is not used in this analysis, but it is the subject of 



another STAR thesis |137|. For more information on the RICH detector see |138|. 

The detector upgrades planned for the upcoming STAR runs include a Silicon Vertex 
Tracker (SVT), an Electromagnetic Calorimeter to cover the barrel at mid-rapidity around 
the TPC and one of the endcaps, a Time-of-Flight patch replacing one of the CTB slats, two 
Forward TPCs (2.5 < \r]\ < 4) and additional Beam-Beam Counters for the trigger system. 
The hermeticity of the STAR detector when all the components are installed is shown in 
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Figure 



5.6. 
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Figure 5.6: Hermeticity of STAR, including all detectors. Note that the RICH occupies real 
estate that will used for the full barrel EMC when completed, so this plot only illustrates 
their coverage. 
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Chapter 6 



Reconstruction and Calibration 



We discuss here the STAR offline reconstruction and analysis software relevant for the re- 
sults presented in this work. The offline simulation and reconstruction software is a major 
component of the experiments, as we rely on it to reconstruct the collision event to a suffi- 
cient extent that the physics goals and the physics capabilities of the sub-detectors can be 
realized, to sufficiently evaluate and visualize the results to determine their adequacy and 
correctness, and to generate acceptance and reconstruction efficiency tables, among other 
things. The offiine software includes all major sub-detectors. We will focus on the recon- 
struction done in the TPC, as this was the main detector used in STAR during the first 
run and is the relevant component of the results discussed here. Then we discuss the STAR 
global chain, where vertex finding and propagation of tracks to find the primary track can- 
didates is done. The other major detector sub-system in the year 2000 run was the RICH, 



and it is discussed in another dissertation |137| . The original design for the STAR offiine 



simulation and analysis software is given in Ref. |139| 



6.1 Event Reconstruction in the TPC 

The task of the TPC reconstruction software is to reduce the raw data taken in an event 
(nearly 7 million ADC values) to lists of meaningful quantities such as space points, particle 
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tracks, vertices, etc. It must employ the best knowledge of the relevant calibration param- 
eters for each run. These include for example the drift velocity of the electrons in the gas, 
trigger time offsets, temperature, pressure, magnetic and electric fields. 

TPC Simulations and Analysis Chain 









Event Gen. 
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Figure 6.1: Schematic of the STAR TPC Software chain. 
The procedure employed here is typical of that employed in many experiments. A flow 



diagram of the different software modules is depicted in Fig. |6.1| . The chain is designed 
to cope with simulated and real data using the same code. We focus on the real data 



reconstruction. The details of the detector simulation (labelled "slow simulator" in Fig. |6jJ) 
are discussed in Chapter |7[ 

One starts by examining the pixel information for every padrow independently in order 
to find the location where a charged particle has crossed the row, which we call a hit. A 
tracking algorithm then links together series of hits to find particle trajectories through the 
TPC volume. Information extracted from the ionization or from the topology of the tracks 
is then used for identification. We now discuss the reconstruction, emphasizing the two 
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parts that are important for the charged hadron distributions discussed in this work: the 
determination of the momenta and the identification via 

6.1.1 Cluster Finding 

The hit reconstruction software converts raw pixel data into reconstructed space points. 
The raw pixel data is arranged into time ordered sequences for each pad in a row, where 
only pixels above a threshold are taken into account. At this stage, one must be careful to 
apply gain corrections and remove channels tagged by DAQ or previous offline runs as bad, 
noisy or dead. Any relative timing corrections between pads are also applied. The cluster 
finding algorithm then looks for contiguous pixels in a 2D space given by the time direction 
and the pad-row direction. After the 2D cluster finder stage, a pass is made to find single 
hits and multiple hits in a cluster. The single-hit finder assumes that the cluster is made 
by the crossing of a single track and estimates the centroid of the cluster. The multiple- 
hit finder tries to find the local maxima in a cluster and may deconvolute the cluster into 
individual hits. Deconvolution of close hits is of critical importance for two-track resolution 
and ^ measurement. The fitted hits are then transformed into TPC space points, taking 
into account the drift velocity, trigger time offsets, sector geometry, and electronics shaper 
response. The space points then contain information on the position of the hit in the STAR 
global coordinate system and on the energy deposited by the track. 

6.1.2 Track Finding 

The detection of charged particle tracks is, of course, the raison d'etre of a Time Projection 
Chamber. The track finding algorithm performs the critical task of converting the recon- 
structed space points into particle tracks, and of determining their 3-momentum. Efiicient 
tracking and good momentum resolution are an essential part of the analysis presented here, 
and especially for analyses such as high-p^ spectra and resonance decay reconstruction. 

The tracking algorithm starts with overlapping hits in the outer pad-rows as the track 
seeds. The first stage of the algorithm consists of a follow- your-nose approach which produces 
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a collection of track segments, i.e. collections of space points associated with a track, either 
complete or partial. Once we have the segments, we perform an initial segment fitting, 
which incorporates multiple Coulomb scattering and energy loss in the gas. Outlier space 
points, i.e. hits that were assigned to the track in the segment finding stage but whose 
position lies far from the allowed track extrapolation error after the fitting stage, may be 
removed at this point. A further extrapolation of the tracks for purposes of segment joining 
is then executed, including also multiple Coulomb scattering and the effects of energy loss. 
The tracking parameters are updated for the joined segments and a clean-up of fragments 
is done in the final filtering step. 

6.1.3 Particle Identification: ^ 

ax 

The TPC allows us to separate pions using their ionization energy loss. The charge collected 
for each hit on a track is proportional to the energy loss of the particle. We parameterize the 
particle's trajectory in order to obtain a path length in the gas which corresponds to each 
segment, i.e. the deposition that corresponds to a specific TPC hit. For a track crossing the 
entire TPC we thus obtain 45 ^ samples, which are distributed according to the Landau 
probability distribution. One of the properties of this distribution is that its tail dies off very 
slowly, and the dispersion of values around the mean is very large (infinite in theory). A 
typical procedure in order to reduce fluctuations from the long Landau tails is to truncate the 
distribution. In the case of STAR, we used 70% truncation, i.e. the highest 30% ionization 
values were discarded. Using the remaining values, a truncated mean is computed, and 
this becomes the basis for any analysis using identified particles in the TPC. The measured 
truncated mean for negatively charged primary particles is illustrated in Fig. |6.2| , shown as 
a function of the momentum. 

The curves are the Bethe-Bloch parameterization used in the analysis for the different 
particle hypotheses. We see that at the lowest momentum, the pions have a greater ionization 
energy loss than the electrons which are already in the saturation region of the curve. 
The pions cross over the electron band at ~ 150 MeV/c reaching a minimum at about 
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Figure 6.2: Particle identification via in the TPC. The plot shows the truncated mean 
value for negative tracks as a function of the reconstructed momentum. The expected 
ionization energy loss curves for e~ , vr", K~ , anti-protons and anti-deuterons are also shown. 

300 — 400 MeV/c. The pions in the relativistic rise merge with the Kaons, which are still 
in the region, at about 1 GeV/c. However, due to the width of the bands, a single vr^ 
peak is only discernible up to ~ 750 MeV. 



^ Resolution 



The ^ resolution for this analysis was found to be 8.8% for tracks with 45 fit points 



(31 ^ samples after truncation). This is illustrated in Fig. 6^. A widely used empirical 



expression ||140|| describes the resolution at minimum ionization, for ^ samples of length 
h cm in argon gas with up to 20% CH4 at a pressure of P atmospheres: 



'^dE/dx 



0.47 



{Ph) 



-0.32 



0.38 



(6.1) 



For the last equality, typical operation of the STAR TPC is at 1 atm and we assume the same 
radial pad length of /i = 1.95 cm (outer sector pads) for all ionization samples. Equation 6J_ 
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Figure 6.3: The ^ resolution as a function of the number of ^ samples. 

then yields an estimate of the best achievable resolution of ~ 7.8% for tracks with 31 ^ 
points. Improvements such as the removal of noisy pads, sector-by-sector calibration, and 
drift-distance dependence are expected to bring the resolution closer to this expected design 
value. These refinements are however not critical for the analysis presented here. 

This is the first stage of the track finding and fitting pertaining to the TPC only. The 
information is then passed to the global reconstruction software, discussed next. 



6.2 Global Event Reconstruction 



There are two additional steps to arrive at the final tracks used in the analysis. The connec- 
tion of the space point information among different detectors is the next step in tracking, i.e. 
a global track fit. This step becomes indispensable when including SVT or any additional 
information from a tracking detector. Finally, after the primary vertex is found, this acts 
as an additional measured point. If a track's trajectory goes near the primary vertex, a 
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subsequent primary track fit is then performed. For the analysis presented here, as we focus 
on production of primary charged hadrons, primary tracks are used. 

The global event reconstruction offline software correlates all the tracking, timing and 
energy deposition information from each detector in STAR and produces the overall event 
characteristics of the triggered events into a Data Summary Tape format for use in all the 
physics analyses. We will discuss the major components here in roughly sequential form, 
although one sometimes performs several iterations among different modules. This is the 
case for example if one uses PID information in global track propagation, but would also 
want tracking information when calculating track segment lengths for the calculation of the 
specific ionization in a PID module. The flow of the global software chain is sketched in 



Fig. p. 4 The SVT is shown in the flgure, although it was not yet installed for the 2000 
run; we will therefore not discuss it here. The main objective of the global part of the chain 
however, is to combine the information from all STAR sub-detectors and was therefore 
exercised with simulated data to ensure readiness for real data. 




Figure 6.4: Schematic of the STAR Offline Software chain. 
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6.2.1 Global Tracking 

The task of a global tracking module is two-fold. At its core, it is in charge of matching 
track and hit information from the different detector sub-systems. Once this is done, it must 
refit the newly matched tracks taking into account all available information. The role of 
the global tracking software was somewhat special for the 2000 run, since only 1 detector 
subsystem provided track information. In principle, one could then simply use the same 
TPC tracking parameters. However, since this would only be the case once, it was decided 
to actually perform a separate fitting routine at this level in order to work in a mode similar 
to that of having several sources of space point information. 

The global tracking then really performed an independent track fit from the TPC. It 
used a 3D helix model and incorporated energy loss and multiple scattering in the TPC gas 
when doing the track propagation. The errors assigned to the space points were given by a 
parameterization, based on real data, of the hit residuals as a function of the track crossing 
angle with respect to the pad-row direction. This is also a case where one needs an iteration 
in order to have information flowing both ways: the tracking needs the errors in order to 
weigh each space point for a minimization routine to obtain the best track parameters, yet 
the hit errors depend on the crossing angle which is known once the track parameters are 
set. 

6.2.2 Primary Vertex Finding 

The offline reconstruction finds a primary vertex for each event by propagating the global 
tracks through the field towards some reference point which is close to the estimated vertex. 
One can choose a start point towards which to propagate the tracks either in the transverse 
plane or along the beam line. Since the interaction point is much better determined in 
the transverse plane (RMS ~ 0.5 mm) than along the beam axis (RMS ~ 90 cm). The 
method at the heart of the routine is a minimization of the perpendicular distances from 
the track vectors to a point. We can see the basis for this method in a simplified scenario 



69 



where all the tracks are straight. In the transverse plane, the real vertex can only be a 
few millimeters away from the central axis. After propagating the helices to their point 
of closest approach to the reference point, we can then approximate them as straight lines 
around this point as in a Taylor series. This has the advantage of leading to an analytical 
solution. We consider then an arbitrary vector Vi, with a unit vector ey. . The distance from 
a point P = {xQ.yQ.zo) (in cartesian coordinates) to the vector Vi is the norm of the vector 
obtained through the cross product of the unit vector ey. with any line connecting the point 
P to any point Q along the line of the vector Vi, Q = {xi,yi,Zi). The line segment P(^ is 
{xq — Xi)x + (yo — yi)y + {zq — Zi)z, and the distance from P to is 

\di\ = \P^x Vi\ . (6.2) 

A least squares minimization must employ a merit function (x^) to assess the quality of the 
fit. We take the partial derivatives of the merit function with respect to the fit variables, 
and set them equal to zero in order to solve for the optimum values of the fit variables. For 
the vertex finder, the merit function is the summation of the squares of the perpendicular 
distances from the vectors to a point {xo,yo, zq), with the appropriate weight given by 
the uncertainty in the propagation of the track, 

By taking the partial derivatives with respect to xq, yo and zq and equating them to zero, we 
are left with a system of three equations for three unknowns, the coordinates of the vertex 

The vertex finder then has two main components. One algorithm is used to extrapolate 
a helix and calculates the coordinates of the point of closest approach from a given point in 
the hend plane, defined relative to the B field which is along the beam axis. The transverse 
plane is defined perpendicular to the beam axis. The other task is to perform the least 
square fit to find the common vertex. Additional iterations can be performed such that the 
linearization is done as close to the actual vertex as possible, and rejecting outlier tracks in 
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the fit. The details of the vertex finding algorithm are found in Ref. | 141 | 

The vertex resolution for high multiplicity events is approximately 150 ^m, both perpen- 
dicular and parallel to the beam axis. The vertex finding efficiency is 100% for events with 
more than 50 primary tracks in the TPC acceptance, decreasing to 60% for those with fewer 



than 5 primary tracks. We discuss the vertex finding efficiency in more detail in Section g.1.3 
as it is an important part of the analysis of the multiplicity distribution. 

6.2.3 Primary Track Fit 

Once the vertex has been found, global tracks whose distance of closest approach to the 
vertex is less than 3 cm are chosen for a refit using the vertex as an additional space point, 
yielding as a result different parameters associated with primary tracks. For high multiplicity 
events, the error associated to the primary vertex is much smaller than the error associated 
to TPC space points, so even though we are only adding one additional point to the track 
it can significantly improve the momentum resolution for tracks that actually come from 
the vertex. This is the reason for choosing the primary track parameters for the analysis 
presented here. The track model in the STAR geometry and magnetic field can be to first 
order represented by a helix, which can be parameterized as 

X = xq + R{cos{^{s)) - cos{^)) 

y = yo + -R(sin($(s)) -sin($)) (6.4) 

where R is the radius of curvature, $ = -(/; + 7r/2 and <I>(s) = <I> — {s/R) * cos(A), and s is the 
path length of the helix (see Fig. |6.5| ). The radius of curvature is related to the transverse 
momentum through 

R = q/c (6.5) 

c = i^— (6.6) 

P± 

where q is the charge of the particle, c is the curvature, and B is the magnetic field. With 
the following units: B in KGauss, p± in GeV/c, and q in units of the proton charge; the 
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(a) 



(b) 



Figure 6.5: Projection of a helix onto the transverse plane (a) and onto the bend plane (b) 

value of K to yield a curvature c in cm is K = 0.000299. The parameters that define the 
helix (xo, Ho, ip, A, q, c) are obtained by two simultaneous 2D fits. We separate the fits 



in the bend plane and in the transverse plane, shown in Fig j6.q . In the transverse plane, 
the space-points of a helix will project onto a circle, so the fit in the transverse plane {x 
and y coordinates in the STAR global coordinate system) is to a circle. In particular, the 
curvature of the track is obtained from such a fit, which relates to the transverse momentum 



as in Eq. |6.6| . With the knowledge of the curvature, we can now focus on the fit in the 
bend plane. We can think of the helix as a string winding around a cylinder. If we cut the 
cylinder and unfold it into a plane, the helix segments will be mapped into straight lines in 
this plane. We then perform a straight line fit in the bend plane z = (tanA)s + zq where 
the coordinates for the fit are the path length (s) of the helix and the beam axis coordinates 
(z). The parameter A is called the dip angle in the helix parameterization, and it is identical 



to the polar angle 9 in the definition of pseudorapidity, Eq. [4.17 . 

It is possible that secondary and decay tracks that are fit with the constraint of going 
through the primary vertex will have worse track parameters than those obtained from the 
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global track fit. This is especially important for high-p_L tracks, since biasing the tracks 
towards the primary vertex can artificially increase or decrease the curvature. For the 
majority of the tracks however, and in particular for true primary tracks, the resolution is 
improved by the primary track fit. Tracks for which the primary fit fails are flagged and not 
used in the analysis. 

Momentum Resolution 



The p± momentum resolution for vr for various multiplicities is shown in Fig. 6.6. This is 
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Figure 6.6: tt~ momentum resolution 



obtained from embedding pions into real events using the TPC simulation discussed in the 
next section. The effects of the momentum resolution in the measurement of the p± spectra 



is discussed in Section 3.2.5 
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6.2.4 Vo's, S's, and Kinks 

A specific part of the main STAR reconstruction code is dedicated to the topological iden- 
tification of neutral and charged particle decays. Neutral particles such as the and A 
can be reconstructed by identifying the secondary vertex, commonly called Vq vertex, of 
their charged daughter decay modes. Kg vr+vr^ and A —>■ p + ■k~ . The more compli- 
cated "cascade" decay schema of the H and O baryons can also be reconstructed by a series 
of topological pattern recognitions using the A's found in the first iteration through the 
decays: H~ — > A + vr" and Q~ — > A + . 

The charged kaon and pion one-prong decays found in the TPC volume are also recon- 
structed at this level {K fi + u). The topology signature of this decay is that of a track 
which decays in flight while propagating in the gas of the TPC, with the daughter charged 
track following a slightly different direction, and is called a kink decay. 

These steps are important for finding strange hadrons, i.e. A and H candidates using the 
Vo-finder algorithm and charged kaons using the kink-finder algorithm. In particular, the 
identification of charged kaons through the kink decay technique extends the available phase 
space to high momenta, where the ^ measurement no longer has any resolving power. For 



recent results from the analysis of strange hadrons in STAR, see e.g. |9C]. 

The secondary vertices are then reconstructed using the global tracks, propagating them 
towards the primary vertex and requiring that the Vq vertex of the track pair should not 
be close to the primary vertex. In addition, when combining the pair momenta to obtain 
the momentum of the parent, one requires that the parent should point back to the primary 
vertex. This is the last important part of the global reconstruction. Additional pieces of the 
code mainly organize the reconstructed information into the STAR Data Summary Tape 
(DST) to be written to disk for further analysis. Although the main tuning of cuts to find 
secondary vertices is done at the analysis level, the initial stages of Vq finding are done at 
the reconstruction level since part of the input are the space points, which are not saved to 
disk to minimize data volume. We turn our attention now to the TPC simulation. 
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Chapter 7 



Detector Simulation 



This analysis relies heavily on a realistic simulation of the response of the TPC. An esti- 
mation of the backgrounds and tracking efficiency is best done through such a simulation, 
and it is therefore necessary to include a significant amount of physics details. The primary 



function of the TPC Response Simulator (TRS) |142| , |143[| is to reproduce the electronic 
signals produced by tracks ionizing the gas of the TPC. The physical processes that are 
implemented are the drift of the ionized electrons in the gas, the amplification of the signal 
in the sense wires, the induction on the readout pads, and the response of the readout elec- 
tronics which produce digitized data. While the software was designed to cope with varying 
levels of detail in each of the different simulation stages, we describe here the implementation 
that was used for the analysis stage in the past year. 

7.1 Detector Response Overview 

The algorithmic basis for the response simulations is that developed in detail by the ALEPH 



TPC group, described in Ref. ||133|| , the standard reference for drift chambers. The TRS 
package contains parameterizations of the TPC response at various levels of detail. 

As part of the requirements, TRS was designed to deal with two sources of ionization 
of the TPC gas: charged tracks and laser events. The charged tracks create clusters whose 
distribution of amplitudes {i.e. number of electrons in a cluster) is Landau-like, whereas 
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the laser creates uniformly distributed ionization, resulting in Poisson distributed cluster 
amplitudes. 

For the simulations used in the analysis, the magnetic field was taken to be uniform with 
a value of 0.25 T in the z direction, and parallel to the electric field. During the analysis 
of the year 2000 data, several corrections needed to be applied to the raw data to take 
into account various distortion effects and modifications to the idealized geometry. These 
modifications were necessary for analyses such as the reconstruction of resonance decays and 
tracking of high particles. However, for the analysis of low momentum charged tracks 
up to 2 GeV/c, these effects were negligible. The analysis was repeated every time new 
distortion and geometry corrections were added to the reconstruction chain and the results 
did not vary significantly. For example, a comparison of the measured raw multiplicity in 
the midrapidity region before and after distortion corrections using different event samples 
showed that they were consistent to better than 1 %. Therefore, it was sufHcient to do the 
simulations in this idealized scenario. 

There are four main processes that are essential for the simulation to reproduce the 
behaviour of a TPC. 

• Ionization Transport - charge transport of the ionization and its deposition in the 
active region of the readout chambers. 

• Charge Collection - electron/ionization collection on the sense wires of the multi-wire 
proportional chamber (MWPC). 

• Analog Signal Generation - charge induction on the pad plane and generation of the 
time evolution of the analog signals on the pads. 

• Digital Signal Generation - conversion of the analog signals into ADC counts. 
We now discuss briefly each of the processes. 
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7.2 Ionization Decomposition and Transport 

The input into the simulation of the TPC response is the charge ionization left by tracks 



in the sensitive volumes of the detector generated by GSTAR p44|| , the GEANT 
implementation of the STAR detector. The ionization transport takes the charge deposited 
in these active areas of the detector and transports it through the field cage structure 
to the read-out plane. The GEANT simulation provides the amount of energy deposited 
(dE) over a given path-length (ds) by a particle with momentum p. Given the average 
ionization potential of the gas, one can calculate the total number of electrons such that 
the transport can be done at the segment level (dE) or to break each segment into parts 
with the limiting case of transporting each ionization electron individually. This provides 
a mechanism which allows the possibility to distinguish between a detailed microscopic 
simulation and a macroscopic parameterization. In addition, one can also vary the length 
of the segment that is transported (and subsequently processed). Thus, the granularity of 
the simulation can be adjusted to provide different levels of detail. 

The ionization is distributed on the pad plane according to distributions which char- 
acterize the effects of diffusion. The role of the charge transporter is to alter the x and y 
positions according to transverse diffusions, to alter the z position to reflect the drift time 
with longitudinal diffusion folded in, and to alter the amount of charge that reaches the 
read-out plane. 

7.3 Charge Collection and Amplification 

Once the electrons arrive at the read-out plane, they must be collected by the individual 
anode/sense wires of the multi-wire proportional chamber (MWPC). It is here that the 
avalanche process multiplies the signal of several tens of electrons to several 10^ — 10^ 
electrons, depending on the potential of the wires. 

The basic principle of avalanche creation in a proportional wire is as follows. As an 
electron drifts toward the wire it travels in an increasing electric field. In the vicinity of the 
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wire, the field E at radius r is given by the hnear charge density A on the wire 

E = ^^-. (7.1) 

27reo r ^ ' 

The electron's path is directed towards the wire. The trajectory also terminates in the wire 
in the presence of a magnetic field if the electric field is sufficiently strong. Once the electric 
field surrounding the electron is strong enough, the electron will have enough energy to 
ionize the gas and still pick up sufficient energy between collisions to ionize the gas further, 
starting an avalanche. The number of electrons continues to grow until all electrons reach 
the wire. The process develops over as many mean-free paths as there are generations of 
electrons, typically 50-100 ;um, and takes less than 1 ns for most gases. Since the signal is 
proportional to the number of electrons collected, the name "proportional wire" was adopted. 



For a review see ||133|1 and |146] which also contains a list of texts. 



The distribution of avalanches started by single electrons can be modelled as a stochastic 



process (see |133|, p. 144). A simple description is given by the Yule-Furry process, after the 
authors who first used it to describe biological population growth. The idea was applied 
to proportional avalanches by Snyder and Prisch (see |147| for references), and is briefly 
described in what follows. 

We start from the idea that as the avalanche (or other stochastic process) develops, the 
variation of the number n of ions (or electrons) is followed as a function of some parameter 
t. This parameter can be thought of as the time, but it is a bit more general in the sense 
that we need only require that its increase describes the progress of the process in the order 
given by the causal sequence of events. Therefore, t need only be a monotonic function of 
time. We are interested in the probability P{n, t) that at some 'time' t there are n electrons 
in the avalanche. The initial condition is that the avalanche starts with a single electron, 
i.e. 

P(1,0) = 1. (7.2) 

We suppose that the probability for the birth of one electron in any interval At is 
proportional to the number of electrons n with some proportionality constant A. The solution 
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to the corresponding differential equation giving the probabihty distribution at a given t is 

P(n,t) = e-^*(l -e-^*f-^ (7.3) 
The mean n and the variance cr^ are 

^nP(n,t) = e^* (7.4) 



CO 

n = 

n=l 



a^ = J2 n^P{n, t) - = e^*(e^* - 1) . (7.5) 



n=l 



Rewriting 7.3 with 7.4 the distribution becomes 

1 / 

P{n) = ^ 1 - ^ , (7.6) 

n \ n J 

where the 'time' variable no longer appears. In the limit n oo, which is appropriate for 
avalanches, we obtain 



P{n) = 4e-"/" , (7.7) 



n 



a' = fi' . (7.8) 

We see that the Yule-Furry process has an exponential signal distribution. Subsequent 
theoretical refinements to this simple expression which were made in order to take into ac- 
count effects like the asymmetric growth of the avalanche profile as well as saturation effects. 
One such refinement is the Byrne process, which introduces the idea that a fluctuation to 
larger n in the first part of the avalanche reduces the rate of development in the second 
part (so an additional parameter has to be introduced, whose effect is to suppress the small 
amplification factors). The Byrne process includes the Yule- Furry process as a special case 
when the additional parameter approaches zero. The solution to the differential equations 
for this process is given by the Polya (or negative binomial) distribution function. 

Exponential behaviour at low to moderate gas gains [i.e. < 10^) is observed experi- 
mentally. In parallel plate geometry however, slight deformation from exponential shape is 
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observed at gas gains above 10^. This is probably due to self-saturation effects which become 
important in the presence of space charge. Since the STAR TPC is generally operated at 
low gas gains, the simple exponential is found to be acceptable for the simulations. 

It is important to note that in any drift chamber operation, the effect of the fluctuations 
in gas gain is to simply degrade the attainable space-point resolution, and for this purpose 
it is not critical to use the exact functional form of the avalanche yields. The effect of 
degradation is the process that the simulation of the gas-gain fluctuation is attempting to 
reproduce, and the simple exponential should be sufficient in this regard. 

Once the charge has been amplified, we can proceed to calculate the amount of charge 
induced on the cathode pad plane. 

7.4 Analog Signal Generation 

This simulation step has three main parts: 

• Determination of the charge induced on single pads from the charge collected on the 
anode wires 

• Sampling of the induced charge signals in time according to the electronics response 
{i.e. pre-amplifier and shaper) 

• Distribution of the analog charge into time bins 
7.4.1 Charge Induction 

The charge induced on a grounded pad plane by a point charge q located a distance d above 
the plane can be calculated by the method of images. The charge density a on the plane is 
given by: 

^^^'^^ = 2^ ((x - + {y-yo)' + d^r/^ ^^-^^ 

where the charge is located at the position {xo,yo)- However, the typical geometry of a 
MWPC is a bit more elaborate. The proportional wires normally form a plane between two 
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grounded cathode planes. In addition, the charge is not point-like, but in the form of a line 
of charge with density A. We can simplify and locate it in the middle of the parallel plates. 
In order to calculate the charge density induced on the pad plane we thus have to include all 
higher-order multi-pole terms. The images are alternately negative and positive, situated 
at Zn = ±(2n -|- l)d, {k = 1,2, • • • ) and z^, = d. The total charge density is obtained by 
integrating over y and doing the sum 



which is called the Endo function. The derivation was done without taking into account 
the limited extent of the pads in the y direction, i.e. for the limiting case of zero width to 
length ratio w/l. 

The effect of finite geometry of segmented cathodes can be taken into account with the 
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addition of another parameter. We can rewrite Eq.7.10 



1 _ tanh^ ( 



where Ci is a normalization constant. Eq. |7.1l] can be used to generalize the Endo function 
to take into account finite geometry effects due to segmented cathode pads, i.e. w/l > 0. 
By introducing the constant C2 

1 _ tanh2 ( ^^i^) 

F{x) = Ci TV (7-12) 

l + Cotanh^^^^i^' 



we arrive at the generalized solution to the distribution of charge induced on a grounded 
pad plane, and is usually dubbed the Gatti function. 

A comparison of the Gatti and Endo functions, along with a Gaussian (for comparison 



purposes) is given in Figure 7.1. For reference, the functions were normalized to 1 at x = 
for the plot and the parameters were a = 1.5 for the Gaussian, d = 2 for the Endo function 
and d = 2, C2 = 0.5 for the Gatti function. 
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Figure 7.1: Comparison of Gatti, Endo and Gaussian functions as profiles of the pad-response 
function. 

For the case of the TPC, the quantity of interest is the total amount of charge (Q) 
induced per pad, so we have to integrate a given pad-response function, which gives the 
charge density, over the area of a pad: 



The advantage of using such functions is that they allow the production of longer tails 
which have non-Gaussian characteristics. The tails are an important feature to understand, 
as they determine the efficiency of the ionization collection which is of relevance in the study 
of ^ resolution. As an additional case, the pad response function was measured for the 
STAR TPC pads in the directions transverse to and parallel to the anode wires. This allows 
us to calculate the pad response as a convolution of the transverse diffusion function with a 
parameterization of the measured pad response function. 




(7.13) 
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7.4.2 Sampling of Signal in Time and Shaper Response 

Once the amount and centroid of the charge distribution on each pad is determined, this 
charge can be sampled in time, thereby emulating the response of the analog electronics. 
The signals are generated by superimposing each analog signal from each of the avalanches 
that induce a signal on the pad plane. This permits one to vary the shaping time of the 
electronics independently of the width of the pad response function. In reality, the time 
evolution of the signal that develops on the wires is almost entirely due to the motion of 
positive ions away from the wire. This produces a signal with a long tail. For the STAR 
geometry, the duration of this signal is on the order of ^ 62 fis. In order to make the 
detector faster, the signal is differentiated after a characteristic time - the shaping time of 
the pre- amplifier. The trade-off is that only a fraction of the total charge is seen by the 
downstream electronics, i.e. the ADC. The fraction F of the charge is given by 



where tm is the length of time the long, undifferentiated signal would persist (~ 62 /xs), to is 
the characteristic time of the development of the signal 1 iis) and ts is the shaping time 
of the pre-amplifier (~ 180 ns). 

For STAR this means that the signal shape is dominated by the shaping properties 
of the electronics. Therefore, the long time response of the chamber is parameterized in 
the electronics processing section of the simulator, instead of modelling the motion of the 
positive ions. The shaper response function for the case of one stage differentiation and two 
stage integration used in STAR is 



where is the step function, and the time constant r is 55 ns for the STAR electronics. 
This function is then convoluted with a Gaussian to parameterize the effect of a longitudi- 
nally diffused cluster to yield the final response of the shaper. 



F = 




(7.14) 




(7.15) 
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Once we obtain the functional form of the signals induced on a pad over the read-out 
period of the TPC electronics, we can distribute the analog signal into discrete time bins. 
This sampling simulates the behaviour of the Switched Capacitor Array (SCA) in the front 
end electronics. Essentially, this entails integrating the amount of charge in a time interval 
At, which is determined by the SCA sampling frequency. Simulated chamber noise and 
electronic noise are also added at this point. 

7.4.3 Digital Signal Generation 

Once the analog charge is distributed into time bins on the pads, we can digitize the signal. 
This step is done via a simple conversion from voltage to ADC counts. This is the final stage 
in the TPC simulation, and the pixel values are then used as input to the reconstruction 
chain in exactly the same way as the raw data. 

7.5 Embedding 

One of the main purposes of the detector simulation response is to obtain estimates of 
tracking corrections to the raw particle yields. It is of great importance to have a reliable 
detector simulation as detailed in the previous section. To make the environment as realistic 
as possible, we use a procedure commonly known as embedding. The idea is to take a real 
event and embed into the raw data file the signal from a few simulated tracks at the level 
of ADC counts. The simulated signal is obtained via a GEANT Monte Carlo simulation of 
the energy deposited by the tracks in the detector volume followed by a simulation of the 
detector response as detailed in the previous section. The final simulated ADC signal is then 
convoluted with the raw data from the real event and then fed into the STAR reconstruction 
software chain. Since we have the information from the GEANT simulation, we can associate 
the reconstructed tracks to the corresponding Monte Carlo track. The association is done 
at the hit level, since from the reconstruction side the hits are the "seed" from which the 
track parameters are built. 
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Figure 7.2: TPC hit position resolution in the drift direction. 



The embedding procedure combines the advantages of having a very reahstic event en- 
vironment along with having a controlled track population. From the comparison of re- 
construction output to simulated input we can obtain information about many interesting 
performance diagnostics of the detector. 



Figure 7.2 shows the position resolution for hits in the TPC in the drift direction. The 
RMS width is found to be 180 ^m, a similar analysis in the padrow direction yields an 
RMS width on the order of 120 fim. Once we have a match between the TPC hits, we can 
use this information to make a matching between tracks. For every reconstructed track, we 
scan through its reconstructed space points and find their corresponding Monte Carlo hits. 
These simulated hits have the necessary information of the track that generated them, and 
so form the basis for the track match. It is this track match information that we use for 
the determination of the reconstruction efficiency. We discuss the corrections to the raw 
spectra obtained through the embedding procedure and through background simulations in 
the next chapter. 
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Chapter 8 



Analysis of Charged Hadron Spectra 



8.1 Event-wise Studies 

In this section we delineate the relevant parts of the analysis that was performed to obtain 
the multiplicity distribution. This is also the baseline for any other analysis. 

8.1.1 Trigger and event selection 

We use the ZDC and the CTB for triggering. All the runs that we used for this analysis were 
recorded with a minimum bias trigger, with no pre-scale. The minimum bias trigger consisted 
of a coincidence requirement between the East and West ZDC's. This trigger proved to 
be better than 99% efficient. At the beginning of the run, however, this had yet to be 
established. The main concern was that there might be an inefficiency in this trigger for the 
most central collisions, where very few spectator nucleons remain to produce a coincidence in 
the calorimeters. Therefore, a high CTB signal, indicative of a high multiplicity in the mid- 
rapidity region, was used as an additional trigger condition. The ZDC trigger thresholds 
were set such that a single nucleon hitting the calorimeter would generate an acceptable 
signal. The CTB threshold was set to accept the highest ~ 30% multiplicities. 

A plot of the ZDC vs. CTB trigger signals for a subset of the events used in this analysis is 



shown in Fig. 8.1. The peripheral collisions are in the lower left corner, where both the CTB 



signal and the sum of the East and West ZDC signals are small. This is indicative of low 
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Figure 8.1: The ZDC (sum of East and West calorimeters) signal is plotted vs. the corre- 
sponding CTB signal. The most central events have the highest CTB signal (high multi- 
plicity) and a low ZDC signal (very few spectator neutrons). 



multiplicity at mid-rapidity and a scarce number of dissociation neutrons. In the collider 
geometry, neither the excited nucleus nor the dissociation protons reach the calorimeters 
because their trajectories are bent by the beam optics. As the overlap of the colliding 
nuclei (i.e. the centrality) increases, more neutrons reach the calorimeters and the ZDC 
signal increases. Likewise, the multiplicity at mid-rapidity increases. After a certain point, 
the collision is sufficiently central that few neutrons reach the ZDC while the multiplicity 



continues to increase. The "boomerang" shape observed in Fig. |8JJ is therefore the result of 
the correlation between impact parameter b and multiplicity on the CTB side, and of a dual 
behaviour on the ZDC side: correlation between b and number of neutrons at high impact 
parameter, and anti-correlation for central collisions. 

An important input into this analysis is the hadronic cross section in Au -|- Au, ctauAu- 
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This value is necessary to normalize the multiplicity distribution. The calculated value is 



o"AuAu = 7.2 barn from Eq. |4^, but it is desirable to have a measurement confirming this 
expectation. The trigger requirement of a ZDC coincidence is not only sensitive to the 
hadronic cross section, but also to the mutual Coulomb dissociation of the two Au nuclei. 
In this process, there will be a pair of correlated forward- and backward-going neutrons, 
but no measurable tracks in the mid-rapidity region. The ZDC's therefore are sensitive to 
the sum of the hadronic + mutual Coulomb dissociation cross sections, which we denote by 
o"xn,xn, the cross section in which at least one neutron is detected in each ZDC detector. 
Since the same ZDC's are used in all interaction zones at RHIC, this combined cross section 
can be measured by all 4 RHIC experiments independently. The STAR preliminary result. 



measured by the van der Meer scan technique [148| is (Txn,xn = 8.9 ± O.Sgtat i 0.7sys barn. 
This value has also been measured in an independent analysis by the RHIC accelerator 



crew |149|. The value reported in a calculation [ |150[ is crxn,xn = 10.9 ib 0.6 barn, and 
seems to be outside of the allowed range of the systematic uncertainty in the measurements. 
Attempts to resolve this discrepancy are underway. In order to arrive at a measured hadronic 
cross section, ctauAu, it is also necessary to disentangle the contributions from the hadronic 
and the Coulomb processes. Since a pp inelastic collision at ^/s = 130 GeV produces on 



average ~ 2.4 (see Eq. |4.1l|) charged particles per unit pseudorapidity at mid-rapidity and we 
expect to identify the vertex for hadronic events, to first order this fraction can be measured 
as 

o"AuAu Events with Vertex 
o'xn,xn Total Triggered Events 

This simple ratio must be corrected for the vertex-finding efficiency and acceptance that 
reduce the counts in the numerator and for backgrounds such as beam+gas that generate 
a trigger and are therefore counted in the denominator. The vertex finding efficiency is 
important in this analysis and for the determination of the shape of the h~ multiplicity 
distribution at low multiplicity. 
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8.1.2 Vertex Acceptance 

The vertex reconstruction for the summer 2000 run was undertaken essentially up to the 
limits of the TPC acceptance, i.e. I^vcrtexl < 2 m where -^vertex denotes the position of the 
primary vertex along the z direction parallel to the beam. Although it is possible to try to 
obtain a vertex as long as there are tracks in the TPC, this was not pursued for the following 
reasons. If an event occurred outside of this region, the TPC sees only the forward- or the 
backward-going particles. This also has the disadvantage of reducing the accuracy of the 
vertex determination, since the z position of tracks going perpendicular to the beam axis 
contribute significantly to the determination of the vertex. In addition, the CTB in this 
case is not triggered on multiplicity at mid-rapidity, but rather on multiplicity forward or 
backward of mid-rapidity with all tracks having a large dip angle. To avoid these systematic 
effects, the vertex finder was set to abort when it determined that the vertex lay outside the 
TPC limits. 

Due to the large size of the beam diamond, the standard deviation of the vertex Z 
position was ~ 100 cm. This results in valid collision events, even at high multiplicity, 
without a vertex determination. In addition, in order to keep a flat acceptance for tracks 
with Ir^l < 0.5, we restricted the multiplicity analysis to events with a vertex between ±95 
cm, about a ±la cut. To correct for the vertex acceptance, we assume that the distribution 
is Gaussian and from the fit parameters estimate the fraction of events that lie outside our 
acceptance cut. As an example, the ^vertex distribution for lOK triggered events from one of 
the runs taken in September 2000 is shown in Fig. 

For the distribution in Fig. ^.21 , the acceptance correction factor is 1.57. Analyzing more 
runs for a total of 166K triggered events, we find an acceptance correction factor of 1.68, (a 
±la cut for a Gaussian distribution would yield a correction of 1.46, so our cut is tighter). 
With the vertex efficiency and acceptance correction, we can then determine the numerator 



in Eq. |8.1| . For the denominator, we must estimate the backgrounds to our minimum-bias 
trigger. 
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Figure 8.2: Distribution of reconstructed vertex z positions, ^vertex, for 10,000 triggered 
events. The limits of the TPC are at ±2 m from the center. The acceptance for events used 
in the multiplicity analysis was placed at ±95 cm. 



8.1.3 Vertex Efficiency 

We know that our offline vertex reconstruction is not 100% efficient, therefore we need to 
correct for this effect when measuring the multiplicity distribution. The vertex reconstruc- 
tion is based on tracing the path of reconstructed tracks back to a common space point. 
Since tracking is done first, before finding the vertex there is no a priori knowledge of which 
tracks are primary and which ones are not. The vertex-finding therefore depends on all 
global tracks (Sec. 6.2.1] ) found in the event. We will then characterize the efficiency as a 
function of the number of global tracks (A'giobai) in the event. The vertex-finding efficiency 
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(cvtx) affects mainly the low multiplicity events. The reason is that these events provide 
the vertex- finding algorithm with very little information to work with, i.e. very few global 
tracks. We find that for events with more than ~ 100 global tracks, the efficiency is ~ 100%. 
To confirm the findings based on software, a visual analysis of ~ 100 events was performed, 
with a resulting lower bound on the efficiency for events with A'^giobai = 100 of evtx > 98%. 

The more important part is the efficiency at low A'giobai- To get a handle on this number, 
we used hijing events generated with large impact parameter 6 = 12 - 20 fm. At these impact 
parameters, hijing should be a reasonable model, since basically only geometry and the pp 
cross section play a role and nuclear specific effects are not expected to influence the results. 
These events were then processed through the STAR reconstruction chain and a record was 
kept of the number of events generated and whether the vertex was or was not found for 
each event. In this way, we obtain a plot of the vertex efficiency correction factor (1/evtx)- 
The correction depends on A'giobai and is illustrated in Fig. 
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Figure 8.3: Vertex finding efficiency correction factor as a function of the number of global 
tracks, iVgiobai- 



To take this correction into account in the h multiplicity distribution, the correction 
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in Fig. ^ was used as a weight for each event depending on the number of global tracks. 
Since the h~ multiplicity distribution refers to negative primary tracks instead of global 
tracks, this correction affects primarily the region A'^^- < 5 and is negligible beyond Nf^- 
= 10. These are the first two bins of the final distribution. However, most of the cross 
section is found in these bins, and the shape of the distribution is sensitive to these values 
at low multiplicity. To figure out what is the fraction of the hadronic cross section that we 
actually see, we proceed as follows. We choose events within a ±95 cm range along the z 
axis of the beam direction. We can then count the the raw number of events with vertex in 
our multiplicity distribution. We compare this to the number of vertex events corrected for 
efficiency, i.e. to the number of events in the multiplicity distribution appropriately weighted 



with the correction from Figure 8^. From these two numbers, we see that the fraction of 
the hadronic cross section available to the offline analysis is 94.9 ± 0.5 ± 4 % of ctauAu, 
where the systematic uncertainty comes from the uncertainty in the estimation of the vertex 
efhciency at the lowest multiplicities based on the simulations. 

8.1.4 Trigger backgrounds 



We now focus on estimating the contribution to the events in the denominator in Eq. |8J_ 
that are not part of the total crxn,xn cross section, e.g. beam-gas events. However, knowing 
which events are background so that we can take them out of our total trigger sample, is 
not trivial. 

To investigate this issue, we focus on events for which the vertex was not found. To study 
systematic effects and variations, from the September data, we used specific minimum bias 
runs with good statistics, and stable detector operation. The run with the most statistics 



in this period has ~ 77K events (the events in Fig. is from this sample). This was used 
along with a few other runs to compare backgrounds and other systematic effects. 

For these background studies, we compared 4 runs with different luminosity conditions 
spread over several weeks. Two of these (August 2000 data) had intensities in RHIC which 
were about a factor of 3 lower than in the other two (September 2000) runs. We are interested 
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in the different intensities because the hadronic interaction rate grows as the product of the 
intensity in each beam, while the beam+gas background grows only as the sum. This type of 
background should then have a different contribution in runs with differing beam intensities. 
To illustrate the difference between the events with and without vertex, a plot of the 



trigger information from lOK events from one of the runs is shown in Fig. 8.4 Events with 
a reconstructed vertex are shown in the left panel; the right panel illustrates the trigger 
signals for events without a reconstructed vertex. 
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Figure 8.4: The trigger signal for events with a valid reconstructed vertex (a) shows the 
expected shape. There are additional event classes that show up when plotting the trigger 
signal for events without a reconstructed vertex (b). 



For events in Fig. |8.4| (a), we see again the very clean "boomerang" band. These represent 
valid hadronic events for the determination of ctauAu- For the events without a reconstructed 
vertex, there are several regions to understand. The events which lie in the "boomerang" are 
most likely good events that did not have a vertex because the collision occurred outside the 
bounds of the TPC, where the acceptance of the offline vertex finder terminates. This does 
not present a problem. However, in the 2001 run a cross check will be done using timing in 
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the ZDC's. There are two other regions that are most probably background and have been 
the subject of further scrutiny: 

• ZDC < 30 (low ZDC region) 

• ZDC ~ 80 (mid ZDC region) 

Both of these regions have events whose trigger signals extend along the CTB axis. If 
the multiplicity is increasing for these events, we should see also an increase in the number 
of global tracks. The same thing is found if we plot the ZDC signal vs. A'giobaifor the events 



with and without vertex, shown in Fig. B.5 
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Figure 8.5: The ZDC signal vs. the number of global tracks for events with a valid recon- 



structed vertex (a) and without one (b). The same structure as in Fig. |8^is seen, indicating 
that there is definitely an increased number of charged tracks in the detector for the events 
in the two regions outside the "boomerang". 



We see that there are events that have ~ 1000 global tracks, do not have a reconstructed 
vertex and lie in a different region than the normal hadronic events. 

We studied the characteristics of these two bands to try to understand their origin. We 
investigated other differences between these bands and the rest of the "good" triggers to 
shed light on their nature. We expect that in the region A'giobai^ 1 most of the events are 
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of electromagnetic origin. In any case, there can be no vertex reconstructed for events with 
less than 2 tracks. By selecting non vertex events with A'giobai>l we focus on the possible 
backgrounds the we're aiming to understand. 

We find that the band at low ZDC has both East and West ZDC triggers just above 
threshold, and the band in the mid ZDC region has one of the calorimeters just above thresh- 
old and the other one high. For the ZDC ~ 80 region, this is just the opposite behaviour 
than is seen in the events with a found vertex, where the ZDC signal is relatively symmetric 
between East and West ZDC's. It is possible that a beam+gas event could generate such an 
asymmetric ZDC signal with some charged tracks seen in the TPC. Asymmetries in tracking 
were also found, but this can also indicate simply an interaction outside the TPC volume. 
Several hypotheses were proposed, but there was insufficient information to unambiguously 
discern the nature of the background. In particular, we would like to know its contribution 
to the region where the valid events lay. Without this knowledge, we can remove only the 
background contribution outside the "boomerang" region. We can therefore only produce a 



lower limit on the fraction in Eq. 8.1. The numerator depends on the efficiency and accep- 
tance corrections already discussed. For the purpose of the calculation of the fraction, we 
used events within ±2 m to reduce the acceptance correction (6%). For the denominator, we 
start with the total triggers. The background is obtained via the number of events without 
vertex that have 1 or more global tracks. Since some of these events are valid events that 
were not found because of efficiency or acceptance, we need to subtract the number of events 
we added to our numerator due to these two corrections from the background estimate. This 
is a lower bound since we don't take into account the background for events without global 
tracks. From lOK events from the low luminosity runs and from 15K events from the high 
luminosity runs we obtain 



Low£ : ^^>70.0±0.7% 

^xn,xn 

High C : > 73.5 ± 0.6 % . (8.2) 
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The errors are statistical, which is enough to compare the two numbers since they were 
obtained in the same way. We can see that there is a difference beyond statistical. However, 
this was expected as the low-luminosity runs systematically have a greater background 
contribution that our procedure does not account for. For the high-luminosity sample, this 
background is reduced and the lower limit on the fraction is higher. Using the fraction from 
the high-luminosity events and the current measured value for axn,xn = 8.9 ± 0.5 barn, we 
obtain a lower limit on the hadronic cross section ctauAu ^ 6.5 barn, consistent with the 
calculated value of 7.2 barn. There should be an error of ±0.4 barn in this estimate which 
comes basically from the cross section measurement (the statistical error on the fraction is 
much smaller). This error in turn is dominated by the uncertainty in the measurement of the 
beam currents in the RHIC ring, which is on the order of 3%. The systematic uncertainty 
on the fraction is the only missing piece, but as discussed, at this point we can only give a 
lower bound. 

It is clear that with only the ZDC and CTB information at the lowest multiplicities, it 
becomes increasingly difficult to separate the background events such as beam+gas collisions. 
It would therefore be of great advantage to obtain an estimate of the interaction point 
without having to rely on tracking, but rather on trigger signals, e.g. on timing between 
the arrival of the East and West ZDC signal. In addition, improved phase space coverage 
for the trigger would drastically reduce uncertainties in the background estimates and help 
to provide a more complete topology of the valid low multiplicity hadronic events. For the 
2001 run ZDC timing will be implemented; and starting with the 2001 pp run, additional 
detectors in the form of scintillator slats covering from the pseudorapidity region 2 < < 4 
will be part of the STAR trigger as well. For the analysis of the multiplicity distribution 



presented in Sec. |9JJ we will use the calculated value of the Au -|- Au hadronic cross section 
of CAuAu = 7.2 barn for the normalization. 

To obtain the multiplicity distribution, one has to obtain the corrected number of nega- 
tive hadrons for every event. The analysis done here relies on tracking. We first obtain a raw 
dN/dp±dr] distribution as a function of r] and p± for every event. Then several corrections 
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are applied to this raw distribution as discussed in the next section. Finally, the corrected 
dN/dp±d'r] distribution is integrated in the range \rj\ < 0.5, 0.1 < p± < 2 GeV/c to obtain 
a corrected multiplicity. One can therefore also obtain p± and ri distributions with this 
procedure. Since the algorithm relies on tracking corrections, we now turn our attention to 
this subject. 

8.2 Tracking Studies 

Particle production was studied through the yield of primary negative hadrons, which are 
mostly 7r~ with an admixture of and p. The distribution includes the products of 
strong and electromagnetic decays. Negatively charged hadrons were the main focus of the 
work in order to exclude effects due to participant nucleons which would show up in the 
positively charged hadron sample. Charged particle tracks reconstructed in the TPC were 
accepted for this analysis if they fulfilled requirements on number of points on the track and 
on pointing accuracy to the event vertex. 

The main goal of this analysis is to determine the corrected yield of primary particles 
coming from the collision. What we have in the final state are measured tracks, both pri- 
mary and secondary, in our detectors. There are several losses and backgrounds that need to 
be corrected: acceptance, decay losses, track reconstruction efficiency, contamination due to 
interactions in material, misidentified non-hadrons, and the products of weak decays. The 
appropriate corrections were obtained mainly by use of the embedding technique. Back- 
grounds were determined by either fully simulated events or by direct measurement when 
possible. All corrections were calculated as a function of the uncorrected event multiplicity. 
The corrections used in the determination of the final spectra are obtained in the following 
order: 

® Geometrical acceptance and decay losses 

® Reconstruction efficiency 

(D Track merging 

® Momentum resolution 
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® Track splitting 

® Electron background 

® Weak decay and secondary interaction background 
® Ghost tracks 

We now discuss the most important corrections in more detail. 
8.2.1 Acceptance 

The acceptance correction takes care of two different things. What we really focus on here 
is whether or not a charged particle makes a measurable signal in the detector such that it is 
possible to reconstruct it afterwards. If the charged particle does not deposit energy in any 
sensitive volume of the detector, it will be lost. This is the more common definition of the 
geometrical acceptance. In addition, a charged particle may also leave no signal if it decays 
in flight before it reaches the detector. In the STAR Monte Carlo simulation implemented for 
this work, these two losses were taken into account simultaneously by adopting the following 
definition of acceptance: a track is accepted if it leaves at least 10 Monte Carlo hits in the 
TPC. The reconstruction code then at least has the possibility to find the track. This 
correction can be calculated in a full Monte Carlo simulation and also using embedding, and 
should be completely independent of the event multiplicity. There should be a dependence 
on particle type, as the decay characteristics are included in our definition. This is illustrated 



in Fig. 8^. The left panel is the acceptance for negative pions in three multiplicity bins. 
The numbers correspond to the mean raw negatively charged multiplicity ((/i~)raw) in each 
of the 3 bins: (/i^)raw = 223 for the high multiplicity, (/i~)raw = 152 for medium, and 
{h~)ra.w = 33 for the low multiplicity bin. The plot shows that the correction is essentially 
independent of the multiplicity, except perhaps at the lowest p±. The right panel is the 
acceptance for the three particle species which make up the distribution. We see that 
the K~ acceptance is lower than that for the pions, which is expected since cr = 3.7 m for 
kaons, and it is 7.8 m for 7r~. 
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Figure 8.6: The acceptance correction for the h analysis. Left panel: vr acceptance for 3 
multiplicity selections. Right panel: Acceptance for vr" K~ and p. 



The acceptance is on average 90% for tracks within the fiducial volume having p± > 200 
MeV/c. We see that the acceptance rapidly drops at low pj_. For the h~ analysis, we accept 
tracks which have 0.1 < p± < 2 GeV/c and |?7| < 1.0. This accounts for the majority of the 
produced particles at mid-rapidity, as we find the yield beyond p± = 2 GeV/c to be only 1% 
of the total yield. A motivation to go as low in p± as possible was driven by the fact that 
for 2001 the magnetic field would be set to the design operating value of 0.5 T for the bulk 
of the data taking, instead of the 0.25 T used in 2000. This raises the low-p_L acceptance of 
the TPC, so essentially the 2000 data would give us access to the lowest p±. Nevertheless, 
to study systematic effects the multiplicity analysis was carried out three times with various 
low-p_L cutoffs: p±>0 (no cutoff, large systematics expected), p±>0.1 GeV/c and p±>0.2 
GeV/c yielding consistent results. 
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8.2.2 Tracking Efficiency 

Once a track has made it into the detector, we focus on the question of how hkely it is 
to be found by the offline software chain. The reconstruction efiiciency was determined by 
embedding simulated tracks into real events at the raw data level, reconstructing the full 
events, and comparing the simulated input to the reconstructed output. This technique 
requires a precise simulation of isolated single tracks, achieved by a detailed simulation of 
the STAR apparatus based on GEANT and a microscopic simulation of the TPC response 
discussed in Chapter 0. The multiplicity of the embedded tracks was limited to 5% of 
the multiplicity of the real event in the same phase space as the simulated data, thereby 
perturbing the real event at a level below the statistical fluctuations within the event sample. 
Fig. ^]7| shows the reconstruction efflciency obtained from embedding as a function of p± for 
a slice at mid-rapidity. The tracking efflciency varies depending on p± and the multiplicity 
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Figure 8.7: The reconstruction efflciency for the h analysis. Left panel: vr efficiency for 3 
multiplicity selections. Right panel: Efficiency for vr", K~ and p at high multiplicity. 



of the event. For the lowest multiplicity events, the vr efflciency is > 95% in the region 
p± > 400 MeV/c. Going to the high multiplicity events degrades the efflciency to ~ 80%. 
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However, since the difference between the multiphcities in the high and low bins is about 
a factor of 7, a reduction of only ~ 15% in the efficiency is a significant achievement for 
the offline reconstruction. In general, for all multiplicities and all 3 particle species, in the 
region p± > 200 MeV/c we always have an efficiency greater than 80%. 



The pseudorapidity dependence of the tracking efficiency for tt is shown in Fig. 8.8. The 
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Figure 8.8: The vr" reconstruction efficiency as a function of pseudorapidity. Two different 
vertex selections are shown, illustrating the influence of the choice of vertex on the efficiency 
due to the detector geometry. 
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two data sets in the figure are for selections of events with vertex z position, ^vertex, within 
±95 cm (blue points) and ±20 cm (red points) of the center of the TPC. The wide vertex 
cut allows us to use most of the available data. The interest in studying such variations 
in efficiency lies in the observation that a particular phase space region, for a given vertex 
position, will be measured in a different region of the detector. For example, tracks with r] 
> 0.5 will not cross all the available padrows in the TPC when the vertex is in a position 
-^vertex > 95 Cm. This is reflected as a drop in the efficiency. By making a tighter 2;vertex cut. 
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one can probe a wider phase space region with the guarantee that tracks have the chance 
to cross all TPC padrows. For a vertex cut of l^^vertexl < 20 cm, tracks within \r]\ < 0.8 will 
cross the entire TPC, and one can see a higher efficiency for these tracks in the figure (red 
points). 

Furthermore, even when tracks cross the entire TPC, we also expect a variation in 
efficiency for a given rj slice as a function of the vertex position. This arises from diffusion 
effects, as a track crossing the detector in the region near the central membrane will produce 
space points that must drift across the entire chamber, and their clusters will be wider than 
those for tracks near the pad plane. We can test the adequacy of our detector simulation by 
studying the raw and corrected yield in a given r] bin as a function of the vertex position, 
and hence as a function of the position of the track in the detector. In addition, since the 
TPC volume is separated into two identical halves separated at z = and set up such that 
the ionization drifts toward the wires (located at ±2.1 m), the behaviour for a given choice 
rj and -Zvertex should be the same as for —rj and — -Zvertex- This is illustrated in Fig. |8]^ which 
shows the raw r/ distribution for 10 cm slices in the choice of -Zvertex interval: [20, 30] cm, [80, 
90] cm; and the corresponding symmetric z^ertex 

selections [-30, -20] cm and [-90, -80] 
cm. The top panels show the raw distributions as a function of r/, left for the [20, 30] cm 
interval and right for the [80, 90] cm region. In the bottom panel a sketch of the detector 
geometry is given, as a slice in the r and z coordinates of the TPC cylinder. The central 
membrane is at z = 0, the drift distance is 2.1 m from the membrane to the wires. The 
first pad row at 60 cm and the last pad row at 200 cm. The lines in this panel represent 
the trajectories of tracks originating from the given vertex position at rj = ±0.5. We can 
see that there are significant systematic variations in the raw yield depending upon where 
the event happened. As expected, the raw yield begins to decrease significantly once tracks 
do not cross the entire TPC. Such systematics must be taken into account in the efficiency 
corrections. 

Since the corrected yield we report must be independent of any detector effects, it is 
easier to study the adequacy of the corrections by focusing on a single t] bin and comparing 
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Figure 8.9: Raw yield dependence on ^vertex- The left panel is for the vertex choice [20, 
30] cm (red) and [—30, —20] cm (blue), right panel is for the interval [80, 90] cm (red) and 
[—90, —80] cm (blue). Dashed lines in the bottom panels represent tracks with rj = ±0.5, 
corresponding to the limits given by the vertical dashed lines in the top panels. 



the raw and corrected yield for the bin as a function of the vertex position. This is illustrated 
in Fig. [8.1C which shows the raw and corrected yields for two choices of rj, —0.3 < t] < —0.2 
(left) and 0.4 < rj < 0.5 (right), as a function of ^vertex- The lower panel shows the geometry 
of the relevant tracks as they cross the detector for the different vertex positions. The steps 
in ^vertex were done in 10 cm bin sizes between ±100 cm (20 bins), and 4 wider bins with 
limits [-190, -140], [-140, -100], [100, 140], [140, 190] in order to get enough statistics (see 
the 2:vertex distribution. Fig. |8.2| ). We again see that there are systematic effects in the raw 
yield. We can identify two general trends. For tracks that do not cross the entire TPC, the 
yield drops rapidly. For the remaining tracks, the raw yield is the lowest for tracks closest 
to the central membrane and increases as the tracks get closer to the pad planes on either 
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Figure 8.10: Raw and corrected yield for —0.3 < r\ < —0.2 (left) and 0.4 < t] < 0.5(right) 
as a function of ^vertex- Data points in each of the upper plots correspond to tracks crossing 
the detector region illustrated schematically in the lower plots. 



side of the TPC. For reference, a fit to the raw data assuming a constant yield is done 
(black points) with a resulting poor x^] 6.8 (left) and 8.8 (right). The corrected yield is 
essentially free of ^vertex systematics, except for two cases. In the events closest to the edges 
of the TPC, where the acceptance is varying rapidly and most of the tracks in the event lay 
outside the TPC, the tracks that do not cross the entire TPC show a corrected yield that 
is lower than the rest. This is expected, as tracks with very little information left in the 
detector will be difficult to reconstruct. We also found an additional systematic effect in the 
efficiency correction traced back to the simulation of tracks crossing the central membrane. 
This effect causes a slight over correction of the yield for such tracks, as can be seen by 
the corrected yield for these tracks which is coloured differently in the figure. Therefore, 
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to prevent these systematic effects from appearing in the final pseudorapidity distribution, 
we included only events within ±100 cm of the center of the TPC, excluded tracks going 
through the membrane, and excluded tracks which did not cross the entire TPC. Each data 
point in the final h~ pseudorapidity and p±_ distributions is then obtained as the average of 
each of the ~ 20 independent measurements of the corrected yield obtained for each of the 

■^-vertex bins. 

Since the most important correction for p± spectra and yields is the efficiency, this is 
where we concentrated a large fraction of the studies of systematic effects. The systematic 
uncertainty due to the corrections was estimated in two ways. In order to ensure consistency 
of the results, we studied the variation in the final spectra due to a large variation in the 



track quality cuts. The distributions of the cut variables are given in Fig. |8.1l| from both 
data and simulations; the left panel shows the distance of closest approach distribution 
and the right panel shows the number of fit points distribution. We varied the selection 
based on number of fit points from 10 to 24 and the selection based on distance of closest 
approach from 3 cm to 1 cm. These choices are labelled "cut 1" and "cut 2" in the figure. 
The different corrections for the more stringent set of cuts (cut 2) were then recalculated. 
The final multiplicity distribution obtained in both cases was compared and they were in 
agreement to ~ 1%. To study the sensitivity to the cuts a small variation in the track quality 
cuts was then made, accepting tracks with 23 fit points. A corrected spectrum using the 
correction factors calculated for tracks with 24 fit points was then applied. This yields then 
a corrected distribution that is systematically higher than the measured value, as one is then 
over-correcting the raw data. The yield at high multiplicity obtained via this systematic 
effect was found to be 6.4% above the measured yield. This is the main contribution to the 
systematic uncertainty in the total particle yield for the analysis presented here. 

8.2.3 Backgrounds 

Instrumental backgrounds due to photon conversions and secondary interactions with de- 
tector material were estimated using the detector response simulations mentioned above. 
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Figure 8.11: Comparisons of the distance of closest approach (a) and number of fit points 
(b) distributions for simulated and real data. 



together with events generated by the HIJING model |102i|l03|] . The simulations were cali- 
brated using data in regions where background processes could be directly identified. The 
measured yield also contains contributions from the products of weak decays, primarily Kg, 
that were incorrectly reconstructed as primary tracks that must also be accounted for in the 
background correction. Figure ^.12| (left panel) illustrates the background fraction of the 
raw signal coming from secondary interactions and decays. The shape of the background 
correction is independent of multiplicity. We see that the hadronic background is approxi- 
mately 15% at p± = 100 MeV/c, decreasing with increasing p± to a constant value of ~ 7%. 
The average fraction of hadronic background tracks in the uncorrected sample is 7%. The 
error on the background correction increases with p± because statistics are limited, but 
there is also a systematic uncertainty associated with this correction. Since the background 
depends on the yields of various particles, in particular Kg and A, differences between the 
model and the data will systematically affect the correction. An estimate of the background 
correction at = 2 GeV assuming the ratios A/p ~ 0.8 and p/vr^ ~ 1 would yield a value 
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Figure 8.12: Background correction to the h~ spectrum coming from weak decays and 
secondary interactions in the detector material (a) and electrons (b) as a function o{p±. The 
corrections were obtained via hijing events processed through the STAR offline simulation 
and reconstruction chain. 



of ~ 1.15 instead of 1.07. This is the main source of systematic uncertainty in the hadronic 
background correction. 



The right panel of Figure 8.12 shows the electron background obtained from ^ data 



(filled points) and from HIJING simulations (hollow points). The shape of the electron 
background is exponential, so we can use the shape and normalize to the data in order 
to obtain the final electron correction. Since the slope of the exponential can affect the 
final slope of the p± distribution, we can also do a correction using the HiJiNG points to 
study the variation in the p± spectrum. We find that the choice of the slope in the electron 
correction changes the corrected value of {p±) by 10 MeV/c and is the main contribution to 
the systematic uncertainty of the determination of {p±). 
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8.2.4 Comparison of TPC Halves and Sector {(/)) Dependence 

From the design of the TPC, the tracking volume is essentially divided into two independent 
halves and it is not a priori obvious that one can combine them without introducing system- 
atic effects. Therefore, one of the necessary studies in the evaluation of the reliability of the 
results is to compare the raw yields in the different halves. However, one has to be careful 
in this comparison to isolate the possible differences. Because of the design of the TPC, we 
expect to see differences in the raw yields even within the same half of the TPC simply due 
to an increased drift length. As an example, we can take tracks close to midrapidity, which 
are emitted at an angle close to 90°. As the event vertex is not a fixed quantity in the z 
axis for the different events, tracks at midrapidity coming from events with different vertex 
z positions will have different drift distances, and thus we expect systematic differences in 
the raw yield at fixed rj for varying vertex z positions, as was discussed in Section B.2.2| . 
To isolate differences between the East {z < 0) and West {z > 0) halves of the TPC, we 
therefore have to make sure that we take event and track samples that are related by the 
transformation r]^^^ = — r^West z^grtex = "-^^rertex- addition, as each TPC half is made 
up of 24 independent sectors to cover the full azimuthal range, it is also interesting to do a 
comparison of the raw azimuthal yields to make sure there are also no systematic effects. In 
doing this comparison, one can also separate the raw yields from positively and negatively 
charged tracks, for there might be distortions that affect these tracks differently. The tracks 



selected must also satisfy requirement of having at least 10 fit points ("cut 1" in Fig. 8. 11 (b)). 
This cut is also used for the h~ analysis. 



Figure |8.13| shows the raw azimuthal distribution for negatively charged tracks in the 
two halves of the TPC. The East half track sample, shown as the hollow square data points, 
was obtained from events with —50 < -Zvcrtcx < —10 cm and from tracks with —0.5 < rj < 0. 
Similarly, for the West half the events had a vertex selection of 10 < 2;vertex < 50 cm and 
tracks were selected according to < r/ < 0.5. The left panel shows the yields for the low 
p± region, 0.2 < p± < 1 GeV/c. A similar plot showing the raw azimuthal distribution in 



108 



I. Negative charge, 2 < pi < 1 | 



^ * ^ 
1^ -♦- ^ 



^ ^ * ^ ^ 
1 * 



i * ^ 



Comparison in Yield, East - West 
- □ -50 < z< -10 , -0.5 < 1] < _ 0.2 <p^ < 1 
• 10 < z< 50 , < 11< 0.5, 0.2 < px< 1 



^ , Negative charge, 1 < pi< 2 



50 100 190 200 250 300 350 



+ t + * + + + 



West - East Comparison, Yield 

□ -50 < z< -10 , 0.5 < n < , 1 < Pi< 2 
. 10 c z< 50 , 0. < T) < 0.5 . 1 < pi< 2 

I I , , , . I , 

50 100 150 200 250 300 



(a) 



(b) 



Figure 8.13: Raw (/> distribution of East and West half of the TPC for negatively charged 
particles. Left panel: low p± (0.2 - 1 GeV/c). Right panel: moderate p± [l - 2 GeV/c) 



the range 1 < p± < 2 GeV/c is shown in the right panel. The samples are normalized per 
event to isolate the differences in the raw yields. 

The size of the bins in the azimuthal direction is 15°, or half a sector. We see that there 
is a periodic structure to the distribution in both cases. This comes about because charged 
tracks will curve in the magnetic field. Therefore, depending upon their charge sign and 
entrance point to a sector {i.e. their azimuthal angle), will have a trajectory that is either 
fully contained in a sector boundary or that crosses a sector boundary. For a specific charge 
sign, tracks in one side of the sector will be more easily reconstructed than on the other side, 
giving the structure seen in Figure This effect is more pronounced for low momentum 

{i.e. large curvature) tracks. The high yield bins for the negative charge tracks should be 
the low yield bins for the positive ones. The low momentum positively charged tracks are 
shown in Figure |8.14 Indeed we see the expected change, for this case the first bin is a high 
yield bin (compare to Fig. |8.13| ). 

To focus on changes between East and West halves of the TPC, we make a ratio of the 



previous histograms. Figure |8.15| shows the ratio for low p± negative (a) and positive (b) 
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Figure 8.14: Raw (p distribution of East and West half of the TPC for low transverse 
momentum (0.2 - 1 GeV/c) positively charged particles. 



tracks. The average difference in both cases is on the order of 1%. 
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Figure 8.15: Ratio of raw (p distribution of East and West half of the TPC for low transverse 
momentum (0.2 - 1 GeV/c) negatively (a) and positively (b) charged particles. 

We can also make a comparison of the ratio in both halves of the TPC. Since we 

know that there will be systematic differences as a function of (p due to the curvature effect 
mentioned before, we can simply make a ratio of the raw yields integrated over (p. For the 
East half of the TPC we find h-/h+ = 0.983 ±0.005 and for the West half a similar analysis 
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yields = 0.988 ±0.005 where the errors are statistical only. We conclude that for the 

purposes of this analysis the two halves of the TPC yield sufficiently similar results. 

However, we will still separate the corrections to the raw yields in other variables as 
mentioned in Sec. |8.2.2| . We expect variations in the raw yields due to different track 
geometries and event topology. In summary, we expect that tracks that do not cross the 
entire tracking volume of the TPC to be more difficult to reconstruct than those tracks that 
do. In addition, it is much easier to reconstruct the tracks when there are not many tracks 
in the detector, so we must see a decrease in the tracking efficiency in a high multiplicity 
environment. Furthermore, low momentum tracks coming from the interaction vertex will 
have a difficult time reaching the tracking volume. The radial distance at the center of 
the first TPC sensitive pad row is 60 cm, which along with the magnetic field of 0.25 T 
places an effective low p± acceptance cut-off for primary tracks of ~ 50 MeV/c. Finally, 
as mentioned before even tracks in identical regions of phase-space will have different raw 
yields depending upon where the event vertex was placed because they will sample a different 
detector geometry. Therefore, we have divided all the tracking corrections according to 

• the phase-space cell occupied by the tracks, 

• the multiplicity of the event, and 

• the position of the primary vertex. 

The previous corrections are the most important ones in the present analysis. The following 
corrections to the spectra and yields were also studied. Each of them was found to produce 
changes of less than 1%. 

8.2.5 Momentum resolution 

The momentum resolution is momentum dependent. It is well known that for a p± distri- 
bution, the resolution of the detector will introduce a change in the slope of the observed 
spectrum. This can be seen from the following simple argument. The effect of a finite mo- 
mentum resolution is that a certain number of tracks will be reconstructed with the wrong 
momentum, and will therefore be counted in an incorrect bin. It is possible for a track to 
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be reconstructed with a lower or a higher p±. So for a given p± bin, there will be a loss 
of particles to adjacent bins and a gain of particles from adjacent bins. The magnitudes of 
these fluxes compared to the yield in the given bin are the important quantities. From the 
rapid decrease in cross section with increasing p_L, we expect that the feeding of particles 
from the lower to the higher p± bins will be higher than the flux in the opposite direction, 
and hence the net effect is to flatten the spectrum to some extent. This effect becomes 
important when the momentum resolution is of the same order as the size of the p± bin. It 
is therefore necessary to quantify this effect in any momentum analysis. 

For the results presented here we concentrate on < 2 GeV/c. We find the momentum 
resolution to be better than 4% as illustrated in Fig. |6.6| . For a p± bin width of 100 MeV/c 
we expect the correction due to momentum resolution to be less than 1% in the full range 
of p±. This effect was quantified in two different analyses. The first approach relies on the 
embedding procedure, and consists of dividing the phase space into bins and using the track 
matching between simulated and reconstructed tracks from. We can keep track of the pj_ 
bin in which a particle was created, and then see if it was reconstructed in the same bin. We 
can thus know, for each phase space bin, what percent of tracks are reconstructed correctly, 
what percent are reconstructed in a different bins and correct for the resolution. 

The only significant correction found this way occurs for the very first p± bin. The 
reason is the following. The raw yield in each bin is composed of two parts: the tracks 
which were correctly reconstructed in the same p± bin as the one they were generated in, 
let's call these the healthy tracks; and the tracks that are found by the reconstruction but 
placed in a different bin than the one they were generated in, let's call these the crippled 
tracks. Now let's focus on the region pj_ < 100 MeV/c where the efficiency drops very 
rapidly with decreasing p±. There will be very few healthy tracks in the first p± bin because 
the efficiency is low. There will be, however, a considerable amount of crippled tracks in the 
first bin coming from the feed-down from the next (high efficiency) p± bins. The crippled 
tracks are almost as numerous as the healthy population for this bin. We therefore must 
apply a correction to obtain a realistic estimate of the initial healthy population in that 
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bin. For the region p±_ > 200 MeV/c the efficiency is independent of p± so we need only 
concentrate on the p± dependence of the parent distribution which we want to measure. 

We can also calculate the expected correction based on a knowledge of the momentum 
resolution, the bin size, and a given p± distribution. One way to treat this problem is 
through an iterative procedure, starting with a given input p± distribution, doing a Monte 
Carlo study by smearing the tracks with the measured p± resolution and looking at the shape 
of the resulting distribution. We repeat the process until the output distribution matches 
the one measured in the experiment. The approach we followed relies on a related method. 
Starting from a given input p± distribution we use the measured p± resolution as a function 
of p± to construct a set of Gaussians (one for each p± bin). We also used different functional 
forms to parameterize the shape of the p± resolution, e.g. Lorenzian and double-Gaussian, 
yielding similar results. The area under each Gaussian is the initial yield for each pj_ bin, 
the mean is the center of the p± bin and the a is obtained from the measured 6p± vs. p± 
curve. We can then figure out what is the contribution of each p± bin to any other bin in 
principle. In practice, each bin only contributes mainly to its nearest and next-to-nearest 
neighbours. The observed yield in the i^^ p± bin, N°^^ is then: 



^ p 1 (px-Px(fc))^ 

^obs = V iVfc / -^e~^^ dp^ (8.3) 



We can think of Eq. 8^ as defining a matrix equation 

oo 

iV°b^ = J^iVfcCfci (8.4) 



fc=0 



where the elements of the matrix C^i are defined by the integral of Eq. and are inter- 
preted as the percent contribution of bin k to bin i, i.e. the mean of the Gaussian is the 
center of bin k and the integration limits are given by the upper and lower limits of the i^^ 
p± bin. The width of the Gaussian is given by the detector resolution at p± bin k. Therefore, 
the integral depends only on quantities that are measurable (the resolution) or defined by 
our analysis (the bin limits and bin center). This matrix then embodies all our knowledge 
of the effects of resolution. 
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For a given p±_ distribution, we can then calculate what the observed distribution will 
be and obtain appropriate correction factors for the effect of resolution. We have performed 
this procedure with different input distributions - power-law function, exponential '\n.p± and 
exponential in m±_ with various slope parameters - and find that in all cases the correction 
for any given bin is less than 1% in the range up to 2 GeV/c with our given bin size. 
This is true for both the h~ /h^ analysis (bin size 100 McV/c, range 0.1 - 2 GeV/c) and for 
the 7r~ analysis (bin size 50 MeV/c, range 0.05 - 0.75 GeV/c). For higher p± the resolution 
plays a more significant role, where we expect a correction of 20% or more at 5 GeV/c 
with a 1 GeV/c bin size and the present magnetic field oi B = 0.25 T. 

An additional advantage of this resolution study is that, once we construct the Cki 
resohition matrix, we can also use it to recover the original yields by inverting the 
matrix, and no iterative procedure is needed: 



8.2.6 Track splitting 

For a loose track selection based on the number of points reconstructed in the TPC, it is 
possible to overestimate the track yield in the presence of split tracks. This effect can come 
about for example when the track-finding algorithm fails to recognize two track segments as 
belonging to a single particle. This can happen typically when there are gaps in the track 
pattern, such as tracks crossing sector boundaries, tracks crossing the central membrane of 
the TPC, or tracks crossing a region of the TPC where the read out pad is noisy or dead. 

To study the effect of splitting in this analysis, we used two approaches that relied 
on the detector simulation. An additional study was made using only real reconstructed 
tracks. All three results are in agreement both in the size of the overall effect and on its 
p± dependence. We will discuss the approach that relies on full simulation, since this is the 
most straight-forward. 

In a full Monte Carlo generated event, we have all the information about the input tracks. 



oo 




(8.5) 
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Running the tracks through the GEANT implementation of the STAR detector, we obtain 
among other things the information on the TPC energy deposition left by the tracks. This is 
the input to the microscopic simulation of the TPC response yielding as output a simulated 
raw-data file that is then passed through the STAR reconstruction chain. A comparison 
of the TPC space points that are found by the reconstruction algorithms to the GEANT 



input is then performed. As discussed in Section 7^, the matching as implemented in this 
analysis is based on spatial proximity. A feature of this matching procedure is that it allows 
a many-to-many matching: if 2 hits are very close together for example there will be a 2-to-2 
hit match. Each of the 2 Monte Carlo hits will be matched to the 2 reconstructed hits. 

The match of simulated point to reconstructed point serves as the footing for a simulated 
track to reconstructed track association. Typically, the association is 1-to-l. A split track 
will have a different topology. There will be one original Monte Carlo track, but it will be 
associated to 2 (or more) reconstructed tracks. The association algorithm allows a many- 
to-many type of matching, so one must be careful to really single out a split track from 
other matching topologies. We look for a single Monte Carlo track matched to more than 
one reconstructed track, and those reconstructed tracks are singly matched to the Monte 
Carlo track. Through this procedure, an estimate of the split track population and was 
found to be also below the 1% level. In addition, the analysis done with the more stringent 



requirement on the number of fit points ("cut 2" in Fig. ^. 11 (b)) guaranteed that there are 



is no double counting from split tracks. The results obtained from both analyses were found 
to be consistent. 

8.2.7 Track merging 

When 2 tracks lie very close together, the cluster finder might not be able to resolve the 2 
ionization peaks and produce a single space point when there were originally 2 for each hit 
in the track trajectory. The losses due to this effect were estimated to be less than 1% for 
the most central collisions and negligible for peripheral collisions. 
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8.2.8 Ghost Tracks 

This correction takes into account possible cases where the track finding algorithm might 
associate space-points from different tracks and give us track parameters from a non-existent 
particle. It was not known if this could be of importance in the high multiphcity environment 
of a heavy ion collision at high energy. It was found that the TPC occupancy was low 
enough that this did not present a problem. Essentially no such tracks were found in all 
the simulation and embedding analyses even at higher simulated multiplicities than those 
observed in the data. 

Finally, all analyses were carried out with 3 different software production versions. As 
we understood the systematics of the detector better, we incorporated our knowledge of 
the calibrations, distortions and corrections to software bugs. The first production was 
done in August while data was still being taken. This was followed almost immediately 
by a new production one month later with improved calibrations. The data presented here 
is produced with the calibrations and distortion corrections processed in early 2001. The 
biggest systematic effect observed throughout this process was a 4.5% change in the efficiency 
obtained from the first simulations compared to the following two software versions. This 
was understood as coming from a more realistic parameterization of the allowed space point 
errors used during the track finding algorithm based on the measured residuals obtained in 
the first analysis of the data. 
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Chapter 9 

Results and Discussion: Charged 
Hadrons 



9.1 Multiplicity Distribution 
9.1.1 Results 



Figure |9JJ shows the corrected, normalized multiplicity distribution within \r]\ < 0.5 and 
p± > 100 MeV/c for minimum bias Au + Au collisions. For this plot, a total of 350K 
minimum bias triggers were analyzed, yielding 119,205 events with a reconstructed vertex 
in the region |zvertex| < 95 cm. The distribution for the 5% most central collisions (360 
mbarn), defined using ZDC coincidence, is shown as the shaded area in Fig. The 
definition adopted here was based on the signals shown in Fig. |8.1| , by selecting events with 
ZDC sum < 66 ADC counts and a high CTB threshold placed at 8500 ADC counts to 
remove the ambiguity between central and peripheral collisions. The actual value of the 



ZDC sum threshold is set by requiring the shaded area in Fig. |0| to equal 5% of ctauAu- 
This is the event selection that was used for the h~ p± and r/ distributions for central events 
presented in the following sections. 

The data were normalized assuming a total hadronic inelastic cross section of 7.2 barn 
for Au + Au collisions at ^/s^ = 130 GeV, derived from Glauber model calculations. The 
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Figure 9.1: Normalized multiplicity distribution of h~ with p± > 100 MeV/c, \r]\ < 0.5 in 
Au + Au collisions at y^s^ = 130 GeV. The shaded area is 5% most central collisions, 
selected by ZDC signal (see text). The solid curve is the prediction from the hijing model. 



first bin (below N^^- = 5) is not shown, due to systematic uncertainties in the vertex re- 
construction efficiency and background contamination. Its relative contribution to the total 
cross section was estimated to be 21% by normalizing the hijing multiplicity distribution 



to the measured data in the region 5 < Nfi- < 25 {i.e. the first three bins of Fig. |9.lD . This 
procedure relies on the assumption that very peripheral interactions are well described by 
the superposition of a few nucleon-nucleon collisions in the geometry of a nuclear collision, 
and can therefore be accurately modelled by hijing. A comparison of the actual yield mea- 
sured in the first bin (after the efficiency correction) to the shape predicted by the hijing 
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model indicates a relative difference of 10%. The systematic uncertainty on the vertical scale 
is therefore estimated to be ~ 10% and is dominated by uncertainties in the total hadronic 
cross section and the relative contribution of the first bin. The systematic uncertainty on 
the horizontal scale is 6.4% for the entire range of multiplicity and is depicted by horizontal 
error bars on 3 data points only to maintain clarity in the figure. 

The shape of the h~ multiplicity distribution is dominated over much of the Nj^- range 
by the nucleus-nucleus collision geometry, consistent with findings from lighter systems and 
lower energies. However, the shape of the tail region at large Nf^- is determined by fluc- 
tuations and detector acceptance. These overall features are also observed in the hijing 
calculation, shown as histogram in Fig. |9]l| (although a more detailed comparison indicates 
that, for the model, the shape of the tail shows a different contribution due to fluctuations 
than that found in the data). 

9.1.2 Discussion 

The negatively charged particle multiplicity distribution reveals that the bulk features of the 
dynamics are dominated by the geometry. The shape shows the same features expected from 
a Glauber picture of a superposition of nucleon-nucleon collisions. The model proposed by 



Kharzeev and Nardi [g3| attempts to extract the contribution to the multiplicity distribution 
that scales as either soft or hard processes, or more precisely, as A^part and A'^coU- Such scaling 
is expected to be p±_ dependent, and therefore this study is better suited for a p±_ distribution 
(see Sec. 9.2). 



It is nevertheless useful to apply such models in order to obtain a statistical determi- 
nation of the number of participants for a given event selection. One implementation of 
the Glauber particle production model is illustrated in Fig. This was done by a Monte 
Carlo implementation, where one generates events with random impact parameter, and then 
samples the da/dh geometrical distribution to find the probability to have an event at the 
given impact parameter. When an interaction takes place, we calculate A^part and N^on 



using Eq,4^. We can thus obtain a da/dN^gj-t distribution which we use as a weight in the 
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Figure 9.2: Fit to the multiplicity distribution da/dN^- using a Glauber model (Sec. |4.1.lD 



final multiplicity distribution (similar to the weight in Eq. [1.13| given by the probability of 
no interaction). We can assume that particle production is proportional to the number of 
participants (Sec. {Nfi-{b)) = g(A^part(&)), or use a linear combination of A'part and 

A'^cou a la Kharzeev-Nardi |Q (Eq. 4.1CI| ). As mentioned before, we also choose Gaussian 



fluctuations, Eq. 4.12. The approach works best from mid-central to central collisions, as 



the fluctuations in peripheral collisions are large. Figure |9^ shows the result of applying 
such a model to the STAR multiplicity distribution. The blue curves are the distribu- 
tions obtained by binning the da/dNpa,rt distribution to get successively — 5%, 5 — 10%, 
10 - 20%, 20 - 30%, 30 - 40%, 40 - 50%, and 50 - 60% of ctauAu, which result in the blue 
curves with near-Gaussian tails when transformed to da/dNf^- because of the fluctuations. 
The red curve is the convolution of all such distributions, giving a fit to the multiplicity 



distribution. The size of the fluctuations in Eq. |4.12| is found to be a = 1.45 ± 0.2. This 
parameter is constrained mostly by the shape of the terminus of the distribution. For a 
Poisson distribution, this value would be 1. The presence of correlations makes this value 
exceed unity. The origin of the correlations from the dynamics of the colliding system is 
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a combination of nucleus-nucleus geometry fluctuations, nucleon-nucleon cross section fluc- 



tuations and hadron resonance correlations |151]. In fact, this parameter for HIJING turns 



out to be rather large, o > 4, as is also apparent in Fig.^. There are also fluctuations due 
to finite detector acceptance. Since the TPC has a large acceptance, the contribution from 
acceptance fluctuations is small. The shape of the terminus of the distribution measured 
with a smaller acceptance apparatus should be reveal a less steep shape of the terminus of 
the distribution {i.e. a larger a), e.g. the PHENIX Pad Chamber measurement |152|. 



It is argued |153| that the centrality dependence of the rapidity density per participant 
pair can distinguish between models based on particle production from gluon saturation 
(EKRT [|154|| ) from those based on fixed scale pQCD such as hijing. In the EKRT model, 
particle production is computed assuming that the pQCD growth of low p± gluons is only 
integrated up to a certain saturation scale. The saturation requirement has as a conse- 
quence that the multiplicity per participant decreases with centrality. For HiJiNG, particle 
production arises from two main contributions: (a) a component arising from low-p± hadron 
production, modelled as beam jet string fragmentation, which is essentially proportional to 
-^part, and (b) a contribution arising from mini-jet production, which is directly proportional 
to A'cou and to the averaged inclusive jet cross section per nucleon-nucleon collision. This 



is essentially the same functional dependence as Eq. |4.10| . The centrality dependence of the 
mid-rapidity density for HIJING is then 



= (A^part>rasoft + fiNcoll) ^jyjV (^'^^ 



where risoft ~ 1-3 and / « 1.2. A plot of dN / dr] / {Npi^j-t / 2) vs. A^'part should then naively 
grow as ~ Apart^''^- Using the Monte Carlo model just discussed, and also the model from 
Kharzeev and Nardi p3| applied to the STAR data, we obtain the centrality dependence 



of dN/drj shown in Fig. |9.3| . The results from PHENIX [|152| and from the Kharzeev and 
Nardi analysis of the PHOBOS results |^] are also shown. The growth with Apart clearly 



disfavours the saturation picture from Ref. [ISJ]. As pointed out in Ref. ||93|, the original 
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Figure 9.3: dN/drj per participant pair vs. A'part- 

saturation ideas pertain to the behaviour of partons in the initial wave function of the nu- 
cleus; and when taking this approach the predictions surprisingly turn out to be quite similar 
to the eikonal Glauber approach discussed in Sec. There are also deviations from the 

simple A'part^''^ dependence which can arise from the dilute edges of the Wood-Saxon nuclear 



distribution (Fig. 4.1) and from other medium effects such as nuclear shadowing of the initial 



parton distributions ||155|| or jet quenching |7^, 77|. The similarity between the conventional 
approach and the high-density QCD approach makes it difficult to distinguish which un- 
derlying picture is correct. These effects are therefore better studied in the differential p± 
distributions. 



9.2 p± Distribution 
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9.2.1 Results 

Since the yield of charged hadrons at midrapidity depends on the extrapolation of the p± 
distribution outside the measured range, we discuss the p± spectrum first. The shape of the 
p± distribution is well suited for more detailed scrutiny. Figure |9!J shows the p± distribution 
for positively and negatively charged hadrons (/i"*" and h~). Statistical errors are smaller 
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Figure 9.4: The p± distributions for h~ (open circles) and h'^ (open "plus" sign) in central 
collisions. The inset shows the ratio h" /h^ vs. pj_. The data are fit to a power law in m± 
to obtain the extrapolated yield dN/dr] (see text). 



than the symbols. The correlated systematic uncertainty is estimated to be below 6%. 
The pj_ distributions in this range of centre-of-mass energies are customarily fit by a QCD 
inspired power-law function of the form d'^ Nf^- / dp^dr] = A{1 +p±/po)~" where A, n, and 
Po are free parameters. This power-law fit is useful in comparison to other reference data 



such as p± spectra from pp collisions [OGj]. For the extrapolation to p± = 0, it is better 
to use a modified version of the power law, and use m± instead of p± in order to better 



reproduce the data at low p±. This is the fit shown in Fig. 9.4, using a mass slightly larger 
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than the vr" mass, since the /i^ spectrum is made not only of pions but includes and p 
as well. The "mean mass" of the h~ distribution must therefore be (m) > m-,^. We studied 
the effect that such variations in the shape of the extrapolating function have on the yield 
by comparing the m_\_ and p±_ power-law fits and also a simple exponential fit at low p_\_. 
The yield in the region — 100 MeV/c was found to be 6.8% using the power law, 
the other fits varied from this by about ±1.5% which we take as an additional systematic 
uncertainty for the extrapolated yield. The choice of mean mass does not influence the 
extrapolated yield obtained from the power law, it was found to vary less than 0.5% 
with mean masses in the range 140 — 200 MeV/c^. The \^/d.o.f. was minimized with a 
value of (m) ~ 187 MeV/c^. In all cases, the yield in the region p_\_> 2 GeV/c was < 1%. 
We therefore use an extrapolation factor of 7 it 1.5% to obtain dN/drj. 

9.2.2 Discussion 

We focus now on the comparison of the measured p_\_ spectrum in Au + Au at RHIC to 



other energies and collision systems. We will focus more on pp collisions. Figure |9.5| , upper 
panel, shows the transverse momentum distribution of negatively charged hadrons for central 
Au + Au collisions at mid-rapidity {\ri\ < 0.1) within 0.1 < p± < 2 GeV/c. The upper panel 



of Fig. |9^ also shows the p^-distributions of negatively charged hadrons for central Pb + Pb 



collisions at ^/s^= 17 GeV from NA49 |10£] and for minimum-bias pp collisions at \/s = 



200 GeV from UAl @, fitted with the same function. The NA49 distribution, which 
was reported in units of pion rapidity, was transformed to units of pseudorapidity. The 
UAl invariant cross section Ed^a/d^p reported in Ref. was scaled by 27r/(Tinei, where 
o"inei = 42 mb for = 200 GeV |^. The power law fits all three datasets well. The mean 
p± can be derived from the fit parameters as {p_i_) = 2pQ/{n — 3). The fit to the STAR data 
gives po = 3.0 ± 0.3 GeV/c, n = 14.8 ± 1.2, and {pi_) = 0.508 ± 0.012 GeV/c. The strong 
correlation of fit parameters pq and n must be taken into account when calculating the error 
on {pa.)- The {p_[_) from STAR is larger than that from both central collisions of heavy 
nuclei at much lower energy ((p±)na49 ~ 0.429 GeV/c) and nucleon-nucleon collisions at a 
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Figure 9.5: Upper panel: h~ p^-spectra for the 5% most central Au + Au collisions at mid- 
rapidity (|?7| < 0.1). Central Pb + Pb data from NA49 at lower energy and pp data from 
UAl at a similar energy are also shown. The curves are power-law fits to the data. Lower 
panel: ratio of STAR and scaled UAl p^-distributions (see text). 



comparable energy ((p_l)uai — 0.392 GeV/c). 



Figure |9^, lower panel, shows the ratio of the STAR and UAl p^-distributions. Since 
the UAl distribution is measured at ^/s= 200 GeV, da/dp± is scaled by two factors for 
quantitative comparison to the STAR data at 130 GeV: (i) i?(130/200), the p^-dependent 
ratio of invariant cross sections for charged particle production in pp collisions at ^/s = 
130 and 200 GeV, and (ii) Taa= 26 it 2 mb~^, the nuclear overlap integral ||2^ for the 
5% most central Au -|- Au collisions. R varies from 0.92 at p± = 0.2 GeV/c to 0.70 at 
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p± = 2.0 GeV/c, and was derived using scaling laws for {p±) and dNch/dq as a function of 
[|^,|9^ together with the extrapolation to 130 GeVof power-law parameterizations at 
V^= 200-900 GeV (9^. The shaded boxes show the total error of the ratio, which is the 
linear sum of the errors of the measured data, depicted by the error bars, and the systematic 
uncertainty due to uncertainties in the scaling with Taa and R. 

There are two simple predictions for the scaled ratio. In lower energy hadronic and 
nuclear collisions, the total pion yield due to soft (low p±) processes scales as the number 
of participants (or "wounded" nucleons) in the collision (see e.g. Ref. |31, 105(| ). The scaled 



ratio in this case is 0.164, assuming 172 participant pairs and a mean number of binary 
collisions of 1050 (= dinei Taa, o"inei = 40.35 mb) for the 5% most central Au + Au events 
(Sec. |4.1.1| ). In contrast, if hadron production is due to hard (high p±) processes and 
there are no nuclear- specific effects, the hadron yield will scale as the number of binary 
nucleon-nucleon interactions in the nuclear collision and the value of the ratio is unity. 
There are important nuclear effects which will alter the scaling as a function of p± from 



these simple predictions, including initial state multiple scattering |121|, shadowing |155| 



jet quenching ||76| , |77|| , and radial flow ||156|| . The scaled ratio exhibits a strong dependence, 
rising monotonically with increasing p± from Wounded Nucleon scaling at low pj_ but not 
reaching Binary Collision scaling at the highest p± reported. This is consistent with the 
presence of radial flow, as well as the onset of hard scattering contributions and initial state 
multiple scattering with rising pj_. Prom pA collisions at lower energy, the data tend to 
show an increase beyond the scaling with binary collisions. This "Cronin" effect is thought 
to arise from initial state multiple scattering in the nuclear medium, which broadens the p± 



spectrum [121|. At SPS, the p± spectra indicate a similar increase beyond binary collision 
scaling for the region p± ^ 2 GeV/c. Although a comparison of the p± dependence of 
the data at different ^/s^ is better done in x± = 2p^l\fs rather than in p_L, there are 



indications from preliminary STAR jOlJ and PHENIX [ll57|| hadron spectra at high pj^that 
the scaling with A'coU is not reached even for 2 - 6 GeV/c. This effect has sparked active 
discussions in the community vis-a-vis the possibility of being related to parton energy loss 
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in a QGP, i.e. jet quenching [p.58| . 



9.3 T] Distribution 
9.3.1 Results 

By integrating the p± spectra for different rj bins, we obtain the dN/dr] vs. rj distribution. 



Figure 9^ shows the normahzed pseudorapidity distribution of h~ and for the 5% most 
central colhsions within \rj\ < 1.0. The black points are the measured yield for 0.1 < p± < 
2 GeV/c and the hollow points are the extrapolation for all p±. The latter was obtained by 
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Figure 9.6: h~ (red circles) and h'^ (blue squared) pseudorapidity distribution from 5% most 
central Au + Au collisions. Data are integrated over pj_ > 100 MeV/c (filled symbols), and 
extrapolated to all p±_ (open symbols). 



fitting a power-law function in the range 0.1 < p± < 2 GeV/c and extrapolating to = 
in order to estimate the content of the first p± bin, as discussed in Section |9]^. The error 
bars indicate the uncorrelated systematic uncertainties. The statistical errors are negligible. 
The correlated systematic uncertainty applied to the overall normalization is estimated to 
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be below 6% for pi_ > 100 MeV/c and 7% for all p_L. 

The h~ density at midrapidity for 0.1 < p_L < 2 GeV/c is dNf^- / dr]\r^=Q = 261±l(stat)ib 
16(syst). Extrapolation to p_L = yields dNf^- / d7]\,^=Q = 280±l(stat)±20(syst). Assuming 
an average of 172 participant pairs per central Au + Au collision, this corresponds to 1.63 ± 
0.12 /i~ per participant nucleon pair per unit pseudorapidity, a 38% increase over the yield 
in pp collisions extrapolated to the same energy (neglecting isospin correction factors 
of order 1-3%) and a 52% increase over Pb + Pb collisions at y/s^= 17 GeV pi. The 



corresponding h'^ density for 0.1 < p_L < 2 GeV/c is dNii+ / dr]\rj=o = 268±l(stat)±16(syst), 
which after the extrapolation to all p± yields dNf^+ / dr]\^=Q = 287 ± l(stat) ± 20(syst). The 
systematic uncertainty is the same as for the h~ since the analysis was done in the same 
way. It is a correlated systematic for both charge sign yields, so it does not mean that the 
yields can be equal within errors. In other words, the ratio h'^ /h~ is 1.025 it 0.004 (where 
most of the systematic uncertainties cancel and we are left with a small statistical error). 
The PHOBOS collaboration has reported a total charged multiplicity density for the 6% 



most central Au + Au collisions of dNch/dri\\^\^i = 555±12(stat)ib35(syst) [159|. Analyzing 
positive charged particles within the framework described above, STAR measures the total 
charged particle density (iA'^c/i/c^^||7;|<o.i = ^67 ± l(stat) ± 34(syst) for the 5% most central 
Au + Au collisions. The PHOBOS centrality selection is based on charged particles within 
3 < |?7| < 4.5, while STAR utilizes spectator neutrons in the ZDCs. A precise equality of 
the two measurements is therefore not expected, due to the difference in centrality selection. 

9.3.2 Discussion 

The rj distribution is almost constant within \r]\ < 1, exhibiting a small rise at larger rj. 
This shape is expected from a boost invariant source (i.e., constant in rapidity), taking into 
account the transformation from y to ry. As a reference, the data are fit by the shape of the 



Jacobian dy/drj given in Eq. \i.2l\ , where the fit parameters are (i) the normalization constant 



(which should be equal to dN^h/dy for a boost invariant source) and {ii) the dimensionless 
ratio p±/m, or rather {p±_) / (m), since the 7] distribution is already integrated over p± and the 
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mass should be a weighted average over all 3 species. Prom the simple fit of the distribution 
to Eq. |4.21j , we obtain the parameters dNj^- /dy = 295 ± 18 and dNi^+ /dy = 304 it 18. Using 
the value of (j)±) from the p±^ distribution, we can convert the other fit parameter into 
an average mass. As discussed previously, the value should be close to the vr" mass since 
pions dominate, but we expect it to be somewhat larger due to the contribution from kaons 
and anti-protons. We obtain (m) = 0.171 GeV/c^, which is consistent with this simple 
picture (and close to the value that minimizes the for the power-law fit in m_\_ to the 
h~ p_L distribution, (m) = 0.187 GeV/c^). This signifies that the system is very close to 
a boost-invariant source. Measurement of the rapidity distribution of identified particles is 
needed to test boost invariance at mid-rapidity, and will be addressed in the n~ analysis. In 
particular, the rapidity distribution and the rapidity dependence of the slopes will be useful 
to address this question, as we find a slight but systematic difference in the tt" slopes as a 
function of rapidity. 

The approach to boost invariance is nevertheless quite striking as compared to lower 



energy Pb + Pb collisions. Figure |9/7| shows the rapidity distribution of charged hadrons 
at SPS energies. The collision systems and beam energies are from central Pb + Pb @ 
Pbeam = 158 AGeV/c (hollow squares), central S+S @ Pbeam = 200 AGeV/c (asterisks), 
and minimum bias nucleon-nucleon collisions are also shown (hollow rhombi). The data 
are shown as dN/dy vs. y where y is calculated assuming the pion mass. The shape of the 
distribution is peaked around mid-rapidity in contrast to the relatively flat region we find 
at ^/s^ = 130 GeV. The forward and backward region should show deviations from this 
behaviour, and it is interesting to study the size of the flat mid-rapidity region. The forward 
TPC's in STAR will be in a position to answer this question. From other experiments. 



results from the PHOBOS collaboration []160|| suggest that the plateau region extends up 
\rj\ < 1.5. The ratio of positive to negative particles /h~ = 1.025±0.004 already indicate 



that the mid-rapidity region is not net-baryon free, as has also been reported by STAR |]161 | 



and BRAHMS [162| for the p/p ratio. 
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Figure 9.7: r/ distribution of negative hadrons at SPS (shown as in the center-of-mass) 



Data are from Ref. 1 105 1 



Lower Limit on Energy Density 

With the h~ multiphcity and {p±) we can return to the question of estimating the energy 
density. As discussed in Section |2.2.1 , we can only obtain a lower limit with these simple 
observables. We already stated that there is no reliable guidance as to the value of the initial 



formation time of the system, r in Eq. p.l4 We assume a value of Ifm/c, similar to the one 
adopted at SPS energies, although the system might equilibrate faster at RHIC energies. 
This is not the only assumption that goes into such an estimate. The initial energy density 
is also modified, possibly by a large amount, by the longitudinal expansion. Transversal 
expansion does not change dE±/dy, but is not known how much longitudinal work is done 
by the system during the expansion. In this sense, it could be misleading to use the initial 
formation time r but a final dE±/dy. Since we measure multiplicity instead of transverse 
energy, this might turn out to be a better estimate since dN/dy is approximately conserved 
during the expansion. However, it is more difficult to estimate the time dependence of {p±)- 
The temperature of the system decreases with time, but the collective flow increases. In 
the best case, both effects would cancel, but it is more likely that the build up of flow does 
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not fully compensate the decrease in temperature. This also results in an estimate of the 
energy density that is a lower limit with respect to {p±). It is also a pre-requisite that the 
behaviour of the system in the initial stages can be modelled via hydrodynamics. Given a set 
of energy-momentum conservation and particle number conservation equations along with an 
equation of state, one can then calculate the energy density, pressure and entropy density 
of the system. If the entropy density is approximately conserved during the expansion, 
and under the assumption that the entropy density is approximately proportional to the 



particle multiplicity in the final state |5S, 108 1, then one can estimate the energy density as 



in Eq. 2.14 



The estimate under these assumptions proceeds as follows. For central collisions, we have 
dNch/dr] = 567 and {p±)= 0.511 GeV/c. We could also use {m±) instead of (p_L)because 
one should not neglect the masses. In the best of all worlds one actually would obtain the 
{m_i_) for each particle and do a sum weighing by the fraction each particle contributes to 
the total h~ distribution. Since we already mentioned the larger uncertainties that go into 
the estimate, such refinements can be ignored. In addition we do not know {p±) for kaons 
and protons yet, nor do we know the ratios /h^ , and p/h~to a degree that we can do 
considerably better than the simple assumption. Since we measure only the energy density 
carried by the charged hadrons, we also use Bjorken's guess e ~ ec/i x 1.5. The resulting 
equation is: 

,>:^5x (9.2) 

The parameters we enter into the formula are the following. We use the same value of r = 1 
fm to compare to SPS (although it probably takes less time to equilibrate the system at 
RHIC energies). This also allows the estimate to be scaled to a different formation time 
once it is better determined. 

The radius (r = tq x A^^^) we use can be simply taken as the one for 6 = 0. We must 
then extrapolate the multiplicity, since it was measured for finite impact parameter. In the 
geometrical model, 5% corresponds to roughly 6 = — 3 fm. For head-on collisions 6 = the 
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/i^ multiplicity increases with respect to the 5% most central by roughly 14% {Nf^- ~ 320 
for & = from Fig. |9.1|) . In addition, we are measuring dN/drj not dN/dy. We can use the 
estimate of dN/dy from the fit to the tj distributions using the Jacobian dy/drj. 

We then have the parameters t = 1 im, A = 197, dNch/dy = (295 + 304) x 1.14 = 682, 
{m±) = \J {pA_f + = \/0.51l2 + 0.1712 = 0.538 which we can insert into Eq. |9l 

obtaining 

The ~ 10% systematic uncertainty is based on variations of (p^) and the , multiplicity 
estimates within their systematic uncertainties. 

This is actually somewhat small when one consider that NA49 reports 3.2 GeV/fm^ for 
Pb+Pb via calorimetry {dE±/dr]). Using tracking they get however 2.6 GeV/fm'^ ||163|| for 
charged particles. The rather small difference between 3.9 and the NA49 result 3.2 also 
points to the fact that r should indeed be smaller for RHIC energies. Other scenarios, such 
as a denser core and a less denser shell, cannot be probed with this simple method since 
this is just the average density of the source volume. It is probably best to obtain estimates 
from actual hydrodynamical simulations that reproduce measured data {e.g. particle yields 
and spectra, elliptic and radial flow) and obtain from these calculation ranges of values for 
the relevant parameters, i.e. excluded regions of e. However, for a mid-rapidity region with 
almost zero baryochemical potential. Figure |2]^ already indicates that at T = 1.5 Tc, the 
energy density in the deconfined phase (for the case of 2+1 flavour QCD) is e « 3.5 GeV/fm^, 
so our simple estimate would indicate a favorable environment to reach the deconfined phase. 

9.4 Centrality dependence 
9.4.1 Centrality selection 

It is interesting to study how the distributions are affected by the centrality of the collision. 
As we mentioned in the discussion of the Glauber model, it is not possible to experimen- 
tally determine the impact parameter of the collision. In practice, any study on centrality 
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dependence relies on selecting ensembles of events based on an observable correlated with 
the impact parameter. Typically these can be transverse energy production, or the charged 
particle multiplicity of the event. That is why the study of charged hadrons is of primary 
concern in further studies of heavy ion collisions, such as strangeness production. The focus 
of the current work however, is the study of charged hadrons, so to study the centrality 
dependence of charged hadron production is somewhat ill-defined. The approach we fol- 
lowed for the centrality dependence relies on a selection of events based on the ZDC and 
CTB trigger signals. Figure The most central collisions are studied by selecting events 
below a ZDC threshold. Figure |9.8| (a). This leads to an ambiguity, since both central and 




5 

olO 



10 



10 



10 



10 



STAR, > 1 00 MeV/c, I T] I < 0.S 

■ Top 5% 

■ 10-5% 
20 - 10% 

■ 30 - 20% 
40-30% 
50 - 40% 




5000 10000 15000 20000 25000 
Central Trigger Barrel Signal (ADC) 



50 100 150 200 250 300 350 
Number Of Negative Hadrons 



(a) 



(b) 



Figure 9.8: Centrality classes for h analysis. Left panel: cuts in ZDC vs. CTB. Right 
panel: Resulting classes in the h~ multiplicity distribution. 

peripheral collisions can have a low ZDC signal, so we set a CTB cut only to break this 
ambiguity (the cut was set at 8500 ADC counts, not shown in the figure in the interest of 
clarity). This approach works up to a point, since the ZDC signal turns around and becomes 
correlated with multiplicity instead of anti-correlated, and we must therefore abandon the 



the simple ZDC threshold cuts in place of cuts in both ZDC and CTB. Figure 9.8, left panel 
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shows the cuts in ZDC vs. CTB signals; the right panel shows the resulting event selection 
in the h~ multiplicity distribution for reference. 

9.4.2 p± and r] vs. Centrality 

With the previous centrality selections, we obtain a (PN/dp±drj distribution for each event 
sample. We again apply the tracking corrections (Sec. |8.2| ) for each event and accumulate the 
corrected distributions for each centrality class. The r] distribution for the different centrality 
bins is shown in Fig. 9^. The shape of the distributions at mid-rapidity (|r/| < 0.5) is very 
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Figure 9.9: The h ij distribution for different centralities. 



similar for all the centrality selections. The p± distribution for the different centrality classes 



is shown in Fig. 9.10. The dN/drj and {p^) values are collected in Table Blj while the p± 



distributions for the various centralities are collected in Table 9.2 



It is instructive to do a comparison between peripheral and central p±_ distributions. It 
is easier to study differences in the distributions by plotting the ratios of the p^ spectra, 



since slight changes in curvature are difficult to see in a log scale such as in Fig. |9.10| . As 
a reference, we take the shape of the measured p± distribution for the most peripheral bin 
(70 — 80% of o"AuAu)- We then take the ratio of each of the distributions in Fig. p. 10 to 



134 



^(%) 




) [GeV/c] 


dN/drj (Meas. 


) dN/drj (Bxtrap.) 


(iVpart) 


(iVcou) 


{b) fm 


- 5% 




0.511 


266.5 


289.8 


345 ± 7 


1050 ± 19 


2.3 ± 0.2 


5-10% 




0.511 


222.5 


242.2 


289 ±9 


826 ± 25 


4.3 ±0.3 


10-20% 




0.509 


169.7 


184.9 


221 ±4 


566 ± 10 


5.8 ±0.3 


20-30% 




0.506 


123.5 


134.7 


152 ±9 


336 ± 25 


7.5 ±0.3 


30-40% 




0.495 


76.1 


83.1 


102 ±4 


209 ± 15 


8.7 ±0.3 


40-50% 




0.482 


44.7 


48.9 


63 ±4 


105 ±4 


10.1 ±0.4 


50-60% 




0.467 


24.4 


26.8 


35 ±3 


36 ± 5 


11.2 ±0.6 


60-70% 




0.453 


12.7 


13.9 


20 ±2 


22 ±4 


12.0 ±0.7 


70-80% 




0.448 


7.1 


7.8 


9±4 


9±4 


12.6 ± 1.1 


Table 9.1: 




) and dNf^- 


-/drj (in the 


region \r]\ < 0.5) for various 


centralities selected as 



fractions of the cross section ctauAu (Fig. The quantities statistically related to the 

multiplicity using a Glauber model calculation are also shown for reference. 



the most peripheral. The simplest exercise is to take the ratio without any further scaling. 
This is shown in Fig. 9.11| (a). If we take the most peripheral distribution as the one which 
most closely should follow the spectrum from pp or pp collisions, we see how the shape 
of the p± distributions change with centrality. Similarly to the pp comparison, we see an 
increase in the h~ yield with increasing p± for a particular centrality selection compared to 
the peripheral distribution. We also see that the excess at fixed increases with centrality. 
This is more easily seen in a double ratio, where we normalize each p± distribution by the 
corresponding dN/drj, i.e. by the integral with respect to p±. This double-ratio is shown in 



Fig. p.ll| . The double-ratio curves then meet at low-p± independent of centrality, and the 
rise with increasing p± is evident. We find that for the most central collisions, the h~ yield 
at 2 GeV/c is almost a factor of 2 larger than expected from a superposition of peripheral 
collisions that would integrate to the same yield. This change in behaviour is consistent with 
the presence of collective effects such as increased radial flow for central collisions. However, 
as mentioned before, the p± spectrum for even the most central collisions still does not 
reach the expected scaling with A'cou at 2 GeV/c. The measurement of the p± distributions 
at higher p± will help be an interesting extension to these studies, as they will shed more 
light on the modification of the pj_ spectra, e.g. whether the "Cronin" effect observed in pA 
collisions and lower energy AA collisions should be seen at these moderate momenta, and if 
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Figure 9.10: The h~ pj_ distribution for different centralities. The data are scaled down by 
successive powers of 10 for display purpose. 

so, whether the jet quenching hypothesis is able to reproduce the centrality dependence of 
the p± ratios such as Fig. ^]5| and Fig. |9.11 . 



Figure |9.12 shows the multiplicity dependence of {p±) for h~ . We use the fitted parame- 
ters to the power law in Fig. |9.10 to obtain {p±), the corresponding multiplicity is obtained 
by averaging the measured yield in the region |r/| < 0.5, 0.1 < p± < 2GeV/c (Fig. p.Op and 
extrapolating to all p± using the power law. This has the effect that the average multiplicity 
measured in the wider slice is higher (267, compared to 261). The data are collected in 



Table 9J., which shows the multiplicity and {p±_) obtained from the spectra for the 
different fractions of cross section. The error bars representing the systematic uncertainty 
are 2% for {p±), 6.4% for measured and 7.1% for extrapolated dN/dr]. By using the power- 
law fits to the pj_ distributions in this range, we obtain the missing fraction in order to 
obtain a total yield (~ 7% for < p± < 0.1 and 1% for p± > 2 GeV/c). The centrality 
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(a) (b) 

Figure 9.11: Ratio of the p_\_ distributions for different centralities to the p_\_ spectrum of 
the most peripheral bin (70-80%). Left paneh simple ratio of the data. Right panel: Ratio 
of data points to the power law fit to the most peripheral bin including a normalization of 
each distribution to the corresponding mid-rapidity density dN/drj of each sample. 



selection is based on cuts in ZDC and CTB trigger signals as shown in Fig.|9.8. 

We see an increase in {p±) in the first peripheral bins and then a saturation. The value 
of increases by 18% compared to SPS energies, when comparing to central collisions. 
To put into perspective the stiffness of the observed Au + Au p_\_ spectra, the measured (p^) 
at = 900 GeV reported in ||9l| is = 0.447 ± 0.003 GeV/c. At the highest Tevatron 



energies, CDF reports |164| {p^) = 0.495 ± 0.014. The saturation of {p±) contrasts the 
structure proposed in Ref. ||65| , where the qualitative dependence for a first order transition 
(under the assumption that the central multiplicity is approximately proportional to the 
sV, where V is the volume of the system and s = e + p is the entropy density in Landau 
hydrodynamics, and that {p±) reflects the temperature of the system) should be an initial 
correlation between dN/drj and followed by a saturation, and then a subsequent rise. 

We find no such structure. This behaviour helps to strengthen the growing evidence against 
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Figure 9.12: The multiplicity dependence of (px) for h . The values obtained from the 



UAl IQ^I and NA49 [[105| p±_ distributions in Fig. |9^ are also shown. 



a strong first order transition. In addition, it is also consistent with current expectations 
of a second order phase transition region at higher baryo-chemical potentials with a critical 
point followed by a region of a smooth cross over between the QGP and hadron gas phase 
below some as yet undetermined value for ^b- Initial estimates from particle ratios at RHIC 
yield values of fiB — 45 MeV, which would then most probably place RHIC collisions at this 
energy in the cross-over region. This will of course need to be studied with more sophisticated 
analyses, in particular, an interesting future direction could be the measurement of the {p±) 
vs. dN/di], along with {p±) fluctuations event-by-event, studied for different centre-of-mass 
energies or collision systems to search for the possible critical point. Theoretical guidance 
in this direction would also be fruitful in this search. 
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p± 




Fraction of crhadronic 




GeV/c) 


0-5% 


5-10% 


10-20% 


20-30% 


30-40% 


0.15 


2226.0 ± 94.8 


2195.9 ± 84.0 


1689.4 ± 53.6 


1244.6 ± 42.5 


787.3 ± 21.6 


0.25 


1729.1 ± 36.9 


1514.3 ± 52.5 


1162.0 ± 33.3 


855.7 ± 26.4 


537.8 ± 13.3 


0.35 


1117.1 ± 26.1 


965.8 ± 34.5 


734.8 ± 21.8 


532.6 ± 17.0 


331.2 ± 8.6 


0.45 


732.2 ± 13.7 


588.7 ± 22.6 


450.1 ± 14.3 


326.3 ± 11.2 


201.4 ± 5.6 


0.55 


478.8 ± 7.3 


381.8 ± 15.9 


291.2 ± 10.1 


211.2 ± 7.9 


130.1 ± 4.0 


0.65 


304.3 ± 6.0 


251.2 ± 11.5 


190.0 ± 7.2 


137.4 ± 5.7 


83.8 ± 2.8 


0.75 


204.7 ± 4.0 


169.6 ± 8.5 


128.6 ± 5.3 


92.7 ± 4.2 


56.4 ± 2.1 


0.85 


139.5 ± 4.6 






UO.U -L O.J. 


38.0 ± 1.5 


0.95 


97 6 ± 2 


80.7 ± 4.8 


60.4 ± 3.0 


43.8 ± 2.4 


26.5 ± 1.2 


1.05 


68.4 ± 1.9 


56.3 ± 3.7 


43.1 ± 2.3 


31.2 ± 1.8 


18.47 ± 0.89 


1.15 


47.0 ± 1.5 


39.8 ± 2.8 


30.3 ± 1.8 


21.7 ± 1.4 


12.71 ± 0.67 


1.25 


34.2 ± 1.2 


28.0 ± 2.2 


21.4 ± 1.4 


15.33 ± 1.07 


9.03 ± 0.51 


1.35 


23.8 ± 1.0 


20.1 ± 1.7 


15.24 ± 1.06 


10.95 ± 0.83 


6.33 ± 0.39 


1.45 


16.6 ± 0.9 


14.20 ± 1.29 


10.99 ± 0.82 


7.69 ± 0.63 


4.50 ± 0.30 


1.55 


11.72 ± 0.74 


10.46 ± 1.02 


7.91 ± 0.64 


5.70 ± 0.50 


3.25 ± 0.23 


1.65 


8.69 ± 0.46 


7.63 ± 0.79 


5.76 ± 0.50 


4.02 ± 0.38 


2.32 ± 0.18 


1.75 


6.20 ± 0.43 


5.61 ± 0.62 


4.18 ± 0.39 


2.96 ± 0.30 


1.68 ± 0.14 


1.85 


4.44 ± 0.39 


3.93 ± 0.47 


2.99 ± 0.30 


2.14 ± 0.24 


1.20 ± 0.11 


1.95 


3.38 ± 0.32 


2.88 ± 0.37 


2.20 ± 0.23 


1.57 ± 0.18 


0.851 ± 0.081 


P± 




Fraction of crhadronic 






GeV/c) 


40-50% 


50-60% 


60-70% 


70-80% 




0.15 


472.0 ± 14.9 


268.5 ± 10.7 


143.5 ± 6.4 


80.6 ± 4.8 




0.25 


325.3 ± 9.3 


183.1 ± 6.6 


97.1 ± 3.9 


54.7 ± 2.9 




0.35 


198.8 ± 5.9 


110.1 ± 4.1 


57.8 ± 2.4 


32.6 ± 1.8 




0.45 


119.1 ± 3.8 


65.3 ± 2.7 


34.5 ± 1.6 


19.3 ± 1.2 




0.55 


76.0 ± 2.7 


41.2 ± 1.8 


21.2 ± 1.0 


11.91 ± 0.78 




0.65 


48.8 ± 1.9 


26.3 ± 1.3 


13.37 ± 0.72 


7.26 ± 0.53 




0.75 


32.3 ± 1.4 


16.90 ± 0.92 


8.68 ± 0.50 


4.67 ± 0.37 




0.85 


21.5 ± 1.0 


11.42 ± 0.68 


5.56 ± 0.36 


2.97 ± 0.26 




0.95 


15.12 ± 0.76 


7.78 ± 0.51 


3.84 ± 0.26 


2.08 ± 0.20 




1.05 


10.03 ± 0.56 


5.18 ± 0.37 


2.52 ± 0.19 


1.35 ± 0.137 




1.15 


6.89 ± 0.41 


3.50 ± 0.27 


1.73 ± 0.14 


0.910 ± 0.098 




1.25 


4.97 ± 0.32 


2.49 ± 0.21 


1.14 ± 0.10 


0.629 ± 0.072 




1.35 


3.50 ± 0.24 


1.71 ± 0.15 


0.828 ± 0.072 


0.449 ± 0.053 




1.45 


2.47 ± 0.18 


1.23 ± 0.12 


0.592 ± 0.053 


0.317 ± 0.039 




1.55 


1.76 ± 0.14 


0.866 ± 0.087 


0.403 ± 0.037 


0.210 ± 0.027 




1.65 


1.24 ± 0.11 


0.630 ± 0.066 


0.279 ± 0.026 


0.154 ± 0.020 




1.75 


0.940 ± 0.082 


0.438 ± 0.049 


0.205 ± 0.019 


0.120 ± 0.014 




1.85 


0.652 ± 0.061 


0.330 ± 0.038 


0.151 ± 0.013 


0.089 ± 0.010 




1.95 


0.474 ± 0.046 


0.226 ± 0.027 


0.115 ± 0.009 


0.061 ± 0.006 





Table 9.2: h yield (PN/{p± dp± drj) at r] = for different centralities. 
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Chapter 10 

Results and Discussion: Identified 
Pions 



The next step in the analysis of charged particles is to use the particle identification ca- 
pabilities of the TPC. This allows us to select a specific particle for further study. The 
first natural candidate is the n meson, the lightest of the hadrons and the most copiously 
produced particle in a high energy heavy ion collision. 

10.1 Raw Yields 

The starting point to obtain the raw yields for the identified particle analyses is Fig. |6.2| . 
One can obtain the raw yields of pions, kaons and protons (and their antiparticles) selecting 
events according to their centrality and fitting the ^ distribution in a given y-p± phase 
space cell. The procedure adopted here is a variation from this which offers a few advantages. 
Tracks were selected according to quality criteria based on number of points on the track 
and on the pointing accuracy to the primary event vertex. The cuts used in the pion analysis 
were more stringent than the ones used for the h~ analysis. This was done mainly to select 
tracks with a good ^ resolution as it depends mainly on the number of valid ionization 
samples, see Fig.|6l^. 



For this analysis, we followed a procedure similar to the one outlined in Ref. [[165| for 
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ionization measurements. We use the truncated mean to estimate the ionization of a given 
track (Im)- The expected ionization for a given momentum and a given mass hypothesis 
(I/i)is known to an accuracy of better than 0.1%, as represented by the curves in Fig. |6.2| . 
We then construct the z variable defined as 

\n{^\ (10.11 



h 

The z variable follows a Gaussian distribution with mean zero for a given particle population 
and with standard deviation given by the measured ^ resolution as a function of the number 
of H samples, Fig.|6]|. The mean will be zero independent of momentum for the given 
particle species under consideration, since the momentum dependence is contained in Ih- 
By incorporating Ih and the measured resolution cfdE/dx-, we can obtain a new distribution, 
Ztt, that has unit width and zero mean for the particle of interest: Z.,^ = ZT^/adE/dx- 



Fig. |10.1| shows the result of fitting to the Z distributions in different phase-space bins. 
Shown is the region 0.2 < y < 0.3 for 6 different p_\_ bins in steps of 50 MeV/c starting from 
p^= 0.3 GeV/c. The abscissa is in log scale, and the ordinate is in linear scale (note that 

is already the log of the ionization). Strictly, only the pion population should follow a 
Gaussian in the z-,^ variable. One can of course construct a z variable for the other species 
{zkj Zp and Ze) in order to extract their yields. The analysis of the other particle species 
is being carried out in STAR also. Here we focus on the pions which are the bulk of the 
produced particles as can be seen from the relative heights of the different Gaussians in the 
plot. The figure shows the result of the different fits using a Gaussian for each of the particle 
species. The pion yield is extracted from the Gaussian centered at zero. In the lower right 
panel, the electron yield is small enough and close enough to the kaon population that one 
can use a single Gaussian function. The fit also stops at = 5 in the bottom right panel, 
thus the anti-proton Gaussian is not shown. 
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Jt Yield 

y = a.2- 0.3 




Figure 10.1: Fits to the Z^^ distribution to extract the yield of tt . Shown is the region 
0.2 < y < 0.3 and 0.3 < p_l < 0.6 GeV/c. 



10.2 Corrections 

The corrections apphed to the raw data follow the same procedure as for the h~ distribution 
with a few important differences. First of all, we chose to do them in cells of {y,p±) instead 
of {ri,p±). The acceptance and efficiency use the embedding procedure as in the h~ analysis. 
Embedding was done in the region .1 < < 2 GeV/c, and |y|<l. For the lowest p± bin, 
50 < p_L < 100 MeV/c, we used a full HIJING simulation. The efficiency obtained from 



embedding is shown in Fig. 10.2. The stringent requirement on 24 or more fit points on the 



track effectively removed the need to undertake a split track correction. For the momentum 
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Figure 10.2: The track reconstruction efficiency as a function of pj_ for vr" mesons. The 
efficiency is lower than in Fig. 8.7 as the requirements for tracks in the ^ analysis are more 



stringent than for the h analysis. 

range p± < 0.6 GeV/c the distortion of the spectrum due to resolution effects was found to 
be negligible. This fact and the good statistics of the sample allowed for a much finer bin 
size in the p± spectrum. This was also desirable since this narrowed the phase space where 
the fits were made to obtain the yield. Event though the use of the z variable fixes the region 
where the pions are found, the other particles will still vary with momentum. Therefore, 
a narrower p± bin will reduce this effect making the fits more stable. The final spectrum 
was also corrected for the products of weak decays and for the expected contribution of 
pions resulting from secondary interactions. Since for this analysis the requirement on 
the distance of closest approach to the primary interaction vertex was set at 1 cm, the 
background correction was smaller than for the h~ analysis, and was found to be < 5%. We 
present now the resulting distributions for identified tt~ mesons. 

10.3 p± and m_L Distributions 
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10.3.1 Results 



We first look at the spectra for central events. Figure |10.3| shows the p_L distribution for ir" 
for the 5% most central collisions (selected by a cut on raw h~ multiplicity). The data are 
grouped for different bins in rapidity, and scaled by successive powers of 2 for the figure. 
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Figure 10.3: vr pj_ distribution for central collisions. The different bins in rapidity are 
scaled by successive powers of 2 for display purpose. 



The more common way to show such a distribution is in transverse mass, m_i, and this 



is shown in Figure 10.4 . The main corrections applied to the raw data are acceptance, 
efficiency and contributions from weak decay background and secondary interactions. The 
effect of energy loss in the determination of the 7r~ momenta was found to be less than 2% at 
100 MeV/c, decreasing rapidly for larger momenta. We also studied the effect of momentum 
resolution as for the h~ data and found that the main effect was the feeding of the lowest 
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Figure 10.4: vr (PN/{2it dy m±dm±) distribution for central collisions. The rapidity bins 
are scaled by successive powers of 2 for display purpose. 



p± bin from its nearest neighbour. The change in the slope for p± > 100 MeV/c was found 
to be <C 1%. The data for the central rapidity bins, \y\ < 0.4, was obtained in the range 
0.05 < p± < 0.75 GeV/c, which accounts for ~ 85% of the total yield. As we go towards 
the forward and backward regions, the Pj_ range where the identification via in the TPC 
can be done shrinks, and we can only measure in the range 0.2 < p± < 0.6 GeV/c. 
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10.3.2 Discussion 

Pions are the bulk of the produced particles, they are the lightest of the hadrons and have 
large interaction cross sections in nuclear matter; they are therefore expected to thermalize 
easily. In addition, their spectra are the least affected by a given collective transverse flow 
velocity, and are thus good probes to study the kinetic freeze-out properties of the system. 
The pion data are fit with a Bose-Einstein distribution of the form A/(exp(m_L/T) — 1). 
Typically, these distributions are fit with a Boltzmann distribution [|120|| . The data at low-p^ 
fit with a Boltzmann distribution always seemed to show an enhancement. This partly was 
attributed to contributions coming from resonance decays. There are also enhancements in 
the high-p^ region which can arise from contributions of pions coming from high momentum 
transfer of partons early in the collision. A consequence of all these structures was that simple 
Boltzmann fits to the same spectrum would lead to different slope parameters depending 
upon the p± region included in the fit. Since we focus here on the low p± part of the spectrum, 
and since pions are bosons, we chose to fit the data with a Bose-Einstein distribution. We 
find that this naturally gives an enhancement in the very low p± region of the spectrum. 
Also, we have already estimated the contribution from weak decays to the distribution and 
corrected for this effect. We have not estimated the contribution from resonance decays 
however. We fit the data over the measured range, instead of excluding the low momentum 
data from the fit, and find no significant enhancement at low p±. The absence of a noticeable 
contribution from resonances might be due to reinteraction of the daughter decay particles. 
A fit using a Boltzmann distribution for the range above p± = 0.3 GeV/c causes a variation 
in the slope parameters of ~ 10 MeV. Since the data at forward and backward rapidity 
have a smaller range in p± as a systematic check we performed a fit restricted to the pj_ 
range which was measured for all the rapidity bins. This resulted in a change in the slope 
parameters of 5-10 MeV for the \y\ < 0.4 spectra. We find the slope parameter for the 
central collisions to be Tgfj = 210 ± 20 MeV, where the systematic uncertainty comes from 
studying the variation in the slope parameter with several changes, including variations in 
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the fit range, choice of fit function, centrahty selection, and analyzing the pions in r] instead 
of y. The error on Tgfj from any particular fit is on the order of 5 MeV. The rapidity 
dependence of the vr" yield and of the slope parameters is discussed in the following section. 

The shape of the low p± pion spectra, in contrast to the h~ distribution, shows very 
little change when going from SPS energies to RHIC. The slope parameters reported by 
NA49 [p^ for vr" near midrapidity are Tes = 188 ± 6 MeV for 7r+ and Tes = 192 ± 3 MeV 



using an exponential fit. WA98 data for vr'' spectra for central collisions yields also similar 
values of Tcff ~ 204 (although fitting the data with a Bessel function, m±Ki{m±/Tcs) yields 
a smaller value of Tofj = 155, so one must be careful when comparing). The pion spectra from 



NA45 1 104] were analyzed in terms of a local inverse slope changing with mj_ to take into 
account the difference at low and high p±. For the low mj_ range (m_L — < 0.8 GeV/c^) 
the local inverse slope is found to be in the range 180 — 200 MeV. The slopes at the 



AGS have smaller values, E895 finds |115| T^s ~ 110 MeV. While the initial temperature 
reached at RHIC is expected to be larger than at the SPS and the AGS, the vr" spectrum 
probes the final state. The saturation of the vr" slope parameter would indicate that the 
freeze-out temperature, assuming that thermal equilibrium is reached, is the same at SPS 
and RHIC. This does not mean that the final state is the same, for there is evidence of an 
increased collective radial flow velocity at RHIC based on the preliminary p± distributions 



of heavier particles from STAR [|45| , |166| . Since the collective velocity affects the heavier 
particles the most, the pion inverse slope is the one numerically closer to the actual freeze- 
out temperature, Tf.o. • The slopes of the and p will be much larger than for the vr" even 
if they all freeze out at the same temperature in the presence of large radial flow. 

In order to gain further insight into the dynamics of the system, a study of the heavier 
particle spectra can be made to address the question of radial flow. In a different approach. 



the combination of the 2-particle correlation results [tl67|| along with the pion spectra is also 
useful: they determine the 6-dimensional pion phase space density, i.e. the dimensionless 
average number of pions per 6-dimensional phase space cell (see e.g. [168[). Such analyses 



can help to disentangle the intrinsic freeze-out temperature and the transverse flow velocity. 
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For example, for the 5% most central data, preliminary analysis of the phase space density 
yield values of Tf o. ~ 120 MeV and ~ 0.53. This is a rather large flow velocity, yet 
since the pions are so light, the slope parameter Tefj is only increased to ~ 200 MeV, while 
the slope parameter for the heavier K~ increases to ~ 300 MeV and is even larger for 
the p The systematics of the centrality dependence of the pion phase space density 

is a currently active area of study in STAR which will certainly help to shed light on the 
determination of the kinetic freeze-out conditions of the system. 

10.4 Rapidity Distribution 
10.4.1 Results 



By integrating the p± spectra in Figure 10.3 for each of the rapidity bins, we are able to 



obtain the dN/dy vs. y distribution for vr". The result is shown in Figure 10.5 . Several 
data points are shown. The black circles are obtained by integrating the p± spectrum in the 
limited p± region 0.2 — 0.6 GeV/c where we have data for \y\ < 0.8. This shows that the 
measured yield in this region is relatively independent of rapidity. The measured yield in 
this region is, however, only about half of the total yield. To illustrate the behaviour of the 
yield in the full available p± region (0.05 — 0.75 GeV/c) for a narrower y region, we plot the 
data shown as the black squares. Finally, the open circles are the extrapolated yield to all 
p± based on the fits to the data using the Bose-Einstein distribution. We find the pion yield 
to be dN/dy\\y\^Q i = 286 ± 10 for the 5% most central events. Not only is it important to 
study the yields as a function of rapidity, but it is equally important to examine the rapidity 
dependence of the slope parameters. At the BPS and at the AGS, they have been found 
to decrease with increasing \y\. The Tefj slope parameter extracted from the Bose-Einstein 
fit to the data is shown in Figure M . The error bars in the figure are the uncorrelated 



point-to-point systematic uncertainties on Tefj, the overall correlated systematic uncertainty 
is ±20 MeV. The rapidity dependence of the yield, the slope parameter Tgfj and {p±) are 



collected in Table 10. 1|. The p± distributions for the most central data for the various y bins 
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Figure 10.5: vr" rapidity distribution. The black circles (squares) are obtained summing the 
yield between p± 0.2 - 0.6 GeV/c (0.05 - 0.75) from Fig. 10.3 . The hollow data points are 
the yields obtained using the fit function extrapolated to all p±. 



are collected in Table 10.3. 



10.4.2 Discussion 

The rapidity distribution has nearly a plateau shape at mid-rapidity. However, there re- 
mains a slight rapidity dependence of the slope parameter Tgff . This behaviour is indicative 
that the idealized boost-invariant mid-rapidity region is not yet reached in Au -|- Au col- 
lisions at y^s^= 130 GeV. A possible cause for the observed effect can be a change in 
the baryon content of the system with rapidity. In the presence of radial flow driven by 
the dominant vr" mesons, a larger fraction of baryons would take away some of the pion's 
kinetic energy, reducing the observed slope parameter. The anti-proton to proton ratios 
reported by STAR |161] and BRAHMS ||162| ] show a drop in the p/p ratio from 0.65 at 
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Figure 10.6: vr slope parameter for the different rapidity bins from Fig. 10.4. The flat line 
at Tgfj = 210 is the average for \y\ < 0.1 



mid-rapidity to 0.41 for the region y ~ 2, which already hints at such effects, but a better 
understanding of the dynamics must wait for the measured baryon yields and the baryon 
rapidity distributions. 

10.5 Centrality Dependence 
10.5.1 Centrality Selection 

For the vr" centrality dependence, we chose to do a selection based on the event multiplicity. 
This is in contrast to the centrality selection for the h~ analysis which was based on ZDC 
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Table 10.1: tt dN/dy, T^q and ?;s. y. 

and CTB trigger signals. The reason for this choice was simply to use a common set of 
cuts which would be useful not just for this analysis, but also for other studies such as 
particle correlations (Hanbury-Brown Twiss) and in particular for the studies of the pion 
phase space density. Having a standard set of values based on the (uncorrected) h~ yield 
facilitates comparison and combination of different observables. There was a worry that since 
the pions are the main component of the charged hadron spectra, a centrality selection based 
on charged multiplicity would introduce a sizable auto-correlation. Certainly, the shape of 



the h~ multiplicity distribution for 5% most central events (the shaded region in Fig. |9.lD 
would look different by selecting events with a straight cut on uncorrected multiplicity. The 
mean of this distribution, i.e. the value of dN/drj for the region |?7| < 0.5, will also shift 
to higher values by doing such a selection. The question is how significant is this shift. 
We studied this by selecting the 5% most central events based on uncorrected multiplicity 
and comparing the mean of the resulting distribution to the mean of the histogram in the 
shaded region of Figure |9.1[ The mean of the histogram using the centrality definition based 
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on multiplicity was 4.5% higher than than the one using the ZDC as centrality definition. 
The difference in the mean of the histograms using these two centrality definitions will also 
decrease for more peripheral collisions. We therefore conclude that the differences are not 
problematic, and when encountered, they can be reconciled by at most a 4.5% effect. 



10.5.2 mj_ Distributions vs. Multiplicity 



Figure |10.7| shows the m_L distribution for 10 different centrality selections. The data are 
taken in the rapidity interval < y < 0.1. Again, the data are fit by a Bose- Einstein 
distribution which agrees well with the spectra even at the lowest m±. The slope parameters 
are very similar for the different centralities. They are collected and plotted vs. the h~ mean 




Figure 10.7: 7r~ m± distribution for different multiplicity selections. The distributions are 
plotted in the rapidity slice < y < 0.1 units. The data are fit to Bose- Einstein distributions 
(curves) from which we extract the slope parameter T^g. 



multiplicity of the corresponding centrality bin in Figure |10.8| . The horizontal error bars are 
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Figure 10.8: Centrality dependence of the slope parameter obtained by a fit to the Tr~ m± 
distribution. 

the systematic uncertainty in the determination of the mean multiplicity of each centrality 
bin. The vertical error bars are the uncorrelated systematic uncertainty in the measurement 
of the slope parameter. For the most peripheral bins, the statistical error is also important, 
and the error bar is therefore larger. 

We see that there is a only a slight dependence of the slope parameter with centrality. 
The main difference happens for the 2nd and 3d most peripheral bins (the most peripheral 
bin has a larger error, so it is hard to see a systematic effect). The 2nd most peripheral bin 
has a slope T^g = 177 =b 8 MeV, while the most central bin has a slope T^s = 210 it 4 MeV 
(where the errors are the uncorrelated systematic uncertainties for the comparison). This 
yields an increase of 19 it 5%. The increase occurs rapidly, as the 60% of (TAuAu bin already 
has a slope parameter of Teg = 196ib9MeV, and from then on the data are consistent with 
having no further centrality dependence. For comparison, the increase in h'^ {p±) from the 
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70 — 80% bin to the most central bin is 14ib3%. In addition, the h~ distribution includes a p 
component, and the p slope has a stronger centrality dependence than the vr" slope. At SPS 
energies, the centrality dependence of the pion slopes was also found to be rather constant. 
This finding is also in agreement with the argument that the slope of the vr" spectrum is 
not changed significantly by the collective expansion. 
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Table 10.2: vr" S N j {li^m ^dm x_dy) at y = 0.05 for various h multiplicity bins. Fig. |10.7| 
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0.375 


1040 ± 36 


1045 ± 30 


991 ± 31 


1047 ± 29 


990 ± 26 


990 ± 24 


1002 ± 31 


964 ± 26 


0.425 


798 ± 35 


839 ± 31 


784 ± 29 


842 ± 24 


749 ± 22 


714 ± 21 


777 ± 19 


737 ± 18 


0.475 


660 ± 20 


657 ± 20 


613 ± 21 


677 ± 21 


594 ± 18 


589 ± 16 


605 ± 15 


574 ± 15 


0.525 


500 ± 18 


494 ± 14 


484 ± 15 


517 ± 16 


457 ± 16 


454 ± 22 


457 ± 14 


452 ± 37 


0.575 


402 ± 12 


364 ± 13 


389 ± 14 


404 ± 13 


357 ± 22 


368 ± 13 


331 ± 14 


328 ± 24 


0.625 


299 ± 10 


299 ± 10 


315 ± 10 


299 ± 10 


273 ± 11 


294 ± 11 


262 ± 14 




0.675 


229 ± 8 


227 ± 9 


231 ± 9 


231 ± 10 


221 ± 11 


224 ± 10 






0.725 


174 ± 8 


171 ± 16 


174 ± 22 


163 ± 7 


166 ± 7 


179 ± 9 







Table 10.3: tt (PN/ {p±dp±dy) for the 5% most central events. 
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Conclusions and Summary 



We have studied the general characteristics of charged particle production in Au + Au colli- 
sions at -^/i = 130 GeV. The gross features of the cross section indicate that the collision is 
dominated by geometry, as evidenced by the shape of the multiplicity distribution. Glauber 
model studies are able to reproduce the measured da/dNf^- distributions and are useful to 
obtain a statistical determination of the number of participants in the collision for a given 
event sample, a quantity related to the impact parameter of the collision. This determi- 
nation of the collision centrality will be used in future studies of the centrality dependence 
of many experimental observables, such as the production of strange mesons and baryons, 
charge fluctuations and identical particle correlations to name a few. 

Prom the yield of charged hadrons, we find that particle production per participant in 
central collisions increases by 38% relative to pp and 52% compared to nuclear collisions at 
V^NN = 1^ GeV. The h~ p± spectrum in Au -|- Au distribution is harder than that of the pp 
reference system at similar centre-of-mass energies for the p± region up to 2 GeV/c. Scaling 
of produced particle yield with number of participants shows a strong dependence onp±, with 
Wounded Nucleon scaling achieved only at the lowest measured p±_ . The h~ pseudorapidity 
distribution is almost constant within \r]\ < 1, indicating that at these energies we are 
approaching a boost invariant mid-rapidity region. This finding contrasts the results from 
lower energies where the 77 distribution was found to peak at mid-rapidity. The shape 
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of the ri distribution in the mid-rapidity region is found to be similar for all centralities. 
Comparing the p± distributions for different centralities to the distribution measured in 
peripheral collisions, we find a rise in the number of produced particles with increasing p± 
up to p± < 2 GeV/c. This is also consistent with the observed rise when comparing to pp 
collisions. However, the ratio of Au + Au to pp indicates that the limit of scaling with the 
number of binary collisions A'coii is not reached. As for future directions with p± spectra, 
an interesting indication coming from preliminary data at high p± is that there a turnover 
in the shape of this ratio beyond 2 GeV/c, which has been intensely debated in the context 
of being a scenario consistent with partonic energy loss in a QGP. 

The measured h" {p±) increases by 14% from peripheral (70 — 80%) to central collisions. 
The shape is found to be relatively smooth. If deconfinement is reached in these collisions, 
this behaviour is consistent with a cross over region between the deconfined and hadronic 
phase, instead of a strong first order phase transition. Prom the combination of the pseudo- 
rapidity and p± distribution results, we discussed the applicability of the Bjorken estimate 
of the initial energy density. The lower limit of e > 3.9 GeV/fm^, taken at face value and 
comparing to fi = lattice results, leads to the conclusion that necessary conditions for QGP 
formation are in fact reached. Stated another way, to the question of whether the QGP can 
be formed in the high energy RHIC collisions, the answer is affirmative. Of course, this 
does not mean that one has a proof of deconfinement, but it helps to set the stage for the 
measurement of the additional probes of QGP formation, which are needed to provide the 
detailed evidence. 

The rapidity distribution of negative pions is found to be relatively flat in the range 
|y| < 0.8. Although this is still a small range compared to the rapidity gap of '^lO units 
between the colliding nuclei, this again signals the approach to the limit of a boost-invariant 
baryon-free mid-rapidity region. The pion transverse mass spectra below 0.6 GeV/c are well 
parameterized by a Bose-Einstein distribution with inverse slope parameter T ~ 190 MeV. 
The slope parameter is relatively independent of centrality, perhaps growing with increasing 
centrality up to a centrality of 60% of uauAu but then saturating. For the most central 
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collision, the slope parameter shows a small but systematic drop with rapidity. This is an 
indication that boost invariance is not fully achieved. A possible cause for this is increase 
in the baryon content with rapidity. The rapidity and transverse momentum distributions 
of protons and anti-protons will help shed more light on this subject. The measurement of 
the m± spectra of more massive hadrons (kaons, protons, As, etc.), as well as a combination 
of identical particle correlations and spectra, are natural extensions to the studies presented 
here in order to determine the conditions at kinetic freeze-out such as transverse flow velocity 
/35°^ and the kinetic freeze-out temperature Tf.o.- 

Identifying deconfinement and a complete characterization of the QGP will require data 
from a wide spectrum of experimental observables, systematic variations of centre-of-mass 
energies and nuclear species, pp and pA reference data and candid guidance from theory. 
We hope that the results presented here provide useful baseline information to both the 
experimental and theoretical efforts in this search. 
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